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Summary 

Helicopters have historically had significant vibration problems that require stringent design speci- 
fications for helicopter subsystems. The primary source of the vibratory loads may be traced to the 
main-rotor system where the rotor blades are excited by the periodic aerodynamic environment. The 
resulting vibratory loads are transmitted down the main-rotor shaft into the helicopter fuselage. Gener- 
ally the effects of the rotor coupling with the fuselage response are ignored in rotor system analysis. 
However, these effects are often important and may significantly alter the helicopter system response. 
At the Langley Research Center, an active mount rotorcraft testbed is being developed for use in the 
Langley Transonic Dynamics Tunnel. This testbed, the second generation version of the Aeroelastic 
Rotor Experimental System (ARES-II), can impose rotor hub motions and measure the response so that 
rotor-body coupling phenomena may be investigated. An analytical method for coupling an aeroelasti- 
cally scaled model rotor system to the ARES-II is developed in the current study. Models of the testbed 
and the rotor system are developed in independent analyses, and an impedance-matching approach is 
used to couple the rotor system to the testbed. The development of the analytical models and the cou- 
pling method is examined, and individual and coupled results are presented for the testbed and rotor 
system. Coupled results are presented with and without applied hub motion, and system loads and dis- 
placements are examined. The results show that a closed-loop control system is necessary to achieve 
desired hub motions, that proper modeling requires including the loads at the rotor hub and rotor control 
system, and that the strain-gauge balance placed in the rotating system of the ARES-II provides the best 
loads results. 
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1. Introduction 


1.1. Motivation 

Aeroelastic testing of model-scale helicopter rotors has been performed in the Langley Transonic Dynamics 
Tunnel (TDT) since the 1960s. Since 1976, the bulk of this testing has been accomplished by using the Aeroelastic 
Rotor Experimental System (ARES), a generic rotorcraft testbed designed for the performance, stability, and 
dynamic loads evaluation of new rotor concepts. Currently, a second generation version of the testbed is under 
development at the TDT. This testbed is called ARES-II (fig. 1.1) and has been designed for the evaluation of 
dynamic loads due to rotor-body coupling and for extended rotor stability testing. To comprehend fully the 
expected capability of the ARES-II, new analysis techniques must be established. These techniques will be used to 
develop pretest predictions, to formulate efficient test plans, and to develop safety-of-flight criteria for safe, 
practical wind tunnel testing. The goal of the current research is to develop analytical techniques for the prediction 
of rotor-body coupling effects. Because the ARES-II may be used for testing of many model rotor system designs, 
developing a coupling method that permits rotor-off treatment of the testbed for evaluation and verification is 




Figure 1.1. ARES-II rotor research testbed. 



necessary. Therefore, a method is employed which permits the development of independent math models of the 
testbed and the rotor system before the rotor-to-body coupling is accomplished. 

1.2. Background 

Helicopters have historically had significant vibration problems. These vibrations can, at times, be severe and 
are responsible for a great number of stringent design requirements for helicopter subsystems. Vibratory loads 
cause crew and passenger fatigue and discomfort, reduce overall system reliability, and increase maintenance 
costs. For military helicopters, vibratory response can make targeting operations difficult. Further, vibration levels 
must be attenuated for the helicopter to be accepted by the public as a viable means of civil transportation. 
Achievements in vibration reduction, however, are quickly being offset by the requirement for higher flight 
speeds. Reductions in flight vibrations, therefore, are becoming more and more difficult to attain. For these rea- 
sons, helicopter vibration reduction has been, and continues to be, a primary area of rotorcraft research in industry, 
academia, and government organizations throughout the world. 

In 1976, a NASA advisory council challenged the helicopter research community with the goal of “the jet 
smooth helicopter,” a maximum vibratory loading of 0.02 g throughout the straight-and-level flight envelope. To 
date, no production helicopter has achieved this goal. Advances in vibration control, however, have been signifi- 
cant. Vibration levels had dropped from 0.3g-0.6g in the mid-1950s to about O.lg in the early 1980s (ref. 1). The 
methods used to attain these reduced levels of vibration have been documented by Reichert in a 1981 paper sur- 
veying the developed techniques (ref. 1). Although somewhat dated, the paper presents the bulk of vibration 
reduction techniques in use on production aircraft today. Furthermore, little additional progress has been made 
toward achieving the jet smooth goal since Reichert's paper was published. 

It is widely accepted that the primary source of helicopter vibratory loads may be traced to the main-rotor sys- 
tem. Other loads exist (e.g., fuselage airloads, rotor-fuselage aerodynamic interference, drive train loads, tail rotor 
loads); however, the effects of these loads are typically small when compared with those produced by the main 
rotor that are passed down the main-rotor shaft into the fuselage. Therefore, the traditional methods of helicopter 
vibration reduction have usually consisted of means to reduce the vibratory loads generated by the main rotor. 

Traditionally, the method used to determine vibration levels throughout the fuselage has been to analyze the 
rotor as an isolated system when calculating the loads imposed on a fixed rotor shaft. These loads were then 
applied independently to an analytical model of the fuselage to determine vibration levels at locations throughout 
the helicopter. However, rotor-body couplings have been known for some time to exist between rotor systems and 
the fuselage that can alter the loads. Several analytical techniques have been implemented to couple the rotor sys- 
tem to the fuselage. One such technique is known as “impedance matching,” where the load and displacement 
fields at the interface between the rotor system and the fuselage are matched. References 2 through 6 are represen- 
tative of rotor-body coupling research using impedance-matching methods. An advantage of the impedance- 
matching method is that the rotor and body may be treated by using different independent analyses. Therefore, one 
is free to choose the analyses that can best model the individual systems. 

1.3. Scope 

The scope of the current research includes the development of an impedance-matching method suitable for 
coupling a generic research rotor to the ARES-II rotorcraft testbed. The ARES-II testbed is capable of six-degree- 
of-freedom actuation of the rotor hub through a motion-based system called a Stewart platform. Therefore, this 
system and its control must be modeled as well. A three-step approach is used for the analysis. First, independent 
analytical models of the ARES-II and the Basic Research Rotor (BRR), the rotor system chosen for this study, are 
developed (sections 4.1 and 4.2). Second, the developed models of the ARES-II and the BRR are used to generate 
a set of impedance matrices representing the load-displacement relationships at the interface joints on each model 
(sections 5.3.1 and 5.3.2). These interface joints are the points at which the two systems are joined during the 
coupling phase of the analysis. Finally, the two systems are coupled together by enforcing displacement and load 
compatibility conditions through simple linear algebraic manipulations of the impedance matrices (section 5.3.3). 
A survey of the available analyses was undertaken to decide which would provide the best independent models of 
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the ARES-II and the BRR, the rotor system chosen for this study. The most recently available, state-of-the-art 
analysis for rotor systems is CAMRAD-II, the second generation version of the Comprehensive Analytical Model 
of Rotorcraft Aerodynamics and Dynamics (refs. 7 and 8). This program was chosen for the Basic Research Rotor 
modeling. To model the ARES-II, the Dynamic Analysis and Design System (DADS), a multibody dynamics 
analysis, was chosen (refs. 9 and 10). Alternatively, a finite-element code could have been used. The ARES-II, 
however, is a relatively stiff system with discrete areas of flexibility — a system a multibody dynamics code is well 
suited to model. Additionally, DADS offers a superior interface for visualizing the testbed motion and for imple- 
menting the control of the motion-based platform. 

The impedance- matching methods proposed in references 2 through 6 have been used for guidance in this 
study. However, two enhancements have been included in the current effort that are not evident in the literature. 
First, the controlled motion-based platform of the ARES-II is unique among rotor testbeds. For this work, the 
impedance-matching method has been extended to include the motions introduced by the platform. Second, two 
load paths have been modeled between the testbed and the rotor system in the current work — the hub and the rotor 
control system. Only the hub load path is evident throughout the literature. The rotor control system load path is 
shown to be necessary to provide an accurate assessment of the system loads. 

1.4. Organization of Study 

The rest of the study is organized in the following manner: 

Section 2. Basic Rotorcraft Analysis: This section provides a brief introduction to rotorcraft analysis for read- 
ers unfamiliar with such techniques. Basic rotor system definitions, control, performance, aerodynamics, and 
dynamics are discussed. 

Section 3. Rotorcraft Wind Tunnel Test Apparatus: This section provides a complete description of the 
ARES-II testbed and the BRR. 

Section 4. Analytical Models of ARES-II and BRR: This section describes the DADS model of the ARES-II 
and the CAMRAD-II model of the BRR. Also described is the development of the equations necessary to pre- 
scribe motion of the ARES-II Stewart platform. 

Section 5. Dynamic Analysis Using Impedance Matching: This section describes the application and limita- 
tions of the impedance-matching approach as applied to the ARES-II and the BRR. A simple spring-mass-damper 
system is used to convey the proper impedance-matching techniques. A complete description of the steps neces- 
sary to develop the impedance models and coupling equations for the ARES-II and the BRR are included. 

Section 6. Presentation of Results: The results obtained with the individual models of the ARES-II and the 
BRR are discussed in this section. Also included are coupled rotor-testbed results for actuated and unactuated 
motion of the ARES-II Stewart platform. 

2. Basic Rotorcraft Analysis 

This section provides a brief introduction to rotorcraft analysis for readers who are unfamiliar with such tech- 
niques. Several good texts (refs. 1 1 through 14) are available which discuss rotor system aerodynamics, dynamics, 
and stability in detail. This discussion emphasizes the major analytical techniques used in analyzing the Basic 
Research Rotor (BRR) system. 

2.1. Definitions 

2.1.1. Azimuth 

For most helicopters, the main-rotor system rotates in a counterclockwise direction (when viewed from above) 
at an angular rate of £2 in radians per second. Azimuthal blade position is denoted by \|/ and is indexed so that 
V|/ = 0° occurs when the blade chosen to be the reference blade is over the aircraft tail and increases in the direction 
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of rotor rotation. Typically, rotor blades are equally spaced azimuthally so that a general relationship exists 
between each blade and the reference blade; that is 


V fc = V + -^- (2.1) 

where 6 is the blade number, N is the number of rotor blades, \| t is the azimuth of the reference blade in radians, 
and \\i b is the azimuth of the 6th blade in radians. Therefore, by using equation (2.1), the reference blade is the 6th 
blade. 

2.1.2. Rotor Speed Parameters 

Several important parameters are associated with the rotational speed of the rotor system. The rotor tip speed 
in hover is defined by 


V tip = OR 


( 2 . 2 ) 


where R is the rotor radius in feet. An associated parameter is the rotor hover tip Mach number given by the fol- 
lowing equation: 



(2.3) 


where a is the test medium speed of sound in feet per second. The hover tip Mach number is typically used to 
ensure proper aerodynamic scaling when testing rotor systems for performance characteristics in the Langley 
Transonic Dynamics Tunnel. 

Rotor advance ratio is typically used to refer to the forward flight speed of a helicopter and is defined by 

V M cos ctj-pp 

H = r; (2.4) 

tip 

where V", is the forward flight speed and a TPP is the rotor tip-path-plane angle of attack with respect to V^. A 
well-designed rotor system with good performance characteristics may be able to reach fi = 0.45. A poor perform- 
ing rotor system may have difficulty in attaining (i = 0.35. 


2.1.3. Articulated Rotor Hubs 

An articulated rotor hub has three hinges for each blade: the flap hinge, the lead-lag hinge, and the pitch bear- 
ing, each of which is shown in figure 2.1 for the rotor hub used with the BRR. The lead-lag and flap hinges permit 
blade motions in the plane and out of the plane of rotation, respectively. The pitch bearing allows each blade to 
change angle of attack under the influence of the rotor control system. 


2.1.4. Coordinate Systems 

Referring to the “fixed system” and the “rotating system” is typical in rotorcraft analysis. Generally, the fixed 
system consists of the helicopter fuselage and its components. The rotating system consists of the rotor and all 
components operating at the same rotational speed as the rotor system such as the main-rotor shaft and rotating 
control system components. For the purposes of this study, however, the rotating system is considered to include 
only those rotor system components that are part of the rotor hub or the blades. All other rotating-system compo- 
nents, including the rotor control system and rotor shaft, are considered to be part of the fixed system. The reasons 
for this consideration become evident in section 4.2.2. 

The coordinate system referred to throughout this report is a typical rotorcraft fixed-system coordinate frame 
oriented as shown in figure 2.2. Unless otherwise specified, all motions and loads are identified with respect to this 
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Lead-lag axis 


Figure 2.1. Basic Research Rotor articulated hub. 


Up 



Aft 


Figure 2.2. Fixed-system coordinate frame. 


coordinate frame. For example, a “positive axial force” is a load oriented along the +X-axis, and a “negative pitch- 
ing moment” is a moment defined by using the right-hand rule along the -K-axis. 

2.2. Rotor Control System 

The rotor control system consists of the swashplate, swashplate actuators, and the pitch links. A close-up view 
of the control system of the ARES-II with the BRR articulated hub is shown in figure 2.3. The swashplate has both 
rotating and stationary sections. The rotating section turns with the hub and serves as the lower attachment point 
for the pitch links. The upper end of each pitch link is attached to the pitch hom which is part of the blade cuff. 
Both ends of the pitch links have spherical rod-end bearings to permit rotational motions at the joints. The 
stationary section of the swashplate is mounted to the shaft with a large spherical bearing that is free to travel up 
and down the rotor shaft. This arrangement eliminates two degrees of freedom from the swashplate (fore and aft 
and side to side), and the nonrotating scissors (fig. 2.3) eliminate a third degree of freedom (rotation about the 
shaft). The remaining three degrees of freedom are constrained by three swashplate actuators. The actuators con- 
trol the swashplate position up and down the rotor shaft and the swashplate tilt from side to side and fore and aft. 
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Pitch horn 



Figure 2.3. ARES-II control system with BRR hub. 


For control systems with the pitch links trailing the blades, a lowering of the swashplate increases the pitch on all 
the blades collectively. Thus, raising and lowering the swashplate controls the “collective” pitch of the blades. By 
tilting the swashplate, a “cyclic” blade pitch schedule may be induced; this means that a sinusoidal oscillation of 
blade pitch is induced as a function of azimuth. Because the swashplate has two axes of tilt (i.e., fore and aft and 
side to side), two cyclic controls are necessary. The cyclic control that induces forward or aft motion of the heli- 
copter is called longitudinal cyclic, and the control that induces right or left sideward motion of the helicopter is 
called lateral cyclic. 

2.3. Rotor Performance 

Because this study focuses on rotor dynamics, little needs to be said regarding rotor performance; however, 
some explanation is necessary because rotor performance parameters are used to define trim targets in section 5. 
Generally, three parameters are of primary concern for rotor performance: rotor-generated lift, propulsive force, 
and torque required. Each of these parameters is normally nondimensionalized so that comparisons may be made 
between rotor systems regardless of size. However, the typical nondimensionalization factors for fixed-wing air- 
craft fail when applied to rotary- wing aircraft in hover flight (V«, = 0). The accepted equations for nondimension- 
alization of rotor systems allow for treatment of both the hover and forward flight conditions as follows: 


C L ~ 


Jt/? 2 p(U/?) 2 


(2.5) 
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D 


( 2 . 6 ) 


nR 2 p{QR) 2 

C Q = , Q — l < 2J > 

nR p(QR) 

where L is lift, D is drag (or propulsive force), Q is torque, p is density, R is rotor radius, and Q is rotor rotational 
speed. Notice that drag D is used for propulsive force. For a wind axis system, lift is normally considered positive 
up and drag positive aft. Therefore, referring to the propulsive force (which is oriented forward) as “negative drag” 
is typical. 

2.4. Rotor Aerodynamics 

The rotor system develops its thrust by inducing a mass flow, or induced velocity, downward through the 
rotor. The total flow through the rotor, or inflow, is a combination of the flight and induced velocities. This inflow 
is a highly nonuniform and volatile flow field that is difficult to model and to analyze. The inflow is, however, a 
critical characteristic in the determination of rotor aerodynamic loads because the local inflow directly affects the 
local blade section angle of attack. Further, the rotor system aerodynamics represents the primary forcing function 
responsible for the dynamic response. Therefore, a good model of the rotor system aerodynamics, and hence, a 
good model of the inflow, is necessary. 

The rotor inflow is influenced by the flight velocity and the rotor wake. Therefore, a representation of the rotor 
wake structure (blade tip vortices and inboard blade vortex sheets) is necessary to calculate the rotor inflow. In 
CAMRAD-D, wake models are available based on the research results of Landgrebe (ref. 15), Kocurek and 
Tangier (ref. 16), Scully (ref. 17), and Johnson (ref. 18). 

References 15 through 18 describe two approaches to model the rotor wake: the prescribed wake and the free 
wake. The simpler of the two is the prescribed wake methodology. In the prescribed wake method, a series of tra- 
jectories is calculated for the tip vortices and the inboard wake sheet based on rotor operating conditions. The 
inflow environment is calculated based on the position and the circulation strength of the wake elements at each 
azimuthal step in the analysis. An example of the wake geometry obtained with this method is shown for hover 
and forward flight in figure 2.4. The second approach is the free wake methodology in which the tip vortex trajec- 
tory is allowed to distort based on the conditions of the surrounding airflow; this results in the most realistic aero- 
dynamic environment. An example of the wake geometry obtained with this method is shown in figure 2.5. For 
this study, a wake model based on the second approach, the Johnson free wake methodology (ref. 1 8), has been 
used for all calculations to obtain the greatest possible model fidelity. 

2.5. Rotor Dynamics 

Main-rotor system dynamic response is the primary contributor to helicopter vibratory loads. The aerody- 
namic environment of the rotor system in forward flight produces periodic blade loads at frequencies that are har- 
monics of the rotor rotational speed. Therefore, these frequencies are the primary concern of the helicopter 
dynamicist. Then, referring to the frequencies as referenced to the rotor speed becomes convenient. Thus, the ter- 
minology n per revolution, n per rev, and nP is introduced. Each refers to n times the rotor rotational speed Q or 
n£l as is sometimes observed in the literature. As a simple example, consider the BRR at its nominal operating 
speed of 660 rpm or 1 1 Hz. A frequency of 22 Hz would be considered 2 per rev; a frequency of 33 Hz, 3 per rev; 
and so forth. Throughout this report, the nomenclature n per rev and nP will be used interchangeably to refer to 
frequencies. As in any structural dynamics problem, the natural frequencies of the rotor blades are of considerable 
concern. However, their calculation is considerably complicated because of the rotational environment of the rotor 
system. Figure 2.6 shows the Southwell or “fan” diagram of the BRR as calculated by CAMRAD-II. This figure 
shows rotor blade frequencies in both Hertz and per rev as a function of rotor speed in revolutions per minute.The 
lowest two frequencies are rigid-body lead-lag and flap modes representing rotation about the lead-lag and flap 
hinges. The higher frequencies are elastic blade modes that have the general appearance of pinned-ffee elastic 
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(a) Hover. 


(b) Forward flight. 


Figure 2.4. Prescribed wake geometry for hover and forward flight. Forward flight geometry shows only tip vortex trajectory. 



Figure 2.5. Free wake geometry in forward flight. Only tip vortex trajectories are shown. 


beam modes. The per-rev lines in figure 2.6 represent the frequencies at which the aerodynamic forcing functions 
excite the blade modes. 

The rotor system dynamics problem is one of forced response in which the aerodynamic loads supply the forc- 
ing functions at discrete frequencies (i.e., per rev). The next concern is the determination of which frequencies 
have the greatest effect on the helicopter vibratory loads. Appendix A contains a detailed development of the man- 
ner in which rotating-blade loads are transferred down the shaft into the fuselage or fixed-system components of 


8 





Figure 2.6. Basic Research Rotor rotating frequencies. Vertical dashed line at 660 rpm represents nominal operating speed 
of BRR. 

the aircraft. Here it is sufficient to point out that the rotor system normally acts as a loads filter for the fixed sys- 
tem, passing down only loads whose frequencies are integer multiples of the blade passage frequency (i.e., fre- 
quencies of kN per rev, where k is an integer and N is the number of blades). For a four-bladed rotor system such 
as the BRR, the blade passage frequency is 4 per rev so that only loads having frequencies of 4P, 8P, 12P, and so 
forth are passed to the fixed system. As shown in appendix A, these fixed-system loads result from loads in the 
rotor system at frequencies of kN per rev and (k ± 1 )N per rev. Therefore, for the BRR, 4-per-rev fixed-system 
loads are generated by 3P, 4P, and 5P loads in the rotor system, and 8-per-rev fixed-system loads are generated by 
7P, 8P, and 9P loads in the rotor system. Furthermore, the harmonic loads generally become smaller with increas- 
ing frequency. Therefore, considering only frequencies at 9 per rev and below is typical in rotor dynamics analysis 
of four-bladed rotor systems. This is a key point in the analysis of rotor-body coupled systems since only two 
fixed-system frequencies (4P and 8P) are considered critical for the analysis of a four-bladed rotor. 


3. Rotorcraft Wind Tunnel Test Apparatus 

The scope of the current research does not include any wind tunnel testing. However, since the objective is to 
develop analytical models of rotorcraft wind tunnel test hardware, descriptions are provided in this section to 
acquaint the reader with the apparatus necessary to perform such testing. This description gives the reader a more 
thorough appreciation of the basis for model-scale aeroelastic rotor testing in the Langley Transonic Dynamics 
Tunnel. 

3.1. Langley Transonic Dynamics Ttinnel 

The Langley Transonic Dynamics Tunnel (TDT) is a continuous-flow pressure tunnel capable of speeds up to 
Mach 1 .2 at stagnation pressures up to 1 atm. The TDT has a slotted test section that is 1 6 ft square with cropped 
comers and a cross-sectional area of 248 ft 2 . Either air or Refrigerant- 12 (R-12), a heavy gas, may be used as the 
test medium. For rotorcraft research, a test medium density of 0.006 slug/ft 3 is typically used. The TDT is particu- 
larly suited for rotorcraft testing because of three advantages associated with the heavy gas. First, the high density 
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of the test medium permits model rotor components to be heavier; thereby, structural design requirements are 
eased while dynamic similarity is maintained. Second, the low speed of sound of R-12 permits much lower rotor 
rotational speeds and forward flight velocities to match full-scale Mach numbers; thus, the centrifugal blade forces 
are eased and the drive system power requirements are reduced. Furthermore, the lower rotational speeds translate 
into a smaller control bandwidth when using active control techniques that depend upon harmonic actuation. 
Third, the high-density environment increases the Reynolds number throughout the test envelope and permits a 
more accurate simulation of the full-scale aerodynamic environment. A more detailed discussion of the advantages 
of the heavy gas test medium and further discussions on rotorcraft testing in the TDT are presented in 
reference 19. 

3.2. Aeroelastic Rotor Experimental System 

The second generation version of the Aeroelastic Rotor Experimental System (ARES-II) model-scale helicop- 
ter rotor testbed is currently under development at the TDT. The stand-mounted testbed is intended for use in 
defining rotor-body coupling phenomena by using aeroelastically scaled model rotor systems. Sections 3.2.1 
through 3.2.3 describe the testbed and its major subsystems. These subsystems include a stand and longeron com- 
ponent (section 3.2.1), an actively controlled hydraulic Stewart platform (section 3.2.2), and a rotor pylon section 
(section 3.2.3). Figures 3.1 through 3.6 provide various views of the testbed. Figure 3.1 shows the ARES-II with 
the BRR installed. The fuselage shell left-side panel has been removed to show the placement of the rotor pylon 
within the shell. The aerodynamic fairing normally installed around the stand has also been removed to show stand 
details. Figure 3.2 shows the testbed with the rotor blades and the fuselage shell removed. A close-up view of the 
upper portion of the ARES-II is shown in figure 3.3. Figures 3.4 through 3.6 show detailed views of the rotor 
pylon and the Stewart platform. For all views, the ARES-II drive belts and miscellaneous hardware have been 
removed for clarity. 

3.2.1. Stand and Longeron 

The ARES-II stand (fig. 3.2) is constructed of a steel baseplate, a steel central column, four steel support col- 
umns, a steel yoke to support the longeron, and two airfoil-shaped aluminum ribs to support an aerodynamic fair- 
ing during wind tunnel testing. The longeron (figs. 3.2 through 3.5) is constructed of three steel box beams each 
2 in. square and a 1 -in-diameter steel bar. Aluminum fittings placed on the longeron (fig. 3.5) provide support for 
the hydraulically controlled Stewart platform. Figures 3.3 and 3.4 clearly identify the longeron and its placement 
on the stand. The connection between the stand and longeron is a revolute joint oriented to allow pitching motions 
of the longeron relative to the stand. An electric motor-driven screw-jack actuator placed behind the stand is used 
to control the pitch motion of the longeron. This arrangement permits the adjustment of rotor pylon pitch (or shaft 
angle of attack), a feature necessary for proper trim of rotor systems in forward flight. 

The ARES-II stand is designed to be soft mounted to the TDT floor with elastomeric pads placed at the four 
comers of the stand baseplate. This arrangement permits limited flexibility and allows tuning of the ARES-II 
rigid-body frequencies. Alternatively, the stand may be hard mounted to the floor. 

3.2.2. Stewart Platform 

A hydraulically actuated Stewart platform is placed between the stand-longeron assembly and the rotor pylon. 
It consists of a six-degree-of-ffeedom baseplate that is driven by six hydraulic actuators attached to the longeron. 
The platform is similar in function to those commonly in use on aircraft motion-based simulators and allows six- 
degree-of-freedom actuation of the rotor pylon which is rigidly attached to the Stewart platform baseplate. 
Figure 3.4 shows the rotor pylon exploded from the longeron and Stewart platform assembly. This view shows the 
orientation of the Stewart platform actuators with respect to the longeron. Figure 3.5 shows a close-up view of the 
Stewart platform. The actuators are connected to the longeron brackets by spherical rod -end bearings. The piston 
rod of each actuator is connected to the Stewart platform baseplate with a spherical rod-end bearing. The act uato r 
rod connections provide a three-point restraint for the baseplate; thereby, six-degree-of-freedom control of the 
platform is allowed. Two of the rod-baseplate joints are labeled in figure 3.5; the third is hidden from view. 
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Figure 3.1. ARES-II with Basic Research Rotor installed. 


The ARES-II Stewart platform provides a limited range of motion about the nominal centered position of the 
platform. Translational motions are limited in range to ±0.4 in. and rotational motions to ±2.0°. For this study, the 
nominal translational and rotational rotor pylon motions applied are ±0.01 in. and ±0.05 , respectively, well 
within the range of motion of the Stewart platform. 

3.2.3. Rotor Pylon 

Figure 3.6 provides a detailed view of the rotor pylon. The rotor pylon consists of the rotor control system, the 
rotor drive system, and the force measurement system. As mentioned earlier, the rotor hub and pitch links are not 
considered part of the rotor pylon. They are shown only to provide a complete view of the assembled system. The 
rotor control system consists of three fly-by-wire hydraulic swashplate actuators and a standard rise-and-fall 
swashplate. The rotor drive system is powered by a variable-frequency synchronous motor rated at 47-hp output at 
12000 rpm. The motor is connected to the rotor shaft through a belt-driven, two-stage speed reduction system. 

Instrumentation on the rotor pylon provides continuous display of swashplate control settings, rotor speed, and 
rotor forces and moments. A 30-channel slip-ring assembly placed in the rotor shaft permits the transfer of 


11 
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Figure 3.4. Exploded view of ARES-IT rotor pylon and Stewart platform. 
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Figure 3.5. Close-up view of Stewart platform components. 

rotating-system data to the fixed system. Fixed-system forces and moments are measured by a six-component 
strain-gauge balance placed at the bottom of the rotor pylon. This balance is called the fixed-system balance. 
Rotating-system forces and moments are measured by a second six-component strain-gauge balance. This balance 
is placed just below the rotor hub between the rotor shaft and the mast and rotates with the shaft assembly. This 
balance is called the rotating balance. 

3.3. Basic Research Rotor 

For the initial investigation on the ARES-II, the BRR is the rotor system chosen. The BRR is a four-bladed, 
articulated rotor system designed for generic rotorcraft wind tunnel research, but it is not intended to represent any 
full-scale flight vehicle hardware. 

3.3.1 . Rotor Hub 

A four-bladed articulated rotor hub is used for the BRR. The hub lead-lag and flap hinges are coincident and 
placed 3.0 in. from the center of rotation, with the pitch bearing placed directly outboard of the hinges. Rotary 
potentiometers placed on the hub and geared to the blade cuffs permit measurement of blade lead-lag, flap, and 
pitch angles. Lead-lag dampers are used on the hub with an effective damping output of 980 in-lb-s/rad. Kinematic 
coupling due to the control system configuration induces a pitch-flap coupling ratio of 0.5, flap up-pitch down; 
that is, the coupling produces a 0.5° nose-down blade pitch for each degree of upward (positive) blade-flap 
displacement. 

3.3.2. Rotor Blades 

The BRR blades are untwisted, rectangular planform research blades with NACA 0012 airfoils. The blades are 
not representative of those used on any full-scale helicopter; however, they have been designed with model-scale 
dynamic properties similar to those used on full-scale helicopters. The blades have a chord of 3.625 in. and, when 
installed on the rotor hub, have a radius of 56.224 in. 
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Intermediate-drive 
lower sheave 


Figure 3.6. ARES -II rotor pylon consisting of rotor drive, control, and force measurement systems. Drive belts and 
miscellaneous hardware are not shown for clarity. 


4. Analytical Models of ARES-II and BRR 

Coupling dynamic systems with an impedance-matching approach requires the development of impedance 
matrices that describe the load-displacement relationships at the interfaces connecting the individual systems. The 
individual systems are then coupled by imposing load and displacement compatibility conditions on the matrices. 
Therefore, the first step in the development of an impedance-matching method for coupling the BRR to the 
ARES-II testbed is the development of individual analytical models of the two systems. This method has the 
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advantage of allowing the selection of the best analysis for each system to be analyzed. For the ARES-II, a finite- 
element code could have been chosen to develop the model. However, the Dynamic Analysis and Design System 
(DADS), a multibody dynamics code, was selected to be used because DADS has a significant advantage over a 
standard purpose finite-element code in the areas of motion visualization, active control interfaces, and hydraulic 
system modeling. For the BRR model, the second generation version of the Comprehensive Analytical Model of 
Rotorcraft Aerodynamics and Dynamics (CAMRAD-II) has been chosen. Each of these analyses offers strengths 
in their respective areas that could not have been matched were a model of the complete system developed in a sin- 
gle analysis. This section presents an overview of the model formulations in each analysis. The development of the 
impedance matrices using results from these models and the coupling analysis is detailed in section 5. Appen- 
dixes B and C provide the reader with detailed program inputs. 

4.1. DADS Model of ARES-II 

DADS (refs. 9 and 10) is a time- marching computer analysis that simulates physical systems by numerically 
integrating a set of differential algebraic equations assembled by the program according to geometry, mass, stiff- 
ness, damping, joint, and external load information provided by the user. The results of the numerical integration 
are supplied as time histories of body positions and system loads. Although the program is primarily a multibody 
dynamics analysis utilizing rigid-body components connected with joints, flexibility may be included at the joints 
and beam elements or component mode synthesis may be used to model elastic structures. For the ARES-II, a 
model consisting of rigid bodies, flexible joints, and a beam element has been developed. For this model, compo- 
nent mass and inertia characteristics have been defined by using the Pro/ENGINEER solid-modeling, computer- 
aided-design software. Component stiffness and damping characteristics have been defined with known experi- 
mentally determined values and estimates obtained from simple analytical models. The construction of the DADS 
model of the ARES-II is discussed from the ground up in section 4.1.1. Figures 3.2 through 3.6 are useful to the 
reader in visualizing the model developed. Section 4.1.2 details the analysis necessary to determine the time- 
varying actuator lengths for control of the ARES-II motion platform. 

4.1.1. ARES-II Model Construction 

The ARES-II stand is shown in figure 4. 1 modeled with two rigid bodies (the base and yoke) and one elastic 
beam (the post). For wind tunnel testing, the stand is constrained to the floor with elastomeric pads placed at each 
of the four comers of the stand baseplate. These pads permit pitch, roll, and normal displacements of the stand but 
restrain axial, side, and yaw displacements. The pads are modeled in DADS by four spring elements with damping 
and constraints to eliminate the axial, side, and yaw degrees of freedom of the stand. 

The ARES-II longeron is modeled as a rigid body. A revolute joint is placed between the longeron and the 
stand yoke to control rotor pylon pitch. The linear screw-jack actuator placed between the stand and the longeron 
to control shaft angles is not modeled; instead a rotational driver constraint is placed at the location of the shaft 
angle revolute joint to control relative pitch motion between the stand and longeron. 

All brackets and fittings for the attachment of the Stewart platform actuators to the longeron are modeled as 
rigid bodies and are rigidly attached to the longeron. The connections between the actuators and the brackets are 
modeled as universal joints. These joints are oriented such that the rotational degree of freedom about the longitu- 
dinal axis of the actuator is eliminated. Thus, the actuators do not spin but still have the necessary freedom of 
motion to actuate the pylon components. The connection of the piston rod portion of the actuators to the actuator 
bodies is modeled by translational joints. Control elements are also modeled to provide means to drive each actua- 
tor according to a specified schedule. The connection between the actuator rods and the Stewart platform baseplate 
is modeled with spherical joints. 

The ARES-II fixed-system balance is modeled as two rigid bodies (top and bottom plates) connected by a six- 
degree-of-freedom elastic joint with damping. The bottom plate of the balance is rigidly attached to the baseplate 
of the Stewart platform. A large rigid body consisting of the main-shaft housing, the intermediate-shaft housing, 
and the drive motor is rigidly mounted to the top of the fixed-system balance. The main-rotor shaft and sheave, 
intermediate shaft and upper and lower sheaves, and the drive motor sheave are all modeled as rigid bodies with 
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Figure 4.1. ARES-II stand. 


rotational degrees of freedom. The main shaft is connected to the main-shaft housing, the intermediate shaft to the 
intermediate-shaft housing, and the drive-motor sheave to the drive motor through revolute joints. Constraint 
equations are used to model the proper drive reduction ratios between the drive motor sheave and the upper sheave 
on the intermediate shaft, and between the lower sheave on the intermediate shaft and the main sheave. 

The mass of the swashplate, pitch links, and swashplate actuators was chosen to be included in the DADS 
model of the ARES-II. Therefore, these components are massless in the CAMRAD-D rotor system model. Their 
mass is included in that of the main-shaft housing. The interface joint for the swashplate, then, is included in the 
rigid-body model of the main-shaft housing. 

The rotating balance, like the fixed-system balance, is modeled as two rigid bodies connected by a six-degree- 
of-ffeedom elastic joint with damping. The bottom plate of the rotating balance is rigidly attached to the top of the 
main shaft. The mast is modeled as a rigid body and rigidly attached to the top plate of the rotating balance. The 
interface joint for the rotor hub is included in the rigid-body model of the mast. 

4.1.2. Actuator Motion Analysis 

A means is necessary to control the ARES-II Stewart platform actuators so that prescribed rotor pylon motions 
may be imposed. To do so one must assume a center of motion for the rotor pylon, or “simulated center of grav- 
ity.” The equations may then be generated to develop schedules of actuator lengths as a function of the desired 
simulated c.g. motions. To simplify the analysis, a rigid structure is assumed so that dynamic response does not 
interfere with the desired motions. 

To begin the analysis, numbering the actuators as shown in figure 4.2 is useful. Two coordinate systems are 
necessary. One is called the “5-frame” (fig. 4.3), a body-fixed coordinate system on the fixed-system balance with 
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Figure 4.2. ARES-II actuator numbering. 


its origin at the balance center. The second is called the “ L -frame” (fig. 4.2) and is fixed to the longeron with the 
same orientation as the fi-frame when the balance is in the nominal centered position. 

A set of vectors may be defined from the balance center to the joints between the Stewart platform actuators 
and the baseplate. These vectors tire defined as rj through r 6 for actuators 1 through 6 and are referenced to the 
balance coordinate frame (fig. 4.3). For the ARES-II configuration these vectors are 


rj = -7.206i fl - 2.237k fl 

r 2 = -6.821 i B - 2.030k B 

r 3 = -7.348ig - 5.724j fi - 2.629kg 
r 4 = -7.348ig + 5.724j fi - 2.629kg 
r 5 =-7.113ig-5.450jg - 2.862k fl 
r 6 = -7.1 13ig + 5.450jg - 2.862k fi 


(4.1) 


where the dimensions are in inches. Another vector is defined from the balance center to the arbitrarily selected 
simulated c.g. (fig. 4.3). This vector is defined as 


r s 


,‘B* r S 




+ r_ k. 


(4.2) 


Therefore, six vectors from the simulated c.g. to the Stewart platform baseplate joints may be generated with 


where the subscript a denotes the referenced actuator (e.g., r gl is the vector from the simulated c.g. to the base- 
plate joint for actuator 1 ). 

Assume now a desired translational displacement vector for the simulated c.g. with respect to the longeron 
coordinate system L as follows: 


d 


8 



(4.4) 


19 



^3 



(a) Top view. 


Simulated c.g. 



(b) Side view. 

Figure 4.3. Vectors for motion analysis. Views are of fixed-system balance mounted on Stewart platform baseplate. 


To obtain the desired translation of the simulated c.g., each Stewart platform baseplate joint must translate an 
equal distance. Therefore, the actuated displacement of each joint for simulated c.g. translation is 

d a = d * (4.5) 

To calculate joint displacement for rotations about the simulated c.g., however, requires the application of a 
rotation of each vector r ga and then a translation of -r ga . To illustrate, assume a single-axis rotation 0 through 

x g 

the simulated c.g. where the rotation is with respect to the longeron coordinate system L. The matrix equation for 
the joint displacement is 
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The general form of equation (4.6) is 

«U = [*, f Hr, a }-{r, a } 


(4.6) 


(4.7) 
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Therefore, in a mathematical sense, all rotations must be completed before the translation occurs. If a rotation 
sequence of y-x-z (pitch-roll-yaw) is assumed, the complete rotation matrix may be written as 


cos 0^ cos 0, 

- sin 0 x sin 0^ sin 0 7 


-cos 0^ sin 0, 


sin 0 cos 0 

y z 

+ sin 0 r cos 0 sin 0, 

x y 


[R g ] = 


cos 0 sin 0 

y 

+ sin 0^ sin 0^ cos 0^ 


cos 0^ cos 0 Z 


sin 0 sin 0 

y 

- sin 0^ cos 0 y cos 0 Z 


-cos 0 V sin 0 V sin 0 cos 0 cos 0 

X y x Ay 


(4.8) 


Combining the results for actuator displacement due to translation and rotation of the simulated c.g. yields 


{d a } = [I]{d,} + [[R g ]-[I]]{r gfl } (4.9) 

Position vectors may be developed for each actuator when the Stewart platform is placed in its nominal or refer- 
ence position. The position vectors are denoted \ a and represent the relative position of the Stewart platform joint 
of the actuator with respect to the actuator joint at the longeron. These vectors are 

= 1.869i L + 6.580j L -3.987k L 

1 2 = 2.254i L - 6.580j L - 3.780k L 

1 3 = 6.346i L + 1.376j L -4.384k L 

1 4 = 6.346i L - 1.376j L -4.384k L * 

1 5 = - 4.825 j L -6.062k L 

1 6 = 4.825j L - 6.062k L 

Then, the actuator lengths may be determined by 

- |i. +d J 

For control of the DADS model Stewart platform actuators, equations (4.9) through (4.1 1) are used to calcu- 
late schedules for each actuator length. These results are applied as time-varying distance constraints between the 
joints at each end of the actuators. 

4.1.3. DADS Model Results 

Because DADS is a time-marching analysis simulating the physical system through numerical integration, the 
output is a set of time histories. Thus, no eigenanalysis is available to determine system frequencies and modes. 
Instead, it is typical to “excite” the DADS model and measure the response in a manner similar to physical-model 
shake testing. Alternatively, a displace-and-release (or “pluck”) excitation may be used to identify the basic free- 
response characteristics of the model. 

The damped natural frequencies of the DADS model of the ARES-II have been identified by displacing the 
hub interface joint, then releasing and allowing free- vibration response. Fast Fourier Transforms of the loads mea- 
sured at the fixed-system balance during the free response are shown in figure 4.4. For the axial and pitch 
response, the ARES-II was excited by plucking the hub interface joint axially. Side and roll responses were gener- 
ated by plucking the hub joint to the side; the normal response, by plucking the hub joint vertically; and the yaw 


(4.10) 


(4.11) 
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response, by plucking with a yawing moment at the hub joint. The FFTs in figure 4.4 have been used to identify 
the damped natural frequencies and to estimate the modal characteristics of the ARES-II system; these are pre- 
sented in table 4. 1 . 

Significant modal coupling exists for the ARES-II elastic modes because the c.g. of the rotor pylon is above 
and forward of the fixed-system balance flexibility. Therefore, axial motions of the ARES-II rotor pylon generate 
pitch response. The associated vertical acceleration of the c.g. due to pitching then generates normal (vertical) 
response. Similarly, roll and yaw motions couple with side motions of the rotor pylon. The information presented 
in table 4.1 are estimates based on the results of figure 4.4 and a physical understanding of the ARES-II. No exper- 
imental data were available that defined the stiffness and damping of the ARES-II elastomeric stand mounting 
pads. Therefore, a trial-and-error approach was used to define these properties to obtain rigid-body pitch and roll 
frequencies of 5 Hz with 8 percent critical damping. The elastic modes are listed with the component estimated to 
be the most dominant first and the least dominant last. For example, the first elastic mode at 58 Hz is considered 
primarily an axial mode with pitch and normal coupled response. To obtain a more accurate assessment of the 
ARES-II frequencies and modal identification would require the development of a finite-element model and an 
eigenanalysis. 

4.2. CAMRAD-II Model of the BRR 

CAMRAD-II is a state-of-the-art rotorcraft analysis that can model a wide range of rotorcraft configurations. 
For this study, CAMRAD-II was used to develop a complete aerodynamic and dynamic model of the BRR. The 
kinematic effects of the swashplate control system were included in the model. However, because the weight of 
the swashplate and control system was included in the DADS model, these components were considered massless 
in the CAMRAD-II model. 

The development of a detailed CAMRAD-II model consists of two stages. The first stage is the development 
of the “rotorcraft shell” model. The CAMRAD-II rotorcraft shell allows the user to input a general set of rotorcraft 
information to construct models of most common rotor system configurations. The shell effectively eliminates 
much of the tedious work involved in generating a model because it constructs much of the detailed system model. 
The second stage is the development of the “core” model. The core model allows the user to supply detailed infor- 
mation regarding the system to be modeled; this permits the generation of configurations that cannot be fully con- 
structed by the shell. The development of the core model is, however, more complicated and may easily require 
input data an order of magnitude more than that required by the shell. 

4.2.1. BRR Shell Model 

The rotor system model of the BRR was initially generated by using rotorcraft shell input. Because each blade 
is, in theory, identical, only one blade is modeled. Phasing routines internal to CAMRAD-II allow for the calcula- 
tion of system response by using all four blades. The hub is modeled as a rigid body inboard of the flap and lead- 
lag hinges. A second rigid body is modeled between the hinges and the pitch bearing. A third rigid body extends 

Table 4.1. Damped Natural Frequencies and Modes of ARES-II 


Frequency, 

Hz 

Mode 

5 

Rigid pitch 

5 

Rigid roll 

17 

Rigid normal 

58 

Elastic axial/pitch/normal 

67 

Elastic side/roll/yaw 

79 

Elastic normal/axial/pitch 

90 

Elastic yaw/side/roll 

131 

Elastic normal/pitch/axial 

163 

Elastic side/roll 


23 




from the pitch bearing to the outboard end of the rotor blade cuff. The blade is modeled by using 5 finite-element 
beams with 19 aerodynamic panels distributed radially. The kinematic effects of the control system are included 
by providing detailed placement of the pitch horn and pitch link attachment points. 

4.2.2. BRR Core Model 

Enhancement of the model generated by the rotorcraft shell is necessary to permit the implementation of pre- 
scribed hub motions. These hub motions, as shown in section 5.3.2, are necessary to develop the rotor system 
impedance matrices for the coupling analysis. Additionally, the hub mass and inertia characteristics calculated by 
the rotorcraft shell are generally inaccurate. Therefore, the core model was used to prescribe the proper hub mass 
and inertia. 

Figure 4.5 shows a schematic of the model developed. The boxes in the figure represent components available 
for use in CAMRAD-II models. The component type is listed at the bottom of each box (e.g., rigid body, Fourier 
series) and the component name, if any, is listed at the top of each box. The circles in figure 4.5 represent joints 
between structural components of the CAMRAD-II model. The value inside the circle specifies the number of 
motion degrees of freedom that the joint permits. Labels next to the joints identify the joint function. 

The CAMRAD-II model that was generated represents the ARES-II/BRR in configuration only. The stand and 
pylon components shown in figure 4.5 are dimensioned to represent the ARES-II but have no stiffness or mass 
properties. Thus, the ARES-II portion of the CAMRAD-II model is strictly a kinematic, but not dynamic, repre- 
sentation of the testbed. It is modeled only to impose the BRR hub motions required for development of the rotor 
impedance matrix described in section 5. 3. 2. 3. The BRR model, however, is a complete dynamic and aerodynamic 


Pitch bearing 



Figure 4.5. Schematic of CAMRAD-II model of BRR. 
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representation of the rotor system. As shown in figure 4.5, the model consists of a rigid body component named 
“stand” which is rigidly attached to the ground. To permit prescribed hub motion, a rigid body component named 
“pylon” is installed at the top of the stand component with a six-degree-of-freedom (6-dof) joint. All rotor system 
components as generated by the rotorcraft shell in section 4.2.1 are attached to the pylon component. The dashed 
box at the left of the figure shows the control input portion of the model. Harmonic input amplitudes may be pre- 
scribed for each of the six possible joint motions. These amplitudes are input to a Fourier series component that 
generates a set of six azimuthally varying signals by using the equation 


8 

7(\|/) = ^ ( H n cos n\\f + H n sin /i\|/) 
n=0 


where H n and H n are the harmonic inputs, T(ij/) is the output signal, and only the mean and first eight harmon- 
ics of motion are considered. The six output signals are combined into a single vector by a “summer” component 
that inputs the signals to the 6-dof joint between the stand and pylon components. Actuator features internal to the 
6-dof joint process the vector and generate linear and angular motions of the pylon and, thus, the hub. 

Output components are included in the model to display hub motion, hub loads at the hub interface joint, and 
swashplate loads at the swashplate interface joint. These output components resolve loads in the fixed-system 
frame for the mean and 4P and 8P harmonic forces and moments. This procedure offered the advantage of having 
to process only 4P and 8P loads in the coupling procedure instead of processing 3P, 4P, 5P, 7P, 8P, and 9P loads as 
would be required if the loads were resolved in the rotating frame. Thus, an effective revolute joint is placed 
between the rotor hub and the mast of the ARES-II model, which although not representative of the configuration 
of the physical model, is effective in producing efficient results due to the reduction in excitation frequencies 
considered. 

5. Dynamic Analysis Using Impedance Matching 

Dynamic analysis of systems through impedance-matching approaches can be useful. However, as with any 
engineering analysis approach, one must be careful to consider the limitations and assumptions of the application. 
For impedance matching, two such limitations stand out. First, impedance matching is a frequency-dependent 
technique. Therefore, impedance characteristics must be defined for each component to be coupled at each fre- 
quency of interest. Second, impedance-matching approaches are developed with linear theory. Therefore, only 
systems that are linear by nature or whose characteristics are linear in a region of interest are suitable for analysis 
by impedance matching. For rotorcraft application, in which one may be interested in coupling a rotor system to a 
fuselage model, the first of these limitations is generally not of considerable concern. Because the rotor system is 
the primary contributor to helicopter vibratory loads and the rotor-generated forcing functions applied to the fuse- 
lage are harmonic in nature, only a few such frequencies of interest exist in a rotorcraft application. Typically, 
rotor dynamics applications are restricted to those frequencies below 9P for a four-bladed rotor system. Therefore, 
for the Basic Research Rotor, which passes to the fixed system only steady loads and loads at 4P and 8P, only two 
frequencies exist for which impedance characteristics must be developed (i.e., 4P and 8P). The second limitation 
of impedance matching, however, is of greater concern in rotorcraft analysis. Although fixed-systems structures 
(i.e., fuselages, testbeds, and whirl stands) are typically modeled by using linear theory, the rotor system is inher- 
ently nonlinear. This nonlinearity is of particular concern at low flight speeds (p, <0.15) where blade tip vortices 
impinge on the following blades, and high flight speeds (p > 0.30) where rotor blade stall produces nonlinear 
response. Therefore, the application of impedance matching to couple a rotor system to a fuselage or test stand 
requires that the rotor impedance model be developed about a steady-state flight condition of interest. Addition- 
ally, the results must be carefully checked to ensure that rotor system response remains within a linear range. 
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The term “impedance” often has different meanings which depend upon the application used. For this study, 
impedance is defined as the force-displacement relationship that exists for a system degree of freedom when an 
interface joint location is being driven by a prescribed oscillatory motion at a particular frequency of interest. Also 
of interest is system “mobility,” which is defined as the displacement-force relationship that exists for a system 
degree of freedom when an interface joint location is being subjected to a prescribed oscillatory force at a particu- 
lar frequency of interest. Both impedance and mobility, as shown in section 5.1, may be represented by square 
matrices. Further, system mobility is shown to be the inverse of system impedance. Therefore, the term “imped- 
ance matching” is appropriate for system coupling using either impedance matrices, mobility matrices, or both. 

Section 4 discussed the DADS model of the ARES-II and the CAMRAD-II model of the BRR. This section 
discusses how these models are used to develop impedance representations of each system and details the 
impedance-matching approach developed to couple the two systems. For this study, two interface joints are used: 
one at the joint between the hub and mast of the main-rotor shaft and the other at the joint between the swashplate 
and the main-rotor shaft. 


5.1. Sample Coupled Model 

A simple model is used to demonstrate the principles of the coupling analysis. This model is shown in 
figure 5.1 . Two systems are shown which are to be coupled at the interface joint /. System A consists of three one- 
degree-of-freedom (1-dof) masses. Spring k\ and damper are placed between the ground and my\ an actuator is 
placed between masses wj and m 2 ; and another spring & 2 and damper c 2 are placed between masses m 2 and m 3 . 
The interface joint for the system is placed on m 3 . System B consists of one 1-dof mass with a spring £3 and 
damper c 3 . The interface joint for system B is placed at one end of the spring-damper so that the spring-damper 
lies between the interface joint and the mass. System A is a 2-dof system due to the constraint that exists because 
of the actuator. This constraint equation may be expressed as 


*2 = *1 + l A 


(5.1) 



System A 



Figure 5.1. Sample model. 
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The remaining two generalized coordinates for system A include x j and xy The equations of motion for these two 
coordinates may be easily generated. They are 


(mj + m 2 )x j + m 2 lA + (Cj + c 2 )i] - c 2 (x 3 - i A ) + (k j + k 2 )x , - k 2 (x 3 - l A ) - 0 
m 3 x 3 + c 2 (x 3 - x ; -l A ) + k 2 (x 3 -x x - l A ) = Fj(t) 

Because the actuator may be used to force system motion, the equations of motion are written as 


( 5 . 2 ) 


(m j + m 2 )x j + (c j + c 2 )jc j — c 2 x 3 + (&] + k 2 )-Xj k 2 x 3 — m 2 l A c 2 l A k 2 l A 
m 3 x 3 + c 2 (x 3 -i,) + k 2 {x 3 -x l ) = F,(t) + c 2 i A + k 2 l A 

Now let 

F ,(t) = F, cos Qt + F, sin Qf 

‘ * c 1 s 

> 

l A (t) = l A cos Q t + l A sin Qt 


( 5 . 3 ) 


( 5 . 4 ) 


where F, F, , / , , and / . are constant and Q is the frequency of the applied force or actuation. If steady-state 

*c *s 

responses are assumed in the form 


x } (t) = Xi cosQf + ^i sin Qf 

J l c l s 

► 

jc 3 (t) = x 3 cosQt + ^3 sin Qt 

a matrix expression for the solution of the steady-state response may be written as 


( 5 . 5 ) 


1- 



*■ 


"V 

2 

-Q (mj + m 2 ) + 

Q(c, +c 2 ) 

~ k 2 



X K 

-£2(cj + c 2 ) 

2 

-£2 (^| + m 2 ) + + &2 

C 2 & 

k 2 

< 


~ k 2 

-Qc 2 

k 2 - £l 2 

Qc 2 



Qc 2 

-h 

—£lc 2 

k 2 - 


k 


Q m 2 l Ac c 2 Q.l A * k 2 l A ' 

^ m 2^A s + c 2^A c ~ k 2^A s 
Fj +k 2 i A +c 2 ai A 


Fj - c 2 ai A+ k 2 i A 


( 5 . 6 ) 
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5.1.1. Mobility Matrix Development 

A numerical example is used to illustrate the development of the sample model mobility matrices. Let 
m \ = m 2 = m 3 = 1 0 slug, k\ = &2 = 10 0 lb/ft, and cj = C 2 = 2.0 lb-s/ft. To develop the mobility matrices for mj and 

m 3 , the force input at the interface is set to a unit cosine force (F. =1 .0, F, = 0.0) and the actuator is held fixed 

*c 1 S 

(l A ^ = l A = 0.0). For an excitation frequency of £2 = 2n rad/s, the vector x resulting from the solution of 
equation (5.6) is 



Similarly, the result when a unit sine force ( Fj 
( l A =l A =0.0) is 

c s 


-0.001343 

-0.006917 

► (5.7) 

-0.025186 

0.012510 

0.0, Fj = 1.0) is applied and the actuator is held fixed 


V 


0.006917 

*1 


-0.001343 

1 J 

► = < 

> 



-0.012510 

*3 


-0.025186 

J J 




With the results from equations (5.7) and (5.8), a mobility matrix may be generated for either m, or m 3 . Because 
the interface joint is on m 3 , it is necessary that the mobility matrix for this mass be developed. The development of 
the mobility matrix for m\ is optional for the coupling analysis but is useful if complete system motion is 
desirable. 

The general method for developing a mobility matrix is to place the cosine and sine components of the degree- 
of-freedom displacement due to a unit cosine forcing function at the interface joint in the first column of the 
mobility matrix. Likewise, the second column of the matrix is filled with the cosine and sine components of the 
degree-of-freedom displacement due to a unit sine forcing function at the interface joint. For m 3 , then, the mobility 
matrix is 


[A 3 ] 


-0.025186 -0.012510 
0.012510 -0.025186 


(5.9) 


Because the system is linear, the mobility matrix may be developed by using only the results from either the cosine 
force or the sine force. Throughout this work, a sine excitation has been used so that the mobility matrix for the 
mth degree of freedom may be generated by 





(5.10) 


Therefore, by using a mobility matrix, the motion of a degree of freedom may be determined if the applied force at 
the interface is known; that is. 
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(5.11) 


< 


> = 






5.1.2. Motion Matrix Development 

Besides generating the mobility matrix for the sample system A, the effect of actuator motion on system per- 
formance is necessary. These effects may be included as a matrix component as well. For this study, this matrix is 
called the “motion matrix” and will be denoted by [B m ] for each degree of freedom m of interest. The motion 
matrix may be generated in a manner similar to the mobility matrix. With the numerical example, equation (5.6) is 
solved for an absence of force input (F / = 0.0, F j = 0.0) while either a unit cosine motion (l A = 1.0, Ia s = 0.0) 
or a unit sine motion (l A = 0.0, l a = TO) is effected by the actuator. The motion matrix is assembled in the same 
manner as the mobility matrix. For the sample system, the resulting motion matrix for m 3 is 


[b 3 1 


-0.047323 -0.220797 
0.220797 -0.047323 


(5.12) 


Note that for this application an “ideal” actuator is assumed so that pure cosine and sine motions are provided by 
the actuator. 

Now, by using superposition, a matrix equation can be developed which defines system motion due to the 
application of an external force at the interface joint and due to actuation. This equation is 


X I 


r 

F, 

>+[B 3 ]< 

' ’ 

U 

1 c 

■ = [A 3 ]- 

C 

c ► 

x l s 

/V 

> 


(5.13) 


5.1.3. Impedance Matrix Development 

Now an impedance matrix is developed for sample system B. The equation of motion for this simple system is 

m A x^ + c ^x A + k.T> x A = Ct ) Xj + k^Xj (5.14) 


If 


x, = x. cos fit + x. sin fl f 

1 ‘c 's 

x 4 = x 4 cos fit + x 4 sin fit 


(5.15) 


a matrix equation may be written for the system as 


- Q 2 m 4 

flc 3 

« 



k 3 Qc 3 

< 

x, ‘ . 

-flc 3 

& 3 - Q m 3 


-* 4 , 


_-flc 3 k 3 


X ls 


For this example, let jfc 4 = 10.0 lb/ft, c 3 = 2.0 lb-s/ft, m 4 = 1.0 slug, and fl = 2k rad/s. Prescribed unit cosine and 
unit sine motions of the interface joint are applied individually and equation (5.16) solved for {x 4 }. For the pre- 
scribed unit cosine (x, = 1 .0, x, = 0.0) motion, 

l C 'i 
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(5.17) 


< 


-0.133286 | 
0.483109 j 


For the prescribed unit sine ( x . = 0.0, x, = 1.0) motion, 

1 c 1 s 


< 


-0.483109 

-0.133286 


(5.18) 


However, to complete the coupling analysis, the force at the joint Fj is necessary, not the motion of m 4 . This force 
is described by 


F l (t) = k 3 (x 4 -x / ) + c 3 (x 4 -x l ) (5.19) 

or 

F I = - */ ) + «c 3 (jc 4 - X , ) 

r (5.20) 

F I S = k 3^ x 4 s - x l) + klc 3^~ x 4 c + x l) 

The impedance matrix [C] may then be developed in a manner identical to the mobility and motion matrices 
described previously. The impedance matrix is thus determined from the results of equations (5.17), (5.18), and 
(5.20) as 


[C] 


-5.2619 -19.0724 
19.0724 -5.2619 


(5.21) 


Finally, the following force equation for the interface joint may be written by using the impedance matrix: 


F I 


X I 

c 

■ = [C]. 

c ► 

x h 


(5.22) 


5.1.4. Coupled Analysis 

To develop the analysis to couple systems A and B of figure 5.1, some compatibility condition must exist at 
the interface joint I of each system. For this analysis, the compatibility conditions include displacement and force 
at the interface joint {x/} and {F/} , respectively. Substituting equation (5.22) into equation (5.13) to ensure that 
force compatibility exists at the interface joint and solving for {x 7 } to ensure that displacement compatibility exists 
as well give 


X I 


r -| 1 

i a 



[i]-[A 3 j[C]J [b 3 ]< 

c ( 


(5.23) 
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Substituting equation (5.23) into equation (5.22) yields the following equation for the interface force: 


Fj 


1 ' 

l \ 

c 

F >s 

> = [C] 

[i]-[A 3 ][C]J [b 3 ]< 

c 

V 


(5.24) 


Equations (5.23) and (5.24) represent the interface displacement and force for the completely assembled and 
coupled system. The results are valid for the frequency for which the mobility matrix [A 3 ], the motion matrix [B 3 ], 
and the impedance matrix [C] were developed. For the coupled sample system at an actuation frequency C2 of 
2jt rad/s, the results for l A = 1 .0, l A = 0.0 are 



■ 


J 0.044853 1 

X 's 

1 0.177474 j 


► 

Fj 

f -3.6209 ] 

> §C 

= \ \ 

f f 

\ -0.0784 J 




(5.25) 


Note also that x 3 = x t so that the motion of m 3 has been obtained as well. Other system motion results may be 
obtained by developing the following mobility matrix for mj: 


[Aj] 


-0.001343 0.006917 
-0.006917 -0.001343 


and by using equations (5.1), (5.1 1), and (5.16). The results are 


= [A,] 


0.004321 

0.025151 


> + < 


1.004321 

0.025151 


U 3 -Q m 4 


Qc-, 


-Qc 3 k^-Q m 4 


- -1 



r A 


k ^ 

< 

X 'c . 


-Qc 3 k 3 


X ls 


(5.26) 


(5.27) 


(5.28) 


-0.091718 

-0.001986 
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(5.29) 



5.2. Comments and Observations Regarding Impedance Matching 

Several characteristics of the impedance-matching method developed in section 5.1 are noteworthy. First, the 
method as presented is only useful for dynamic systems undergoing steady-state forced response at discrete fre- 
quencies. Second, a 1-dof response of a single location on a system is representable by a 2 x 2 matrix for each fre- 
quency of interest. Therefore, for complete 6-dof response of an interface joint or other location, a 12 x 12 matrix 
is necessary. Multiple 12x12 matrices are necessary for multiple excitation frequencies. Third, all response matri- 
ces may be developed with similar techniques. 

It is evident as well that a relationship exists between a mobility and an impedance matrix. To verify this rela- 
tionship, one need only recall variations of equations (5.11) and (5.22) such as 

{x} = [A]{F}1 
{F} = [C]{x}J 

In practice, the only reason to use both a mobility matrix and an impedance matrix is to identify the manner in 
which the response matrices were developed. A mobility matrix is defined by applying a forcing function and 
examining the motion response. The impedance matrix is developed by applying a motion and examining the 
force response. 

5.3. Application of Impedance Matching to ARES-II and Basic Research Rotor 

Most techniques developed in section 5.1 are directly applicable to the analysis of the ARES-II testbed and the 
BRR. However, some modifications are required to model the rotor system properly. The following sections detail 
the procedures necessary to couple the two systems. 

5.3.7. ARES-II Impedance Modeling 

Two load paths exist for the transfer of loads between the ARES-II and the BRR. These paths include the rotor 
hub connection to the mast and the swashplate connection to the main-rotor shaft and main-shaft housing. The 
largest load transfers occur at the rotor hub connection. The swashplate, however, restrains the pitch links and thus 
carries the loads generated by blade pitching moments and control. As shown, the swashplate loads may, at times, 
represent a significant portion of the harmonic load transfer from the rotating system to the fixed system. There- 
fore, to neglect these loads would introduce an uncertainty as to the validity of the model. Furthermore, because 
one goal of this study is to define the loads sensed by both the rotating and the nonrotating balances, including the 
swashplate load path is necessary because the rotating balance is unable to measure these loads. Note that includ- 
ing the swashplate loads in the analysis represents a departure from the standard practice evident in the literature 
where, almost universally, the swashplate loads are neglected. 

The DADS model of the ARES-II described in section 4.1 is used to develop the necessary impedance charac- 
teristics of the hub and swashplate interface joints. Within the DADS model, the two interface joints at the hub and 
swashplate are included as points of interest to ease load application and to identify system response. These two 
interface joints and their orientation are shown in figure 5.2. The main-rotor shaft assembly in the DADS model 
can rotate. However, for this study resolving all forces and moments in the nonrotating system is convenient 
because the number of frequencies required for analysis is reduced to two: 4 per rev and 8 per rev. Therefore, the 
main-rotor shaft is held fixed, and a rotational joint is assumed at the interface between the mast and the rotor hub. 
As a result, none of the rotational inertia of the main-rotor shaft assembly or drive system is included in the model. 
These effects, however, are considered minimal when compared with the rotational inertia of the rotor system. 
Additionally, because no drive system flexibility is modeled, including the drive system inertias would be rela- 
tively meaningless. 

5.3.7.7. Harmonic analysis . In contrast with the dynamic analysis of the ARES-II testbed presented in sec- 
tion 4.1.3 in which all frequencies were excited with an initial displacement, the development of an impedance- 
matching analysis requires examination of the system steady-state response due to harmonic excitation. Therefore, 


32 




Figure 5.2. Location and orientation of interface joints on ARES-II pylon. 


a significant amount of time-history data is generated with DADS throughout the analysis of the ARES-II testbed 
to relate the excitation and response at discrete frequencies. Additionally, a reduction of the steady-state responses 
to constant coefficients is necessary to develop the required impedance matrices. Because DADS is a simulation 
analysis that uses a time-marching integration solution, all transient response must decay before an accurate 
assessment of the steady-state amplitude of motion may be obtained. These steady-state amplitudes are often 
obtained with an FFT algorithm that generally requires a substantial length of time-history data for accurate 
analysis. However, since the impedance-matching analysis requires discrete frequencies of interest, an alternative 
approach is suggested. 

A harmonic analysis of the time-history data generated by DADS is used that is similar to the conditional sam- 
pling method used in reference 7. With this method, only one complete cycle of data is necessary for analysis. The 
length of the cycle must be equal to the period of the rotor rotational frequency. For the BRR, a rotor system oper- 
ating at 660 rpm, this period is 0.09091 s. For accurate results, this method also requires that most transient 
response has decayed. However, the harmonic analysis is much less affected by any remaining transients than the 
FFT. 
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For the analysis, the time-history window W is divided into 7 time steps (A t = W/f) and evaluation takes place 
over 7 + 1 samples. Sample 1 of the window is aligned in the time history with y = 0° and sample 7 + 1 with 
\)/ = 360° (reference blade over the tail for the first and last time steps). The index j is used to denote the time step 
= 1 to 7 + 1), and a weighting factor wj is necessary for the analysis. At the window endpoints (j = 1 and 
j = 7 + \),wj= MU, and elsewhere, Wj = 1/7. The components of each harmonic may then be obtained by 


J + 1 

fc n = X 2w j fj cos n Vj 


j= 1 
J + 1 


f S n = X 2w j fj 
j = 1 


sin n\\r ■ 


( 5 . 30 ) 


where n is the harmonic number, yj is the amplitude of the signal at the jth time step in the time history, \\ij is the 
azimuth of they'th time step, and f c and f s are the cosine and sine amplitude components for the nth harmonic 
of the time history. 

For this study, it has been found that a simulation time of only 0.5 s is required for the decay of transient 
response before good quality harmonic analysis results are obtained. Therefore, the time-history window from 
which results have been obtained runs from 0.5 to 0.59091 s; this, compared with the 2.5 s of data required for FFT 
analysis, represents a significant savings in computation time. This harmonic analysis approach has been used for 
all cases in which time-history data analysis was required throughout this study. 

S.3.1.2. ARES-II mobility matrix. Because applying a harmonic force and observing the associated displace- 
ments throughout the DADS model is simpler than applying a pure harmonic displacement and observing the 
forces, the necessary mobility matrices are generated directly with the DADS model of the ARES-II. To generate 
a full set of mobility matrices at the two frequencies of interest, 4P and 8P, requires 12 independent executions of 
the analysis. Six runs are required in which each force and moment is applied to the hub interface joint using 
superimposed 4-per-rev and 8-per-rev forces and moments. An identical set of force and moment runs are required 
for application of load to the swashplate interface joint. Because the analytical model of the ARES-II is com- 
pletely linear, only sine excitations were used to develop the mobility matrices. 

The general form for the mobility equation (eq. (5.1 1)) may be written for a 6-dof interface joint as 
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where the displacement and rotation column vectors and the force and moment column vectors are each cosine and 
sine pairs (i.e., {x} = {x x } T , {F x } = { F F } T ), and the mobility matrices are each 2x2 and represent the 

displacement due to the applied force as described by the subscripts (i.e., [ A J0 / y 1 is the mobility matrix for the 
x-displacement due to a force applied in the y-direction). Therefore, a single DADS run in which a forcing 
function is applied to an interface joint allows the generation of one matrix column of the large matrix in 
equation (5.31). 

The eight 12x12 mobility matrices generated for the ARES-II are given in the following equations: 


For hub interface joint response due to 4P hub interface joint force and moment excitation. 
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(5.32a) 


For hub interface joint response due to 4P swashplate interface joint force and moment excitation. 
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For swashplate interface joint response due to 4P hub interface joint force and moment excitation. 
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For swashplate interface joint response due to 4P swashplate interface joint force and moment excitation, 


t A W =10 ^ x 
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For interface joint response due to 8P hub interface joint force and moment excitation, 


(5.32d) 


f A W = 10 ^ x 
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For hub interface joint response due to 8P swashplate interface joint force and moment excitation, 


f A «/j= 10 * 
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For swashplate interface joint response due to 8P hub interface joint force and moment excitation, 


(5.320 


[ a W = 10_4x 
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For swashplate interface joint response due to 8P swashplate interface joint force and moment excitation, 


I A W = 104 * 
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The subscripts of the mobility matrix designate the type of mobility information that the matrix contains. The sub- 
script to the left of the slash indicates the joint at which the response is measured; the subscript to the right of the 
slash indicates the joint at which the forcing function is applied. The numeric subscript indicates the frequency of 
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excitation in per rev. For example, equation (5.32b) represents the mobility of the hub joint due to a 4P forcing 
function applied at the swashplate joint. 

Significant physical insight regarding the ARES-II response may be obtained by studying the mobility matri- 
ces presented in equations (5.32). The reader is referred to equation (5.32a), which presents the hub joint response 
due to 4P hub joint excitations. The first two columns of this matrix identify the response of the hub joint when 
subjected to a 4P axially directed load. As shown in rows 1 and 2, the axial response is dominant. Rows 5, 6, 9, and 
10, however, indicate that normal (vertical) and pitch responses are excited as well. Similarly columns 3 and 4 
identify the response of the hub joint when subjected to a 4P side-directed load. The response is primarily to the 
side; however, coupling exists both in roll (rows 7 and 8) and in yaw (rows 1 1 and 12). 

Out-of-phase results are evident in the 8-per-rev mobility matrices (eqs. (5.32e) through (5.32h)). In 
equation (5.32e), for instance, the 8P side-directed load at the hub results in a side response that is 70° out of phase 
with the force; this is due to the proximity of the forcing function frequency (8P or 88 Hz) to the elastic yaw-side- 
roll mode at 90 Hz (refer to table 4.1 and fig. 4.4). 

Two further comments regarding the mobility matrices of equations (5.32) are warranted. First, these matrices 
have not been generated with a unit sine forcing function as used in section 5.1.1. Since the ARES-II is a relatively 
stiff structure, minimal response would have been generated by a unit load. Therefore, larger excitation loads were 
used and the results presented in equations (5.32) normalized to obtain the displacement-per-unit-load values. Sec- 
ond, the results are not entirely indicative of a strictly linear system. Note the submatrices labeled A in 
equation (5.32a) and the submatrices labeled B in equations (5.32b) and (5.32c). Ideally, the two submatrices 
labeled A would be identical as would the submatrices labeled B. However, these minor differences exist because 
of small errors inherent in the DADS and harmonic analysis solutions. The differences are considered minor and 
no action has been taken to adjust them or to prevent their occurrence. 

Finally, a mobility matrix is generated for the entire ARES-II system with the matrices in equations (5.32). A 
two-step process is used for clarity. First, a mobility matrix is generated for both the 4P and the 8P equations as 
follows: 


I a s/// 4 ] I a s/s 4 1 

[ A S/tf g ] t A S/Sg] 


Second, the final assembly of the system mobility matrix occurs as 


[A] 


[A 4 ] [0] 
[0] [A g ] 


(5.33) 


(5.34) 
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The completed mobility equation, then, is 


t X * 4 > 


< F * 4 > 

{X* > 


{ P 5 4 > 


= [A] 

4 

< x * 8 > 




A 

_ {F V 


where 

i*H } = 


Cm x Tj y ii y it z ii Z.IJ ^7 

H nc H ns H nc H ns H nc H ns X H„ f x H mg ^ X, "nc "n*J 


(5.35) 


< X *„> = 
<*V = 

{F S > = 


5., x s„. x, x. z s mf z s 


e e _ e .. 0 .. 0 , e . 


nr ns 


nc X, V X, X S_„ *S_„ >S, 


nc ns nc ns nc ns J 


F F. F„ F„ F, F, M r M r M u M v M r M 7 


X H ‘ X H " >« ~ y H * Z H Z H X H X H yH >H_. Z H 

nc ns nc ns nc ns nc ns nc 


nc ns 


F F F F F. F, M r Af Af , M v Af M 


*s ys y$ 


nc ns nc ns nc ns nc 


y s ys z s z s 

nc ns nc ns 


5.3.I.3. ARES-II motion matrix. The motion matrix is developed by actuating the ARES-II model in each of 
the three translational and three rotational motions individually, then processing the response at the hub and 
swashplate interface joints. Thus, only six DADS runs are required, one for each motion. The actuated motion 
analysis developed in section 4.1.2 is used to generate the actuator schedules for the motions. For translational 
motions a nominal amplitude of 0.01 in. is used, and for rotational motions a nominal amplitude of 0.05° is used. 
All motions are actuated at 4 per rev by using sine excitations. The resulting response at each of the interface joints 
is analyzed by using the harmonic analysis, then scaled by the appropriate input amplitude to provide a 
displacement-per-unit-displacement relationship. These relationships are used to generate a motion matrix for the 
hub [ ] and the swashplate [B 5 ] interface joints. These two component matrices are shown in the following 

equations: 

For hub interface joint response due to 4P pylon actuation, 
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0.0088 

1.3329 


(5.36a) 
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For swashplate interface joint response due to 4P pylon actuation. 


f»J = 


0.6000 

- 0.0573 

0.0006 

0.0000 

0.3162 

- 0.0617 

0.0013 

- 0.0013 

14.1799 

- 0.5953 

- 0.0003 

0.0009 

0.0573 

0.6000 

0.0000 

0.0006 

0.0617 

0.3162 

0.0013 

0.0013 

0.5953 

14.1799 

- 0.0009 

- 0.0003 

0.0107 

- 0.0010 

0.5498 

- 0.0446 

0.0107 

- 0.0010 

- 8.9523 

- 0.1156 

0.1219 

- 0.0116 

0.2184 

2.5199 

0.0010 

0.0107 

0.0446 

0.5498 

0.0010 

0.0107 

0.1156 

- 8.9523 

0.0116 

0.1219 

- 2.5199 

0.2184 

0.0038 

- 0.0054 

0.0000 

0.0000 

1.0302 

- 0.1050 

0.0007 

- 0.0006 

0.7257 

- 0.6732 

- 0.0002 

0.0010 

0.0054 

0.0038 

0.0000 

0.0000 

0.1050 

1.0302 

0.0006 

0.0007 

0.6732 

0.7257 

- 0.0010 

- 0.0002 

- 0.0001 

0.0000 

- 0.0081 

0.0019 

- 0.0001 

0.0000 

1.0808 

0.0035 

- 0.0014 

0.0001 

- 0.0105 

0.0487 

0.0000 

- 0.0001 

- 0.0019 

- 0.0081 

0.0000 

- 0.0001 

- 0.0035 

1.0808 

- 0.0001 

- 0.0014 

- 0.0487 

- 0.0105 

0.0112 

- 0.0029 

0.0000 

0.0000 

0.0266 

- 0.0045 

0.0001 

- 0.0001 

1.5185 

- 0.0430 

0.0000 

0.0001 

0.0029 

0.0112 

0.0000 

0.0000 

0.0045 

0.0266 

0.0001 

0.0001 

0.0430 

1.5185 

- 0.0001 

0.0000 

0.0000 

0.0000 

- 0.0112 

0.0007 

0.0000 

0.0000 

- 0.0395 

0.0022 

- 0.0001 

0.0000 

1.3326 

- 0.0088 

0.0000 

0.0000 

- 0.0007 

- 0.0112 

0.0000 

0.0000 

- 0.0022 

- 0.0395 

0.0000 

- 0.0001 

0.0088 

1.3326 


The motion matrix defining the 4-per-rev response is 


(5.36b) 


[B 4 ] = 


|[B W J [0] 

[0] [Bj 


(5.37) 


Therefore, since 8-per-rev motions are not actuated, developing additional matrices for the 8P response is unnec- 
essary. The completed ARES-II system motion matrix is 


[B] 


[B 4 ] [0] 

,[0] [II 


(5.38) 


where the identity matrix is necessary to expand the matrix size to 48 x 48 like the mobility matrix [A]. 

To drive the ARES-II by using the motion matrix [B] a control vector {U} is necessary. This control vector is 
a 12-element vector defining the cosine and sine components of the desired motions of the simulated c.g. referred 
to in section 4.1.2. However, because of the manner in which the motion matrix has been developed (i.e., [B w ] 
and [B 5 J on the diagonal), it is necessary to use the vector twice. Therefore, the completed motion equation is 4 


< X tf 4 > 

< x s 4 > 

< x v 


[B] 


{U} 

{U} 

{ 0 } 

.{0}. 


where the vectors X were defined in equation (5.35), and 


{U} 





U. 








(5.39) 


5.3.2. Basic Research Rotor Impedance Model 

The CAMRAD-II model discussed in section 4.2 is used to develop the impedance model of the BRR. The 
development of the impedance matrix is, however, more complicated than that required for the sample system of 
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section 5.1.3. Several factors must be considered. First, the rotor system is an inherently nonlinear system. 
Therefore, to apply an impedance-matching approach to a rotor system, care must be used to ensure that the 
response is within a linear range. This requires the establishment of a steady-state, fixed-hub (or motionless hub) 
condition about which small perturbations may be made. Therefore, the loads at the rotor system hub and swash- 
plate interface joints represent the loads due to the operation at the fixed-hub condition plus an incremental load 
due to the motion of the rotor system. Second, because the rotor control system (i.e., pitch links and swashplate) is 
massless in the CAMRAD-II model of the BRR, the swashplate interface joint loads are dependent only upon the 
aerodynamic and dynamic forces in the rotor system. Therefore, the impedance equation must generate forces and 
moments at both the hub and swashplate interface joints based on the fixed-hub condition loads and motion at the 
hub joint. No additional loads caused by motion of the swashplate are introduced. Finally, it must be recognized 
that 4P hub motions generate changes in both 4P and 8P loads because of the changing aerodynamic environment. 
The same may be said of 8P hub motions. This cross coupling must be considered in the development of the 
impedance matrix. 

5.3.2.I. Fixed-hub load conditions. The choice of meaningful fixed-hub, or motionless hub, load conditions is 
an important aspect of the impedance-matching approach because, in general, the developed impedance matrices 
are valid only for a small motion about the fixed-hub condition due to system nonlinearity. For this study, a set of 
flight conditions that are representative of lg steady flight for a helicopter with gross weight of 8000 lb has been 
chosen. The performance parameters required to match these conditions are provided as a function of advance 
ratio p in table 5.1. CAMRAD-II runs are made for each advance ratio; this allows trim procedures internal to 
CAMRAD-II to converge to control positions that satisfy the performance parameters. For each of these runs, the 
rotor is trimmed so that first-harmonic flapping is removed with respect to the rotor shaft, a procedure typical in 
model rotor system analysis. At each of the flight conditions, the resulting control positions are recorded and a 
fixed-hub load vector consisting of the 4- and 8-per-rev hub and swashplate loads is defined. This vector is 
denoted by {F 0 ) and is generated as shown in the following equation: 


1 F S 4 > 


53.2.2 . Hub motion bad conditions. Because rotor systems are inherently nonlinear, no assurance exists that 
an impedance matrix can be developed by using only a cosine or a sine excitation as for fully linear systems. This 
study determined that it is necessary, in general, to consider both cosine and sine excitations in the development of 
an impedance model of a rotor system. A similar approach has been observed in the literature (e.g., ref. 6). This 
requirement is most likely because of the nonlinear effects introduced to the rotor system by the vortex wake; this 
is supported by the rotor impedance matrices developed for the flight speed range in which the rotor system 


Table 5.1. Performance Parameters for Fixed-Hub 
Load Conditions 


p 

C L 

C D 

0.05 

0.0058 

-0.000011 

0.10 

0.0058 

-0.000043 

0.15 

0.0058 

-0.000098 

0.20 

0.0058 

-0.000174 

0.25 

0.0058 

-0.000256 

0.30 

0.0058 

-0.000392 

0.35 

0.0058 

-0.000533 
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behavior is considered to be the most linear 0.15 < p < 0.30. For these flight conditions, the rotor impedance matri- 
ces were found to be more typical in form to those developed for fully linear systems. Additionally, CAMRAD-II 
runs made with a uniform inflow distribution undisturbed by a vortex wake were observed to produce impedance 
matrices typical in form to those for linear systems. Comprehending the nature of this effect would require study 
beyond the scope of this research. For this study, it was necessary to excite a rotor system with both cosine and 
sine hub motions for all flight conditions to develop an appropriate rotor impedance model. 

A series of 24 CAMRAD-II runs is made for each advance ratio to define the load characteristics of the BRR 
under hub motion. These runs include individual cosine and sine excitations for the three translations and the three 
rotations at both 4P and 8 P. For translational hub excitations, an amplitude of 0.01 in. is used, and for rotational 
hub excitations, an amplitude of 0.05° is used. Because pure hub motions are required to develop these loads, 
inputting both translational and rotational motions is necessary in order to obtain pure rotational hub motion with 
the CAMRAD-II model of the BRR. The actuated pylon joint (fig. 4.5) is placed below the hub so that input rota- 
tional motion induces both rotation and translation at the hub. For each of the runs, the control positions generated 
for the fixed-hub conditions were used. In no case did the actuated motion cause the first-harmonic flapping with 
respect to the shaft to go out of trim by more than 0 . 1 °. 

5.3.2.3. Generating impedance matrix. Because the impedance matrix should permit the calculation of inter- 
face loads due to hub motion only, accounting for the fixed-hub loads is necessary when the impedance matrix is 
developed. Mathematically, this is written as 

{F/} = {F 0 } + [C]{X ; } (5.41) 

Therefore, to develop the impedance matrix [C], the equation 

{F V- {F o b } 

t c abl = 7 y \ ( a = 1, 2, ..., 6; b = 1, 2, ..., 6 ) (5.42) 


is used where the subscripts a and b denote matrix element location. 

For clarity, a two-step process is employed to develop the rotor system impedance matrix. First, matrices are 
developed based on harmonic response due to harmonic excitation. These matrices are shown in the following 
equations: 


[c 4/4 ] = 

[C 8/4 ] = 


r 



' 

[C„ 4/4 ] [0] 

c 

C„ [0] 

" 4/8 


.1%,) [°! 

W/8 “ 

t°i 


1 C *,J <°i 

f _ 

’ c »„„ i°i 


JSJ [°1 

^ 8/8 “ 

(0 ! 



(5.43) 


The subscripts denote the interface joint for the response (H for hub and S for swashplate), the response harmonic 
to the left of the slash, and the excitation harmonic to the right of the slash. The matrices [0] are used because the 
impedance load at the interface joints depends on hub joint motion but not swashplate joint motion. Second, the 
final assembly of the system impedance matrix occurs as follows: 


[C 4/4 ] [C 4/8 ] 

[*^8/4] [c 8/8 | 


(5.44) 
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A set of the component matrices required to develop the impedance matrix for p. - 0.10, a highly nonlinear 
flight condition, is shown in the following equations: 

For 4P hub interface joint loads due to 4P hub actuation. 


I 0.8256 

0.0125 

- 0.0013 

- 0.0246 

0.0114 

- 0.0073 

- 0.0087 

0 . 0204 

- 0.0785 

- 0.0176 

- 0.0993 

0 

•- 0.0102 

0.8286 

0.0202 

0.0037 

0.0178 

0.0024 

0.0375 

- 0.0282 

- 0.0376 

- 0.0751 

- 0.1081 

0 

- 0.0046 

0.0253 

0.8226 

0.0106 

0.0097 

0.0017 

0.0702 

0.0055 

- 0.0448 

- 0.0176 

- 0.0757 

0 

- 0.0238 

- 0.0028 

- 0.0129 

0.8262 

0.0127 

0.0100 

0.0212 

0.0545 

- 0.0139 

- 0.0288 

- 0.0967 

-0 

- 0.0011 

0.0022 

- 0.0014 

- 0.0016 

0.7547 

0.0335 

0.0280 

0.0207 

- 0.0193 

- 0.0262 

- 0.0803 

-0 

- 0.0015 

0.0008 

- 0.0019 

0.0002 

- 0.0203 

0.7363 

- 0.0008 

- 0.0048 

- 0.0605 

- 0.0221 

0.0345 

0 

- 0.0377 

- 0.0310 

- 0.0398 

- 0.0118 

0.0652 

- 0.0781 

4.7546 

1.1236 

0.0286 

- 0.9574 

- 0.5711 

0 

0 . 0259 

- 0.0198 

- 0.0031 

- 0.0337 

0.1344 

- 0.0843 

- 0.5825 

4.6405 

0 . 6463 

- 0.4314 

- 0.8440 

0 

- 0.0501 

0.0509 

- 0.1226 

- 0.1342 

0.2925 

0.1655 

- 0.1202 

1 .2477 

3.5257 

0.0922 

- 2.8765 

0 

- 0.0037 

0 . 0081 

0.0302 

- 0.0526 

0.2203 

0.0835 

- 0.5749 

- 0.0619 

- 0.8938 

4.2094 

- 1.5096 

-0 

- 0.0288 

- 0.0112 

- 0.0342 

0.0002 

0.4766 

- 0.1810 

- 0.0143 

0.3993 

0.2083 

- 0.0580 

- 0.1922 

6 

0.0035 

- 0.0055 

- 0.0127 

- 0.0312 

0.4258 

0.4290 

- 0.1431 

- 0.2952 

- 0.4960 

- 0.0067 

- 8.1871 

0 


0.0703 

0.0144 

0.0274 

0.0243 

0.0848 

0.0518 

0.6210 

0.5540 

0.2625 

0.0426 

6.7245 

0.5865 


For 4P hub interface joint loads due to 8P hub actuation, 


(5.45a) 


[C ] = 10 4 x 


0.0001 - 0.0021 - 0.0005 0.0009 

- 0.0012 0.0009 - 0.0001 0.0016 

- 0.0004 - 0.0031 - 0.0018 - 0.0031 
0.0000 - 0.0016 0.0010 0.0026 
0.0123 0.0137 0.0178 0.0163 

0.0081 0.0144 0.0001 0.0062 

- 0.0002 0.0328 0.0094 - 0.5996 

- 0.0008 0.0017 0.0148 0.1827 

- 0.0061 - 0.0342 - 0.0050 - 1.3330 
- 0.0018 - 0.0342 - 0.0011 - 0.2651 
- 0.0433 - 0.1293 - 0.0736 3.7321 

- 0.0225 - 0.0490 0.0145 0.4054 


0.0003 - 0.0001 
- 0.0003 0.0000 


0.0025 
- 0.0004 
0.0083 
0.0004 
0 . 0089 


0.0001 

0.0065 

0.0228 

0.0076 

0.0690 


- 0.0237 - 0.0190 
- 0.0319 - 0.0337 
0.0178 - 0.0391 
- 0.0290 - 0.1371 
0.0243 - 0.0209 


0.0017 - 0.0073 
0.0036 - 0.0176 
- 0.0120 - 0.0072 
0.0147 - 0.0138 
- 0.0216 0.0508 

- 0.0580 0.0743 

0.0967 - 0.1441 
0.0333 0.1566 

- 0.0673 0.0509 

0.0763 0.1101 

0.1951 - 0.3260 
0.3279 - 0.4412 


0.0213 - 0.0067 
0.0048 0.0216 
0.1349 0.1138 
0.0399 - 0.0009 
0.3468 0.4883 
0.1226 - 0.1036 
- 0.1552 - 0.0565 
- 0.4218 - 0.3658 
- 0.8648 - 0.6110 
0.0135 0.5080 


0.0021 
- 0.0048 
- 0.0060 
0.0067 
0.1316 
- 0.0408 
0 . 1241 
- 0.0639 


- 0.0052 

0.0128 

- 0.0248 

0.0022 

0.1348 

0.0519 

0.3411 

0.0890 


- 0.2739 - 0.2546 
0.0437 - 0.4112 
- 0.5392 - 0.9326 


For 8P hub interface joint loads due to 4P hub actuation, 


0.0026 0.0041 0.0026 0.0002 0.0065 

0.0031 0.0035 0.0041 0.0040 0.0061 

- 0.0059 - 0.0068 - 0.0072 - 0.0047 - 0.0067 
- 0.0001 - 0.0005 - 0.0010 - 0.0047 - 0.0041 
0.0089 0.0056 0.0067 - 0.0081 0.0021 

0.0421 0.0353 0.0361 0.0263 0.0658 

- 0.0190 - 0.0275 - 0.0222 - 0.0099 - 0.0319 
0.0182 0.0201 0.0046 0.0041 0.0254 

- 0.0055 0.0020 0.0090 0.0205 0.0029 

- 0.0316 - 0.0295 - 0.0300 - 0.0100 - 0.0455 
- 0.1271 - 0.0838 - 0.0992 - 0.0847 - 0.1850 
0.1380 0.1093 0.1083 0.0473 0.1829 


0.0035 

0.0036 

- 0.0042 


0.0180 

0.0464 

- 0.0494 


0.0035 - 0.0285 
0.0129 0.1674 

0.0224 0.5152 

- 0.0152 - 0.1804 
0.0081 - 0.0127 
0.0179 0.1391 

- 0.0468 - 0.3246 
- 0.0398 - 1.4289 
0.1200 1.6127 


0.0582 0.0600 0.0303 

0.0447 0.0952 0.0823 

- 0.0587 - 0.1429 - 0.1259 
0.0238 - 0.0346 - 0.0195 
0.1394 0.0401 0.0552 

0.6610 0.7478 0.8181 

- 0.1851 - 0.3536 - 0.3439 
0.3161 0.3670 0.1848 

- 0.0041 0.5911 0.1475 

- 0.2295 - 0.7601 - 0.4087 
- 1.4209 - 2.1158 - 2.2583 
2.0162 2.3352 2.6169 


0.0019 

0.0438 


0.0297 

0.0204 


- 0.0532 - 0.0539 
- 0.0543 - 0.0140 
- 0.1283 0.0938 
0.2448 0.2691 
0.1058 0.0397 
0.0291 0.4987 
0.3062 - 0.0986 
0.1365 - 0.2660 
- 1.1741 - 0.8416 
0.5789 1.1202 


(5.45c) 
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For 8P hub interface joint loads due to 8P hub actuation, 


[C„ ] = 10 3 x 

w 8/8 


0.3497 

- 0.0107 

0.0051 

0.0005 

0.0004 

0.0008 

0.0020 

0.0027 

0.0291 

- 0.0054 

0.0066 

0.0029 

0.0113 

0.3502 

- 0.0070 

0.0154 

0.0150 

0.0010 

- 0.0046 

0.0052 

0.0164 

0.0300 

0.0065 

0.0057 

- 0.0053 

- 0.0084 

0.3488 

- 0.0264 

- 0.0251 

- 0.0031 

- 0.0231 

0.0041 

- 0.0109 

- 0.0009 

- 0.0112 

- 0.0099 

0.0072 

- 0.0052 

0.0108 

0.0073 

0.0066 

- 0.0010 

- 0.0068 

- 0.0287 

- 0.0057 

0.0005 

- 0.0022 

- 0.0041 

- 0.0007 

- 0.0005 

- 0.0004 

0.0559 

0.0540 

- 0.0410 

0.0082 

- 0.0129 

0.0012 

0.0056 

- 0.0436 

- 0.1086 

0.0034 

0.0061 

0.0049 

0.0050 

0.0036 

0.3381 

- 0.0233 

0.0084 

0.0747 

0.0503 

0.1639 

0.0015 

0.0025 

- 0.0052 

- 0.0055 

0.0079 

- 0.0040 

0.0004 

2.1830 

- 0.0639 

- 0.0627 

- 0.0535 

- 0.0332 

- 0.0610 

0.0039 

0.0065 

- 0.0008 

0.2161 

- 0.0013 

0.0001 

0.0557 

2.1796 

0.0532 

- 0.0209 

0.0467 

- 0.0032 

0.0058 

- 0.0019 

0.0032 

- 0.3078 

- 0.0007 

- 0.0002 

0.0233 

0.0659 

2.1922 

- 0.0437 

0.0088 

0.0295 

0.0011 

0.0008 

- 0.0026 

0.0968 

0.0006 

- 0.0048 

- 0.0323 

0.0179 

0.0009 

2.1769 

- 0.0492 

- 0.0363 

- 0.0136 

- 0.0186 

- 0.0139 

0.6445 

- 0.0649 

- 0.1811 

0.1166 

- 0.0568 

- 0.2078 

- 0.1365 

6.1916 

- 3.2330 

0.0057 

0.0144 

0.0120 

0.2088 

0.1645 

- 0.0453 

- 0.0409 

0.0240 

0.2184 

0.1697 

3.2380 

6.4803 


For 4P swashplate interface joint loads due to 4P hub actuation, 


(5.45d) 


[C„ ] = l(r x 

A 4/4 


- 0.0138 

0.0021 

- 0.0051 

- 0.0018 

- 0.0075 

- 0.0019 

- 0.0543 

0.1119 

0.1151 

- 0.1799 

0.0212 

- 0.0745 

- 0.0060 

- 0.0093 

- 0.0010 

- 0.0047 

- 0.0032 

0.0046 

- 0.2261 

- 0.0119 

0.0514 

0.1905 

0.0058 

0.0565 

0.0010 

- 0.0033 

- 0.0036 

0.0006 

- 0.0035 

0.0017 

- 0.3062 

0.0576 

- 0.0228 

0.1821 

0.0442 

0.0526 

0.0089 

0.0079 

0.0020 

- 0.0019 

0.0073 

0.0045 

- 0.0474 

- 0.1903 

- 0.0988 

0.0913 

- 0.0157 

0.1211 

C . 01 94 

- 0.0190 

0.0300 

- 0.0326 

- 0.0696 

- 0.0342 

- 0.0882 

0.1642 

- 0.2141 

- 0.0552 

- 3.6747 

- 1.4332 

0 . 0254 

0.0325 

0.0304 

0.0532 

0.0583 

- 0.0568 

- 0.0284 

0.0396 

0.0575 

0.5948 

1.7857 

- 3.0352 

0.2049 

- 0.0415 

0.0668 

0.0273 

0.0981 

0.0300 

- 0.0712 

- 1.4626 

- 1.7935 

3 . 1471 

- 0.3432 

1.0809 

0.1129 

0.1581 

0.0204 

0.0663 

0.0707 

- 0.0579 

3.1639 

- 0.3705 

- 1.0518 

- 2.5459 

0.0485 

- 0.6286 

0.0236 

0.0437 

0.0662 

- 0.0022 

0.0759 

- 0.0184 

4.6688 

- 1.1728 

0.0187 

- 2.1509 

- 0 . 4502 

- 0.5214 

- 0.1143 

- 0.0912 

- 0.0286 

0.0388 

- 0.1001 

- 0.0745 

1.3385 

2.8577 

1.3036 

- 1 . 8980 

0.1669 

- 1.7056 

0.0048 

- 0.0051 

0.0084 

- 0.0093 

- 0.0197 

- 0.0098 

- 0.0085 

0.0476 

- 0.0538 

- 0.0340 

- 1.0384 

- 0.4082 

0.0064 

0.0083 

0.0084 

0.0149 

0.0159 

- 0.0161 

- 0.0150 

0.0221 

0.0248 

0.1705 

0.5031 

- 0.8615 


For 4P swashplate interface joint loads due to 8P hub actuation. 


(5.45e) 


[C_ ] = lCr x 


- 0.0082 

- 0.0079 

- 0.0032 

- 0.0042 

- 0.0033 

- 0.0069 

0.0169 

- 0.0245 

- 0.1007 

- 0.0731 

- 0.1062 

- 0.0763 

- 0.0016 

- 0.0036 

0.0015 

0.0005 

0.0017 

- 0.0043 

0.0225 

- 0.0337 

- 0.0182 

- 0.0499 

0.0066 

- 0.0812 

- 0.0014 

- 0.0049 

0.0014 

0.0029 

0.0027 

- 0.0027 

0.0174 

- 0.0378 

- 0.0272 

- 0.0183 

- 0.0081 

- 0.1002 

0 . 0057 

0.0070 

0.0010 

0.0061 

0.0041 

0.0075 

0.0053 

- 0.0121 

0.0850 

0.0699 

0.1088 

0.0464 

- 0.0477 

- 0.0106 

- 0.0460 

- 0.0313 

- 0.0075 

- 0.0234 

0.0949 

- 0.5579 

- 0.2422 

0.1665 

- 0.1210 

- 0.1945 

- 0.0140 

0.0367 

0.0490 

- 0.0029 

0.0168 

0.0613 

- 0.0048 

0.1159 

0.2004 

0.1439 

0.2055 

0.1151 

0.1124 

0.1027 

0.0498 

0.0672 

0.0554 

0.0927 

- 0.1926 

0.2197 

1.3875 

1.0121 

1.5158 

0.8231 

0.0400 

0.0732 

- 0.0167 

0.0103 

- 0.0129 

0.0873 

- 0.3101 

0.4783 

0.5247 

0.9226 

0.2210 

1.3114 

0.0452 

0.0939 

- 0.0085 

- 0.0289 

- 0.0288 

0.0603 

- 0.3016 

0.6659 

0.6928 

0.4599 

0.4292 

1.6796 

- 0.0771 

- 0.0944 

- 0.0199 

- 0.0908 

- 0.0648 

- 0.0997 

- 0.1600 

0.2726 

- 1.1909 

- 0.8928 

- 1.6105 

- 0.4436 

- 0.0137 

- 0.0030 

- 0.0132 

- 0.0091 

- 0.0024 

- 0.0067 

0.0264 

- 0.1559 

- 0.0706 

0.0450 

- 0.0379 

- 0.0520 

- 0.0044 

0.0098 

0.0138 

- 0.0011 

0.0045 

0.0167 

- 0.0006 

0.0319 

0.0505 

0.0343 

0.0518 

0.0265 


(5.45f) 
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For 8P swashplate interface joint loads due to 4P hub actuation, 


[C s ] = 10^ 

■* 8/4 


- 0.0055 

- 0.0076 

- 0.0038 

- 0.0004 

- 0.0095 

- 0.0018 

- 0.0046 

- 0.0594 

- 0.0860 

- 0.0755 

0.0408 

- 0.1493 

0.0024 

- 0.0045 

- 0.0042 

- 0.0038 

- 0.0035 

- 0.0062 

- 0.0057 

0.0173 

- 0.0468 

- 0.0387 

0.0865 

0.0350 

- 0.0007 

0.0063 

0.0055 

0.0031 

0.0038 

0.0060 

0.0092 

- 0.0061 

0.0457 

0.0559 

- 0.0891 

- 0.0549 

- 0.0041 

- 0.0045 

- 0.0013 

0.0019 

- 0.0058 

- 0.0013 

0.0197 

- 0.0425 

- 0.0630 

- 0.0448 

0.0827 

- 0.1365 

- 0.0267 

- 0.0297 

- 0.0412 

- 0.0098 

- 0.0306 

- 0.0216 

- 0.3239 

- 0.4213 

- 0.3923 

- 0.4334 

0.2449 

- 0.4621 

- 0.0183 

- 0.0274 

- 0.0274 

- 0.0437 

- 0.0380 

- 0.0147 

- 0.2268 

- 0.3427 

- 0.4945 

- 0.5183 

0.1661 

0.0989 

0.0733 

0.1213 

0.0655 

0.0108 

0.1410 

0 . 0409 

0.0316 

0.7768 

1.2959 

1.2321 

- 0.8264 

1.9072 

- 0.0458 

0.0513 

0.0559 

0.0568 

0.0305 

0.0851 

0.1220 

- 0.4033 

0.3934 

0.3921 

- 0.9636 

- 0 . 8724 

0.0302 

- 0.0635 

- 0.0604 

- 0.0423 

- 0.0209 

- 0.0773 

- 0.0652 

0.2957 

- 0.3321 

- 0.4998 

1.0959 

1.3223 

0.0560 

0.0829 

0.0382 

- 0.0086 

0.1003 

0.0397 

- 0 . 1789 

0 . 6268 

1 . 1259 

0.8686 

- 1.4917 

1.8268 

- 0.0077 

- 0 . 0102 

- 0.0130 

- 0.0034 

- 0.0100 

- 0.0075 

- 0.0936 

- 0.1219 

- 0 . 1252 

- 0.1384 

0.0910 

- 0.1275 

- 0.0041 

- 0.0070 

- 0.0077 

- 0.0129 

- 0.0097 

- 0.0043 

- 0.0673 

- 0.0866 

- 0.1286 

- 0.1386 

0.0346 

0.0595 


For 8P swashplate interface joint loads due to 8P hub actuation. 


(5.45g) 


[C. ] = 10 x 

** 8/8 


0 .0083 
■ 0.0029 
0.0019 
0.0000 
0.0247 
0.1782 
0.1261 
0.0516 
0.0106 
- 0.0047 
0.0066 
0.0501 


- 0.0013 
- 0.0005 
0.0009 
- 0.0008 
- 0.0032 
- 0.0031 
0.0210 
0.0044 
- 0.0082 
0.0142 
- 0 . 0012 
- 0.0007 


0.0000 

0.0002 

0.0000 

0.0000 

- 0.0020 

- 0.0051 

- 0.0008 

- 0.0030 

0.0006 

0.0005 

- 0.0006 

- 0.0014 


■ 0.0007 

0.0002 

0.0000 

- 0.0006 

0.0010 

- 0.0024 

0.0108 

- 0.0049 

0.0015 

0.0077 

0.0002 

- 0.0006 


- 0.0002 

- 0.0003 

0.0009 

0.0006 

- 0.0217 

- 0.0209 

0.0028 

0.0052 

- 0.0087 

- 0.0035 

- 0.0063 

- 0.0060 


- 0.0004 

- 0.0002 

- 0.0004 

0.0003 

0.0156 

- 0.0270 

0.0049 

- 0.0002 

0.0038 

0.0005 

0.0045 

- 0.0077 


0.0084 

- 0.0168 

0.0191 

- 0.0003 

0.0067 

0.0019 

- 0.0682 

0.2407 

- 0.2735 

0.0536 

- 0.0013 

- 0.0004 


0.0167 

0.0060 

0.0022 

0.0193 

0.0088 

- 0.0067 

- 0.2340 

- 0.0351 

- 0.0817 

- 0.2678 

0.0030 

- 0.0052 


- 0.0263 

- 0.0016 

0.0102 

- 0.0211 

- 0.0366 

- 0.0622 

0.3730 

- 0.0443 

- 0.0699 

0.3212 

- 0.0140 

- 0.0131 


- 0.0029 
- 0.0230 
0.0212 
0.0067 
- 0.0109 
- 0 . 0408 
0 .1005 
0.3184 
- 0.2952 
- 0.0260 
- 0.0077 
- 0 . 0141 


- 0.0065 

- 0.0039 

0.0021 

- 0.0029 

- 0.0578 

- 0.0896 

0.0891 

0.0405 

- 0.0036 

0.0714 

- 0.0173 

- 0.0247 


0.0989 
- 0.0477 
- 0.0418 
- 0 . 0049 
- 0.1989 
2.0997 
- 1.5497 
0.7920 
0.3968 
0.0510 
- 0.0507 
0 . 5904 


(5.45h) 

As shown, the form of the matrices are not indicative of a linear system. The reader is referred to equation (5.45a), 
which presents the impedance matrix component for the 4P hub interface joint response due to 4P actuation. The 
submatrix labeled A shows that, unlike the mobility and motion matrices, the diagonal elements of the submatrices 
are unequal for the impedance matrix. Similarly, the off-diagonal submatrix terms are not equal and opposite as 
for the mobility and motion matrices. 


53.2.4. System linearity. Because the impedance matrices generated for the BRR are not typical of those for 
fully linear systems, immediate concern is raised as to the applicability of the impedance-matching approach to the 
rotor system. Therefore, a study has been performed which validates the impedance matrices for translational 4P 
hub motions with amplitudes of 0.015 in. and rotational 4P hub motions with amplitudes of 0.10°. Figure 5.3 
shows representative results for 4P axial hub motions ranging from —0.015 to 0.015 in. As shown, the axial hub 
load response for the range of motion shown remains linear with respect to the imposed hub motions. This type of 
response indicates that the impedance-matching method is applicable for the range of motions expected from this 
study. 


5.3.3. System Coupling 

The final step necessary to couple the BRR system to the ARES-II is the generation and application of the sys- 
tem coupling equations. These equations are generated in a manner similar to that used for the sample model (sec- 
tion 5.1.4). The motion of the interface joints is most easily defined by observing the system from the point of 
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Amplitude of 4P axial cosine hub motion, in. 



Amplitude of 4P axial sine hub motion, in. 


Figure 5.3. BRR model linearity for axial cosine and sine hub motion excitations. 


view of ARES-IL Utilizing the mobility equation (eq. (5.35)), the motion equation (eq. (5.39)), and the linear the- 
ory of superposition, the motion of the interface joints on the ARES-II may be expressed as 

{X ; } = [A]{F / } + [B]{U} (5.46) 

Load compatibility is enforced at the interface joints by substituting the rotor impedance equation (eq. (5.41)) into 
equation (5.46). Displacement compatibility is then enforced by solving for the interface joint displacements {X/} 
as follows: 


{X 7 } = [[I]-[A][cf 


[A]{F 0 } + [B]{U} 


(5.47) 


The interface joint forces {F 7 } are determined by substituting equation (5.47) into equation (5.41) which gives 


{F 7 } = {F 0 } + [C][[I] - [A][C]J j[A]{F 0 } + [B]{U}| 


(5.48) 
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5.3.4. System Loads 

The development of system loads throughout the ARES-II testbed is of interest. Of particular importance for 
this study are the loads generated at each of the system balances. Therefore, a general set of equations has been 
developed to define the loads in the ARES-II. The approach taken has been to generate two sets of matrices. The 
first set [D] defines the loads transfer from the interface joints to locations throughout the ARES-II, and the second 
set [E] defines the effect of actuation on ARES-II component loads. 

The load transfer matrices are developed in a manner identical to that used to develop the system mobility 
matrices in section 5. 3. 1.2. The only difference is that a load response due to load excitation is developed instead 
of a motion response due to load excitation. The DADS runs used to develop the mobility matrix are used to define 
the load response; hence, no additional runs are required to develop the matrix [D], The matrices defining the 
effects of actuation on system loads are developed in a manner identical to the motion matrix development 
described in section 5. 3. 1.3. Again, the DADS runs that have already been made supply all information necessary 
to develop the matrix [E], The loads vector {L} may then be established by using 


{L} = [D]{F 7 } + [E]{U} 


(5.49) 


6. Presentation of Results 

Results from independent, uncoupled analysis of the ARES-II and the BRR models are presented in 
sections 6.1 and 6.2; coupled analysis results are presented in section 6.3. Because of the volume of results pre- 
sented, tables and figures for section 6 are placed at the end of the section rather than integrated in the text. 

6.1. ARES-II Motion Actuation 

The actuator motion analysis of section 4.1.2 was used to develop time- varying actuator length schedules for 
all ARES-II actuation in DADS. For all results presented, a simulated center of gravity was chosen which was 
placed on the rotor shaft centerline and 10.0 in. above the fixed-system balance center. Figure 6.1 presents a repre- 
sentative actuator length schedule for 0.01 in. 4P axial sine motion of an inflexible ARES-II model by using the 
equations developed in section 4.1.2. However, application of motions to the DADS elastic model of the ARES-II 
produces a different response. Table 6.1 presents the response at the hub interface joint using actuator schedules 
that command 0.01 in. 4P sine translation and 0.05° 4P sine rotation. The columns of table 6.1 represent the 4P 
sine motion commanded. The rows represent the cosine and sine components of the hub interface joint response of 
the elastic model due to the commanded motion. As shown by the table, the desired motion is either attenuated or 
amplified, phase shift is evident, and some significant coupling of pylon motion exists. 

When the elastic model of the ARES-II is considered, the Stewart platform must account for the dynamic 
response of the system to produce the desired motion. The motion matrix for the hub [ B } from equation (5.36a) 

may be used to calculate a control vector that provides a specified hub interface joint motion as follows: 

{U} = [B^r’lX^} (6.1) 


Figure 6.2 presents the actuator schedule required to obtain a pure 0.01 in. 4P axial (x direction) translational 
motion of the rotor hub interface joint when flexibility is included in the model. A comparison of figure 6.2 with 
figure 6.1 indicates the difference in actuator motion required to obtain pure hub motions for the flexible model. 
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6.2. BRR Fixed-Hub Loads 

A fixed-hub load vector {F 0 } was obtained with the CAMRAD-II model of the BRR for each forward flight 
condition listed in table 5.1. The results for the 4P components of the load are presented in figure 6.3 along with 
the total magnitude of the load transmitted to the fixed system by the hub and swashplate interface joint loads. For 
this figure, the total magnitude has been calculated with 


F «ut - i< f H t + 'X > 2 + < f H, + 'V 2 


This formulation neglects the effect of lateral swashplate loads on the total hub moment; however, these loads are 
small for all cases. 

The trends observed for total magnitude as a function of flight speed are typical of those encountered with all 
rotor systems. Typically, low-speed flight (p = 0.05) vibrations are generally minimal but increase rapidly through 
p = 0.10 because of a phenomenon called transition in which many blade- vortex interactions occur. Vibratory 
loads through the cruising range (p = 0.20 to 0.25) are generally low, but loads increase again at the higher flight 
speeds (p > 0.30). Of particular note is the magnitude of the fixed-system swashplate loads at the high speeds 
(p > 0.30) for normal force and pitching and rolling moments. If these loads were neglected as has been common 
practice throughout much of the literature, an inaccurate assessment (up to 30 percent error) of the total fixed- 
system loads would result. Additionally, for the ARES-II configuration in which the loads sensed by both the 
rotating- and fixed-system balances are important; neglecting the swashplate loads would result in a serious defi- 
ciency in the balance loads prediction. 

6.3. BRR/ARES-II Coupled Results 

Because the system coupling method developed in section 5.3.3 requires only simple matrix operations, cou- 
pled rotor-testbed solutions for a wide range of ARES-II motions may be obtained in a very short time. Because of 
this, a large volume of results have been obtained from this study. Selected results have been chosen for presenta- 
tion in this section. Tables 6.2 and 6.3 present representative interface joint motions and loads and allow the com- 
parison of the uncoupled rotor loads with the coupled rotor loads. Tables 6.4 and 6.5 present the loads at the fixed- 
system and rotating-system balances. 

The notation used throughout the tables is straightforward. Displacements are identified as translational (or 
linear) motions along an axis (L x , L Y , L^, or as rotational (or angular) motions about an axis (A x , A y , A z ). 
Interface joint locations are identified with either H for the hub or S for the swashplate. Amplitudes are identified 
by their frequency (4P or 8P) and the component (cosine or sine). For example, H8S refers to an 8P sine compo- 
nent amplitude at the hub interface joint. Forces and moments are presented in units of pounds and inch-pounds. 
Displacements are presented in units of inches and degrees. 

The rotor-testbed coupled results in tables 6.2 and 6.3 are obtained with equations (5.47) and (5.48). Table 6.2 
presents the motion and loads for the three forward flight speeds: (X = 0.10 (transition), p = 0.20 (cruise), and 
p = 0.30 (high speed). No pylon actuation was used. Table 6.3 presents the motion and loads for representative 
pylon actuation schedules at a low-speed advance ratio p = 0.10. Results are included for 0.01 in. 4P sine actuation 
of the three translational motions (tables 6.3(b) through (d)), and for 0.05° 4P sine actuation of the pitch motion 
(table 6.3(e)). Also presented are results for 0.01 in. 4P cosine actuation in the axial direction (table 6.3(a)) to dem- 
onstrate representative differences in results due to cosine versus sine actuation. 

The balance loads results for the coupled system are obtained with equation (5.49). Table 6.4 presents the bal- 
ance loads for p = 0.10 without pylon actuation. Therefore, these results may be compared with those presented in 
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table 6.2(a). Table 6.5 presents the balance loads for the five pylon actuation schedules used to generate the results 
in table 6.3. 

The fidelity of the balance measurements for normal force and pitch and roll moments are presented in 
figure 6.4 for the coupled system without actuated motion. In this figure, the balance moments shown are 
corrected for the moments produced by the interface joint lateral forces and presented throughout the forward 
flight speed range. The correction equations are 


M = M +4.213 F„ 

RB con RB y 


Af =Af -4.213F„ 

yRB y RB "x 

con 

Af = M r +23.538F„ + 15.318F. 

X FB X FB “ V v 


M , =M -23.538F„ - 15.318F V 

yFB ^FB .t x 


(6.2) 


Table 6.6 presents the results obtained by using the control vector inputs necessary to obtain 0.01 in. 4P axial 
sine motion at the hub interface joint for the coupled rotor-testbed. The control inputs were established by using an 
iterative solution technique to minimize the phase shift and coupling evident in the results presented in 
table 6.3(b). No closed-form method of calculating the proper control inputs as for the stand-alone ARES-II 
(section 6.1) is possible. As shown, significant changes in control inputs are necessary to achieve proper hub 
motion, particularly for the sine components of the axial and pitch actuation. 

Several observations regarding the results obtained for the coupled system are offered as follows: 

1. Although the loads due to coupling (i.e., the impedance forces and moments) [C]{X/} are generally small for 
the unactuated cases (table 6.2), they are significant enough that considering their effects is useful, particu- 
larly for the axial and side loads where the loads due to coupling can be as large as or larger than the uncou- 
pled loads {F 0 }. A case could be made for neglecting the loads due to coupling for the normal force and the 
moments because their contributions are small for the unactuated results. 

2. When the pylon is actuated, including the loads due to coupling becomes critical in the analysis. Some of the 
load may be attributed to the dynamic load associated with actuating the rotor system. However, a significant 
portion of the load is due to the change in the aerodynamic response of the rotor. 

3. The motion response is greatly affected by the rotor system for the actuated cases. Comparison of the results 
of the testbed alone (table 6.1) with the results for the coupled rotor-testbed (table 6.3) indicates that signifi- 
cantly different hub motions may be expected when the rotor is coupled to the ARES-II. All motions have 
been observed to grow for the coupled system. For the axial motion, the growth is fortuitous because a 
response close to that requested by the control vector results. This response is, however, dependent upon the 
rotor system and flight parameters selected. 

4. The motion response is typically amplified compared with the input control vector {U} for the actuated 
cases. The exception to this is for side motion in which the motion is attenuated. Some motions can be signif- 
icantly larger than expected based on the input control vector. For example, the displacements associated 
with pitch actuation (table 6.3(e)) are nearly three times larger than the input. Additionally, the axial response 
(L x ) for these conditions exceeds twice the range established for the linear response range of the rotor 
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system. To ensure validity of the results, more CAMRAD-H runs would have to be made to determine the 
boundaries of the linear range. Furthermore, as indicated in table 6.5(e), the pitching moment at the fixed- 
system balance is excessive for the actuated pitch motion. These loads are in excess of the balance limits and 
would damage the balance. 


5. Because of the observations in paragraph 4 and the phase shift and coupled rotor pylon motions evident 
throughout the results, it is apparent that a closed-loop control system must be implemented to account for 
the dynamic response of the system and provide the desired motion at the rotor hub. A comparison of 
table 6.6, in which the inputs are chosen to provide the desired hub interface joint motion, and table 6.3(b), in 
which the desired motions are input for the control vector, illustrates the importance of a closed-loop control 
system for proper performance of the ARES-II. 


6. The measurement potential of the rotating balance loads is much better than that of the fixed-system balance. 
Significant loads due to the dynamic response of the ARES-II are apparent throughout the results of the 
fixed-system balance loads. The rotating balance shows minimal effect due to dynamic response. Therefore, 
if a reliable means of measuring the swashplate interface loads during wind tunnel testing is identified (e.g., 
strain-gauge pitch links), better system loads may be obtained by using the rotating balance rather than the 
fixed-system balance 


Table 6. 1 . Hub Interface Joint Response Due to Actuated 4P Sine Motion 


[No rotor model] 


Type of hub 
joint response 

Commanded 4P sine motion for — 

0.01 in. 
axial 

0.01 in. 
side 

0.01 in. 
normal 

0.05° 

roll 

0.05° 

pitch 

0.05° 

yaw 

Axial cos 

-0.0008 

0.0000 

-0.0010 

0.0000 

-0.0008 

0.0000 

Axial sin 

0.0069 

0.0000 

0.0054 

0.0000 

0.0233 

0.0000 

Side cos 

0.0000 

-0.0006 

0.0000 

-0.0001 

0.0000 

0.0003 

Side sin 

0.0001 

0.0062 

0.0001 

-0.0156 

0.0001 

0.0019 

Normal cos 

-0.0001 

0.0000 

-0.0011 

0.0000 

-0.0006 

0.0000 

Normal sin 

0.0000 

0.0000 

0.0103 

0.0000 

0.0006 

0.0000 

Roll cos 

0.0000 

0.0011 

0.0000 

0.0002 

0.0000 

-0.0005 

Roll sin 

-0.0001 

-0.0048 

-0.0001 

0.0545 

-0.0001 

0.0024 

Pitch cos 

-0.0017 

0.0000 

-0.0026 

0.0000 

-0.0022 

0.0000 

Pitch sin 

0.0066 

0.0000 

0.0154 

0.0000 

0.0766 

0.0000 

Yaw cos 

0.0000 

0.0004 

0.0000 

0.0001 

0.0000 

-0.0004 

Yaw sin 

0.0000 

-0.0064 

0.0000 

-0.0020 

0.0000 

0.0666 
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Table 6.2. Coupled Rotor-Testbed Results for No Actuated Motion 
(a) p = 0.10 


Component 

L x 

Ly 

L z 


Ay 

Az 

Input control vector { U } 

4C 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

4S 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

Hub and swashplate interface displacements { X 7 } 

H4C 

-0.000172 

-0.000594 

0.000003 

0.004305 

-0.000708 

-0.003719 

H4S 

0.000866 

-0.000251 

0.000092 

0.001885 

0.005145 

0.001361 

S4C 

-0.000082 

-0.000140 

-0.000004 

0.001663 

-0.000437 

-0.000990 

S4S 

0.000296 

-0.000080 

0.000089 

0.000391 

0.002483 

0.000365 

H8C 

0.000000 

0.000055 

0.000029 

-0.000380 

-0.000032 

0.000074 

H8S 

0.000002 

-0.000187 

-0.000064 

0.000186 

-0.000077 

-0.001580 

S8C 

0.000001 

0.000021 

0.000030 

-0.000060 

0.000035 

0.000556 

S8S 

0.000009 

-0.000154 

-0.000065 

0.000294 

-0.000023 

-0.001419 


Component 

F x 

Fy 

F Z 

M x 

My 

M z 

Hub and swashplate interface forces and moments {F /} 

H4C 


-15.74 

19.28 

62.01 

-2.34 

-132.93 

H4S 


0.76 

8.36 

71.31 

67.46 

52.75 

S4C 

0.42 

-0.11 

-0.15 

-6.38 

0.52 

-0.02 

S4S 

-0.07 

-0.04 

7.16 

1.34 

0.26 

2.02 

H8C 

-0.26 

1.98 

-2.58 

-7.05 

-2.27 

-24.14 

H8S 

0.75 

0.07 

5.01 

-5.13 

-5.59 

-6.57 

S8C 

-0.04 

-0.01 

0.17 

0.56 

0.24 

0.05 

S8S 

0.04 

0.00 

-0.95 

-0.61 

-0.03 

-0.27 

Hub and swashplate forces and moments for fixed-hub condition {F 0 } 

H4C 


- 11.11 

19.25 

58.31 

-3.71 

-122.66 

H4S 

7.22 

2.89 

7.52 

70.28 

64.34 

57.23 

S4C 


-0.10 


—6.53 



S4S 


-0.04 

7.32 

1.31 

0.33 


H8C 


-0.31 

-2.58 

—4.98 

-7.73 


H8S 


0.09 


-2.23 

-2.84 


S8C 


-0.04 


-0.12 

0.44 


S8S 


0.02 


-0.46 

-0.25 


Impedance forces and moments [C] {X/} 

H4C 

-1.42 

—4.64 

0.03 

3.70 

1.37 

-10.27 

H4S 

6.98 

-2.13 

0.84 

1.03 

3.12 

—4.49 

S4C 

-0.01 

0.00 

0.18 

0.15 

0.09 

0.05 

S4S 

0.00 

0.01 

-0.15 

0.03 

-0.07 

-0.04 

H8C 

0.08 

2.29 

0.00 

-2.07 

5.46 

—4.00 

H8S 

-0.15 

-0.02 

-1.06 

-2.90 

-2.75 

-16.63 

S8C 

-0.04 

0.02 

-0.06 

0.68 

-0.19 

-0.02 

S8S 

0.01 

-0.01 

-0.64 

-0.15 

0.22 

-0.18 
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Table 6.2. Continued 


(b) )J. = 0.20 


Component 

L x 

Ly 

L z 

A x 

Ay 

A z 

Input control vector {U} 

4C 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

4S 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

Hub and swashplate interface displacements {X 7 } 

H4C 

-0.000096 

-0.000353 




-0.000670 

H4S 

0.000632 

-0.000059 




0.002128 

S4C 

-0.000049 

-0.000075 




-0.000218 

S4S 

0.000214 

-0.000019 




0.000529 

H8C 

-0.000017 

0.000039 




0.000468 

H8S 

-0.000010 

0.000135 




0.000386 

S8C 

-0.000013 

0.000031 




0.000494 

S8S 

0.000006 

0.000053 




0.001194 


Component 

F X 

Fy 

F Z 

M x 

My 

M z 

Hub and swashplate interface forces and moments {Fj} 

H4C 

0.05 

-9.02 

6.78 


-9.08 


H4S 

10.70 

1.73 

-0.95 

24.78 

48.09 

79.31 

S4C 

0.23 

-0.16 

3.76 

-3.75 

1.18 

1.08 

S4S 

-0.26 

0.05 

2.63 

3.76 

-0.23 

0.74 

H8C 

-1.27 

1.25 

-8.23 

3.12 

-1.21 

-3.52 

H8S 

-1.13 

4.59 

-6.42 

-13.28 

-11.21 

-39.74 

S8C 

-0.01 

-0.06 

1.33 

0.05 

0.54 

0.39 

S8S 

-0.01 

-0.04 

0.63 

-0.07 


0.19 

Hub and swashplate forces and moments for fixed-hub condition {F 0 } 

H4C 

0.89 

-6.29 

6.66 

39.81 

-8.27 

-17.03 

H4S 


2.24 

-1.34 

24.43 

45.76 

79.99 

S4C 

0.24 

-0.16 

3.78 

-3.81 

1.21 

1.09 

S4S 

-0.26 

0.05 

2.75 


-0.22 

0.78 

H8C 

-0.80 

0.06 

-11.49 

3.27 

-0.81 

—4.95 

H8S 

-0.79 

-0.19 

-9.36 

-9.98 

-10.16 

—47.82 

S8C 

-0.01 

-0.06 

1.33 

-0.01 

0.46 

0.39 

S8S 

-0.01 

-0.03 



0.35 

0.20 

Impedance forces and moments [C] {X 7 } 

Ha 

-0.84 

-2.74 

0.12 

1.26 

-0.81 

-3.05 

H3S9B 

5.10 

-0.51 



2.33 

-0.69 

BBS 

0.00 

0.00 

-0.01 

0.06 

-0.02 

0.00 


0.00 

0.00 



-0.01 

-0.04 

H8C 

-0.47 

1.20 

3.26 

-0.15 

-0.41 

1.43 

H8S 

-0.34 

4,78 

2.94 

-3.30 

-1.05 

8.08 

S8C 

-0.01 

0.00 



0.08 

0.00 

S8S 

0.00 

0.00 



0.06 

-0.01 
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Table 6.2. Concluded 


(c) (i = 0.30 


Component 

Lx 

Ly 

L z 

A x 

Ay 

A z 

Input control vector { U } 

4C 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

4S 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

Hub and swashplate interface displacements (X/} 

rnrnsmm 

-0.000248 

-0.000408 


0.003031 




0.000589 

0.000212 


-0.001620 


0.002437 


-0.000108 

-0.000094 





S4S 

0.000200 

0.000033 


-0.000729 



H8C 

-0.000035 

0.000007 





H8S 

-0.000010 

0.000201 


-0.000956 



S8C 

-0.000011 

0.000017 





S8S 

0.000010 

0.000102 




0.001979 


Component 

F X 

Fy 

F Z 

M x 

My 

M z 

Hub and swashplate interface forces and moments {F/J 

H4C 

-1.85 

-9.50 

3.91 

52.23 

-24.10 

1 

H4S 

11.00 

9.27 

8.30 

-6.10 

37.29 

88.27 

S4C 

0.97 

0.30 

7.95 

-14.08 

-8.16 

2.16 

S4S 

-0.93 

0.02 

0.39 

13.41 

1.57 


H8C 

-0.49 

-0.16 

-4.59 

1.76 

-17.05 

11.94 

H8S 

-1.99 

7.23 

-8.11 

-0.69 

-10.63 

-53.83 

S8C 

-0.43 

0.03 

-1.33 

6.42 

1.29 

-0.29 

S8S 

0.01 

-0,30 

0.87 

-0.79 

3.92 


Hub and swashplate forces and moments for fixed-hub condition {F 0 } 

H4C 

0.49 

-6.19 

4.04 

55.23 

-24.67 

-15.62 

H4S 

6.33 

1.76 

7.73 

-3.88 

40.23 

104.27 

S4C 

0.98 

0.31 

8.02 

-14.16 

-8.28 

2.18 

S4S 

-0.93 

0.01 

0.53 

13.29 

1.58 

0.17 

H8C 

0.67 

-0.18 

-7.08 

1.14 

-16.33 

15.64 

H8S 

-1.67 

0.26 

-11.60 

3.25 

-9.54 

-66.27 


-0.45 

0.03 

-1.48 

6.79 

1.35 

-0.33 


0.00 

-0.31 

0.86 

-0.78 

4.10 

0.36 

Impedance forces and moments [C] {X;} 

H4C 

-2.34 

-3.31 

-0.14 

-3.01 

0.57 

^.43 

H4S 

4.68 

7.50 

0.57 

-2.22 

-2.95 


S4C 

-0.01 

-0.01 

-0.07 

0.08 

0.12 


S4S 

-0.01 

0.00 

-0.14 

0.11 

-0.01 


H8C 

-1.16 

0.02 

2.49 

0.61 

-0.72 

-3.70 

H8S 

-0.33 

6.97 

3.49 

-3.94 

-1.09 

12.44 

S8C 

0.02 

0.00 

0.14 

-0.37 

-0.07 


S8S 

0.00 

0.01 

0.01 

-0.01 

-0.19 
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Table 6.3. Coupled Rotor-Testbed Results for Actuated Motion 
(a) 0.01 in. 4P axial cosine actuation; jj. = 0.10 


Component 


Ly 

Lz 

A x 

Ay 

A z 

Input control vector 

{U} 

4C 

0.010000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

4S 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

Hub and swashplate interface displacements {X 7 } 

H4C 

0.010146 

-0.000448 

0.000095 

0.004010 

0.025044 

-0.004780 

H4S 

0.002971 

-0.000358 

0.000265 

0.002768 

0.012598 

0.001454 

S4C 

0.006981 

-0.000032 

0.000086 

0.001466 

0.017423 

-0.001232 

S4S 

0.001412 

-0.000090 

0.000261 

0.000777 

0.007917 

0.000356 

H8C 

0.000017 

0.000137 

0.000024 

-0.000819 

0.000401 

0.000242 

H8S 

0.000008 

-0.000531 

-0.000148 

0.000745 

-0.000271 

-0.004271 

S8C 

-0.000016 

0.000060 

0.000025 

-0.000155 

0.000043 

0.001358 

S8S 

0.000029 

-0.000420 

-0.000152 

0.000822 

-0.000043 

-0.003970 


Component 

Fx 

Fy 

Fz 

M x 

My 

M z 

Hub and swashplate interface forces and moments { F/} 

H4C 

83.11 

-14.72 

22.92 

62.23 

8.06 


H4S 


-3.26 

9.81 

78.27 

71.01 


S4C 


-0.09 


-4.80 

0.50 

0.01 

S4S 

-0.10 

0.03 

7.55 

1.94 

-0.72 

2.12 

H8C 

1.04 

4.24 

-2.35 

-13.54 

12.23 


H8S 

1.30 

-0.27 

11.83 

-6.11 

-17.48 


S8C 

-0.21 

0.04 

-0.35 

3.05 

0.00 

-0.11 

S8S 

0.05 



-0.08 

-2.43 

-1.10 

1.35 

-0.66 

Hub and swashplate forces and moments for fixed-hub condition {F 0 } 

H4C 

-0.22 

-11.11 

19.25 

58.31 

-3.71 

-122.66 

H4S 

7.22 

2.89 

7.52 

70.28 

64.34 

57.23 

S4C 

0.43 

-0.10 

-0.33 

-6.53 

0.43 

-0.07 

S4S 

-0.07 

-0.04 

7.32 

1.31 

0.33 

2.06 

H8C 

-0.34 

-0.31 

-2.58 

—4.98 

-7.73 

-20.14 

H8S 

0.90 

0.09 

6.06 

-2.23 

-2.84 

10.06 

S8C 

0.00 

-0.04 

0.23 

-0.12 

0.44 

0.07 

S8S 

0.03 

0.02 

-0.31 

-0.46 

-0.25 

-0.09 

Impedance forces and moments [C] {X/} 

H4C 

83.34 

-3.61 


3.92 

11.77 

-50.39 

H4S 

25.98 

-6.16 


7.99 

6.67 

1.61 

S4C 

-0.12 

0.02 

0.31 

1.74 

0.07 

0.08 

S4S 

-0.03 

0.08 

0.24 

0.63 

-1.05 

0.06 

H8C 1 

1.38 

4.55 

0.24 

-8.56 

19.96 

-35.63 

H8S 

0.40 

-0.36 

5.76 

-3.88 

-14.64 

-21.67 

S8C 

-0.21 

0.08 

-0.58 

3.16 

-0.44 

-0.19 

S8S 

0.02 

-0.10 

-2.12 


1.60 

-0.58 
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Table 6.3. Continued 

(b) 0.01 in. 4P axial sine actuation; p = 0.10 


Component 

L x 

L y 

L z 

A x 

Ay 

A z 

Input control vector 

{U} 

4C 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

4S 

0.010000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

Hub and swashplate interface displacements {X/} 


-0.002280 

-0.000559 

-0.000160 

0.003915 

-0.008200 

-0.004179 


0.011221 

-0.000041 

0.000184 

0.001161 

0.031117 

0.001325 


-0.001201 

-0.000151 

-0.000167 

0.001450 

-0.005849 

-0.001072 

■ 

0.007370 

0.000047 

0.000180 

0.000029 

0.020460 

0.000376 

H8C 

0.000007 

0.000095 

0.000030 

-0.000605 

0.000174 

0.000145 

H8S 

0.000007 

-0.000370 

-0.000109 

0.000477 

-0.000149 

-0.003004 


-0.000008 

0.000040 

0.000031 

-0.000104 

0.000040 

0.000936 


0.000019 

-0.000295 

-0.000112 

0.000576 

-0.000034 

-0.002776 


Component 

F x 

Fy 

F Z 

M x 

My 

M z 

Hub and swashplate interface forces and moments {F/} 

H4C 

-20.11 1 

-13.06 


64.38 

-10.89 

■ 

H4S 

99.70 

3.34 

9.58 


79.69 


S4C 

0.38 

-0.05 



-0.11 


S4S 

-0.07 

0.08 

7.69 

1.67 

-1.48 

2.16 

H8C 

0.42 

2.93 


-11.48 

4.37 

-39.57 

H8S 

1.13 

-0.10 


-5.97 


-8.98 

S8C 

-0.15 

0.08 


2.39 

-0.78 


S8S 

-0.01 

-0.04 

-1.86 

-0.03 

0.92 


Hub and swashplate forces and moments for fixed-hub condition {Fq} 

H4C 

-0.22 

- 11.11 

H9| 

58.31 

-3.71 

-122.66 

H4S 

7.22 

2.89 


70.28 

64.34 

57.23 

S4C 

0.43 

-0.10 

-0.33 

-6.53 



S4S 

-0.07 

-0.04 

7.32 

1.31 

0.33 


H8C 

-0.34 

-0.31 

-2.58 

—4.98 

-7.73 


H8S 

0.90 

0.09 

6.06 

-2.23 

-2.84 


S8C 

0.00 

-0.04 

0.23 




S8S 

0.03 

0.02 

-0.31 




Impedance forces and moments [C] {X/} 

H4C 

-19.89 

-1.95 

0.91 

6.08 

-7.18 

-28.94 

H4S 

92.48 

0.44 

2.06 

-5.18 

15.36 

-6.31 

S4C 


0.05 

-0.02 

0.94 

-0.55 

-0.01 

S4S 


0.12 

0.38 

0.36 

-1.80 

0.10 

H8C 


3.24 

-0.02 

-6.50 

12.09 

-19.43 

H8S 


-0.19 

2.74 

-3.74 

-7.88 

-19.04 

S8C 


0.12 

-0.42 

2.51 

-1.22 

-0.15 

S8S 


-0.06 

-1.55 

0.43 

1.16 

-0.43 
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Table 6.3. Continued 


(c) 0.01 in. 4P side sine actuation; p = 0.10 


Component 

L X 

Ly 

L Z 

A x 

Ay 

A z 

Input control vector (U) 

4C 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 


fS , 

0.000000 

0.010000 

0.000000 

0.000000 

0.000000 


Hub and swashplate interface displacements {X/} 

H4C 

-0.000297 

-0.001664 




-0.003178 

H4S 

0.000810 

0.007533 


-0.014586 


-0.004463 

S4C 

-0.000121 

-0.000739 

-0.000008 



-0.000581 

S4S 

0.000277 

0.005673 

0.000087 



-0.005397 

■ 

-0.000002 

0.000067 

0.000040 



0.000160 


0.000002 

-0.000127 




-0.001 1 14 

S8C 

0.000001 

0.000028 

0.000041 



0.000624 

S8S 

0.000005 

-0.000105 

-0.000062 



-0.001018 


Component 

Fx 

Fy 

Fz 

M x 

My 

M z 

Hub and swashplate interface forces and moments {F/} 

H4C 


-25.56 


67.41 

-5.37 


H4S 


65.13 


63.14 

65.69 


S4C 


-0.14 

-0.38 

-5.94 

1.13 

-0.08 

S4S 


-0.02 

7.85 

2.22 

0.03 


H8C 

-0.53 

2.12 

-3.42 

-7.25 

-2.56 


H8S 

0.96 

-0.29 

4.81 

-5.37 

-4.44 

-3.41 

S8C 

0.00 

0.01 

0.29 

0.00 

-0.19 

0.08 

S8S 

0.00 

0.05 

-1.01 

0.04 

-0.55 

-0.29 

Hub and swashplate forces and moments for fixed-hub condition {F 0 } 

H4C 

-0.22 

- 11.11 


58.31 

-3.71 


H4S 

7.22 

2.89 


70.28 

64.34 

■ ■ 

S4C 

0.43 

-0.10 

-0.33 

-6.53 

0.43 

-0.07 

S4S 

-0.07 

-0.04 

7.32 

1.31 

0.33 

2.06 

H8C 

-0.34 

-0.31 

-2.58 

-4.98 

-7.73 

-20.14 

H8S 

0.90 

0.09 

6.06 

-2.23 

-2.84 

10.06 

S8C 

0.00 

-0.04 

0.23 

-0.12 

0.44 

0.07 

S8S 

0.03 

0.02 

-0.31 

-0.46 

-0.25 

-0.09 

Impedance forces and moments [C] {X ; } 

H4C 

—4.35 

-14.45 




-3.84 

H4S 

6.05 

62.24 




-7.70 

S4C 

-0.05 

-0.04 


0.59 

0.69 

-0.01 

S4S 

-0.06 

0.03 

0.54 

0.91 

-0.30 

0.15 

H8C 

-0.19 

2.43 

-0.84 

-2.27 

5.17 

-2.99 

H8S 

0.07 

-0.39 

-1.25 

-3.14 

-1.60 

-13.48 

S8C 

0.00 

0.04 

0.05 

0.11 

-0.63 

0.01 

S8S 

-0.03 

0.03 

-0.70 

0.49 

-0.30 

-0.21 
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Table 6.3. Continued 

(d) 0.01 in. 4P normal sine actuation; |X = 0.10 


Component 

Lx 

Ly 

Lz 


Ay 

A z 

Input control vector {U} 

4C 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

4S 

0.000000 

0.000000 

0.010000 

0.000000 

0.000000 

0.000000 

Hub and swashplate interface displacements {X/} 


-0.002530 

-0.000574 

-0.001324 

0.004019 

-0.009618 

-0.004126 

H4S 

0.008906 

-0.000037 

0.010624 

0.001136 

0.035705 

0.001467 

S4C 

-0.001268 

-0.000153 

-0.001325 

0.001511 

-0.00703 1 

-0.001065 

S4S 

0.004307 

0.000049 

0.010590 

0.000021 

0.026734 

0.000408 

H8C 

0.000014 

0.000089 

0.000020 

-0.000590 

0.000335 

0.000058 

H8S 

0.000004 

-0.000486 

-0.000116 

0.000605 

-0.000247 

-0.003919 

S8C 

-0.000013 

0.000036 

0.000021 

-0.000089 

0.000034 

0.000929 

S8S 

0.000023 

-0.000389 

-0.000119 

0.000757 

-0.000033 

-0.003599 


Component 

F x 

Fy 

F Z 

M x 

My 

M z 

Hub and swashplate interface forces and moments { F/} 

H4C 

-22.07 

-13.76 

8.89 

63.41 

-15.92 

-149.19 

H4S 

80.54 

3.38 


64.39 

83.18 

56.28 

S4C 

0.35 

-0.01 

-0.59 

-5.08 



S4S 

0.02 

0.11 

6.99 

0.40 

-2.07 

1.96 

H8C 

0.85 

2.82 

-1.89 

-12.15 

10.31 

-43.74 

H8S 

1.12 

0.22 

9.63 

-7.86 

-15.89 

-13.26 

S8C 

-0.16 

0.14 

-0.29 

2.65 

-1.52 


S8S 

-0.06 

-0.04 

-2.25 



■kBHI 


Hub and swashplate forces and moments for fixed-hub condition {F 0 } 

H4C 


- 11.11 

19.25 

58.31 

-3.71 

-122.66 

H4S 

7.22 

2.89 

7.52 

70.28 

64.34 

57.23 

S4C 


-0.10 


-6.53 

0.43 


S4S 


-0.04 

7.32 

1.31 

0.33 


H8C 


-0.31 

-2.58 

—4.98 

-7.73 


H8S 


0.09 


-2.23 

-2.84 


S8C 


-0.04 


-0.12 

0.44 


S8S 


0.02 





Impedance forces and moments [C] {X/} 

H4C 

-21.84 

-2.66 


5.10 

-12.21 

-26.53 

H4S 

73.32 

0.48 

78.52 

-5.89 

18.85 

-0.95 

S4C 

-0.08 

0.09 


1.46 

-1.12 

-0.08 

S4S 

0.09 

0.15 



-2.39 

-0.10 

H8C 

1.19 

3.13 


-7.17 

18.03 

-23.60 

H8S 

0.23 

0.13 

3.56 

-5.63 

-13.04 

-23.32 

S8C 

-0.16 

0.17 


2.77 

-1.95 

-0.19 

S8S 

-0.09 

-0.05 

| -1.94 


1.28 

-0.54 
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Table 6.3. Concluded 
(e) 0.05° 4P pitch sine actuation; p = 0.10 


Component 

** 

Ly 

Lz 

A x 

Ay 

A z 

Input control vector {U} 

4C 


0.000000 

0.000000 

0.000000 


0.000000 

4S 


0.000000 

0.000000 

0.000000 


0.000000 

Hub and swashplate interface displacements {X/} 

H4C 

-0.004462 

-0.000414 

-0.000951 

0.002686 

-0.017315 

-0.006197 

H4S 

0.036147 

0.000062 

0.000951 

0.000475 

0.149276 

0.001169 

S4C 

-0.002167 

-0.000150 

-0.000958 

0.000812 

-0.012095 

-0.001479 

S4S 

0.016474 

0.000076 

0.000941 

-0.000235 

0.118428 

0.000350 

H8C 

0.000043 

0.000272 

0.000035 

-0.001572 

0.001102 

0.000468 

H8S 

0.000025 

-0.001179 

-0.000306 

0.001741 

-0.000468 

-0.009324 

S8C 

-0.000045 

0.000124 

0.000036 

-0.000303 

0.000065 

0.002660 

S8S 

0.000061 

-0.000924 

-0.000315 

0.001822 

-0.000082 

-0.008778 i 


Component 


Fy 

F z 

M x 

My 

M z 

Hub and swashplate interface forces and moments {F /} 

H4C 

-40.69 

-5.94 

21.87 

65.56 

-32.03 

-231.03 

H4S 

306.05 

4.78 

17.19 

51.84 

150.92 

46.72 

S4C 

0.08 

0.26 

-0.76 


-3.79 

0.23 

S4S 

0.11 

0.45 

9.52 


-7.55 

2.66 

H8C 

3.16 

7.48 

-2.83 

-28.32 

35.42 

-109.40 

H8S 

2.73 

-0.86 

25.11 

-10.89 

-33.12 

-20.56 

S8C 

-0.57 

0.40 

-1.53 

9.12 

-3.93 

-0.54 

S8S 

-0.16 

-0.22 

-5.69 

1.30 

4.34 

-1.56 

Hub and swashplate forces and moments for fixed-hub condition {Fq} 

H4C 

-0.22 

- 11.11 

19.25 

58.31 

-3.71 


H4S 

7.22 

2.89 

7.52 


64.34 

■ 

S4C 

0.43 

-0.10 


-6.53 


-0.07 

S4S 

-0.07 

-0.04 

7.32 

1.31 

I 

2.06 

H8C 

-0.34 

-0.31 

-2.58 

-4.98 

-7.73 

-20.14 

H8S 

0.90 

0.09 


-2.23 

-2.84 

10.06 

S8C 

0.00 

-0.04 


-0.12 

0.44 

0.07 

S8S 

0.03 

0.02 

-0.31 

-0.46 


-0.09 

Impedance forces and moments [C] {X;} 

H4C 

—40.47 

5.17 

2.62 


-28.32 

-108.37 

H4S 

298.83 

1.88 

9.66 


86.59 

-10.52 

S4C 

-0.35 

0.36 



-4.22 

-0.16 

S4S 

0.18 

0.50 

2.21 

-1.33 

-7.87 

0.60 

H8C 

3.51 

7.79 

-0.25 

-23.34 

43.14 

-89.26 

H8S 

1.84 

-0.95 

19.04 

-8.66 


-30.63 

S8C 

-0.57 

0.44 

-1.77 

9.24 

-437 

-0.61 

S8S 

-0.19 

-0.23 

-5.38 

1.76 

4.59 

-1.48 
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Table 6.4. Coupled Rotor-Testbed Results for Balance Loads With No Actuation and p = 0.10 



L x 

Ly 

L z 

A x 

Ay 

A z 

Input control vector 

{U} 

4C 

0.000000 

0.000000 

0.000000 



0.000000 

4S 

0.000000 

0.000000 

0.000000 



0.000000 

Fixed- system balance forces and moments 


»? 

Fr 

Fz 

M x 

My 

M z 

4C 

“2.51 

-10.75 

17.18 

431.29 

-76.44 


4 S 

13.60 

-1.94 

29.26 

78.72 


81.62 

8C 

-1.19 

-4.45 

3.10 

-26.03 

22.29 

99.79 

8S 

3.01 

4.46 

-7.26 

9.16 

-39.08 

-259.68 

Rotating-system balance forces and moments 

4C 

-1.76 

-16.13 

19.28 

129.29 

-9.53 

-132.95 

4S 

14.74 

0.68 

8.43 

68.44 

128.59 

52.77 

8C 

-0.26 

1.82 

-2.50 

-15.72 

-3.39 

-24.14 

8S 

0.76 

-0.37 

4.83 

-4.34 

-2.41 

-6.63 
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Table 6.5. Coupled Rotor-Testbed Results for Balance Loads With Actuated Motion 
(a) 0.01 in. 4P axial cosine actuation; p = 0.10 


Component 

L x 

Ly 

L z 


Ay 

A z 

Input control vector 

{U} 

4C 


0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

4S 


0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

Fixed-system balance forces and moments 



Fy 

F Z 

M x 

My 

M z 

4C 

281.67 


96.53 

403.60 

4573.81 


4S 



60.29 

182.04 

1616.29 


8C 

-2.73 

-10.12 

2.25 

-53.92 

46.87 

240.02 

8S 

7.49 

13.13 

-16.67 

38.19 

-95.22 

-727.70 

Rotating-system balance forces and moments 

4C 

89.98 

-15.33 

22.98 

125.14 

373.36 


4S 

35.14 

-3.36 

9.99 

92.55 

215.35 

58.87 

8C 

1.05 

3.95 

2.25 

-32.23 

16.67 

-55.76 

8S 

1.35 

-1.49 

-16.67 

-1.91 

-11.94 

-11.77 


(b) 0.01 in. 4P axial sine actuation; p = 0.10 


Component 

L X 

Ly 

Lz 

A x 

Ay 

A z 

Input control vector 

{U} 

4C 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

4S 

0.010000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

Fixed-system balance forces and moments 


F x 

Fy 

F z 

M x 

My 

M z 

4C 

-38.90 

-8.60 

-10.85 


-1080.83 

-225.16 

4S 

298.53 

0.41 

107.09 

6.98 

5282.36 

83.55 

8C 

-2.06 

-7,14 

3.06 

-38.51 

36.22 

165.62 

8S 

5.30 

9.19 

-12.38 

24.94 

-61.96 

-508.44 

Rotating-system balance forces and moments 

4C 

-21.64 

-13.47 


120.26 


-151.62 

4S 

107.24 

3.08 

9.70 

51.20 

516.49 

50.94 

8C 

0.42 

2.72 

-2.52 

-24.41 

6.16 

-39.57 

8S 

1.17 

-0.96 

8.49 

-3.41 

-5.88 

-9.10 
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Table 6.5. Continued 


(c) 0.01 in. 4P side sine actuation; ji = 0.10 


Component 

L x 

Ly 

L z 

A x 

Ay 

A z 

Input control vector { U } 

4C 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

4S 

0.000000 

0.010000 

0.000000 

0.000000 

0.000000 

0.000000 

Fixed-system balance forces and moments 


Fx 

Fy 

F z 

M x 

My 

M z 

4C 

-5.07 

-34.87 

15.21 

803.95 

-176.86 


4S 

12.46 

271.96 

29.62 

-2967.45 

564.00 


8C 

—1.55 

-4.71 

4.43 

-30.34 

26.98 

1 14.29 

8S 

2.65 

3.44 

-6.91 

8.63 

-34.14 

-187.94 

Rotating- system balance forces and moments 

4C 

- 4.77 

-28.03 


■ EZEEH 


-126.51 

4S 

13.76 

69.94 




49.46 

8C 

-0.53 

2.01 

-3.31 

-16.58 

—4.80 

-23.12 

■ 

0.97 

-0.53 

4.64 

-3.39 

-0.36 

-3.46 


(d) 0.01 in. 4P normal sine actuation; p. = 0.10 


Component 

Lx 

Ly 

L Z 

A x 

Ay 

A z 

Input control vector { U } 

4C 

0.000000 1 

0.000000 


0.000000 

0.000000 

0.000000 

4S 

0.000000 

0.000000 


0.000000 

0.000000 

0.000000 

Fixed-system balance forces and moments 


Fx 

Fy 

F Z 

M x 

My 

M z 

4C 

-36.06 

-9.15 

-76.10 

386.21 

-1454.16 

-223.49 

4S 

104.74 

0.38 

643.67 

4.13 

6517.03 

90.51 

8C 

-2.17 

-7.33 

1.88 

-35.06 

37.39 

160.87 

8S 

5.96 

11.59 

-13.09 

29.54 

-75.54 

-657.14 

Rotating-system balance forces and moments 

4C 

-23.74 

-14.12 

7.95 

122.30 

-1 12.69 


4S 

86.21 

3.90 

93.52 


435.56 


8C 


2.56 

-1.83 

-24.58 

13.95 

-43.74 

8S 

1.17 

-1.00 

9.30 

-5.99 

-11.10 

-13.41 
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Table 6.5. Concluded 
(e) 0.05° 4P pitch sine actuation; p = 0.10 


Component 

L x 

Ly 

L z 


Ay 

A z 

Input control vector {U} 

4C 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

4S 

0.000000 

0.000000 

0.000000 

0.000000 

0.050000 

0.000000 

Fixed-system balance forces and moments 


F X 

F r 

Fz 

Mx 

My 

M z 

4C 

-69.47 

-2.04 

-68.55 

199.58 

-2308.30 

-313.26 

4 S 

294.59 

2.46 

538.40 

-53.98 

16243.61 

78.74 

8C 

-5.86 

-19.37 

2.74 

-95.82 


464.69 

8S 

15.30 

29.99 

-34.56 

92.89 



Rotating-system balance forces and moments 

4C 


-5.91 

21.20 


-210.40 


4S 

|f|K£ 

4.79 

17.85 

1 

1492.91 

46.72 I 

8C 

3.18 

6.99 

-2.73 

-61.48 

48.93 


8S 

2.86 

-3.61 

24.25 

-0.45 

-21.39 

-20.92 
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Table 6.6. Coupled Rotor-Testbed Results for Actuated Motion Corrected for 
Proper Motion at Hub Interface Joint 

[0.01 in. 4P axial sine actuation; p = 0.10] 




Ly 

L z 

A x 

Ay 

A z 

Input control vector {U) 

4C 

0.001434 

-0.001213 

-0.000122 

-0.003473 

0.002847 

0.002258 

4S 

0.025596 

-0.000581 

0.000042 

-0.001358 

-0.025542 

-0.001503 

Hub and swashplate interface displacements (X/) 

H4C 

-0.000014 

0.000008 





H4S 

0.009910 

0.000006 





S4C 

-0.000008 

-0.000123 





S4S 

0.010421 

-0.000076 





H8C 

0.000002 

0.000077 





H8S 

0.000004 

-0.000262 





S8C 

-0.000005 

0.000031 





S8S 

0.000013 

-0.000211 






Component 

Fx 

Fy 

Fz 

M x 

My 

M z 

Hub and swashplate interface forces and moments {F/} 

H4C 

-1.35 

-9.03 


63.41 

-3.83 

-142.85 

H4S 

89.32 

3.30 

7.73 

67.01 

57.86 

57.54 

S4C 

0.46 

-0.13 

-0.51 

-6.99 

0.81 

HU 

S4S 

-0.16 

0.03 

7.64 

2.82 

-0.53 

2.14 

H8C 

0.16 

2.39 

-2.61 

-9.68 

0.76 

-30.76 

H8S 

0.94 

-0.03 

6.70 

-4.71 

-8.58 

-7.22 

S8C 

-0.1 1 

0.04 

-0.01 

1.64 

-0.33 


S8S 

0.01 

-0.03 

-1.35 

-0.26 

0.52 

-0.38 1 

Hub and swashplate forces and moments for fixed-hub condition {F 0 J 

H4C 


-11.11 

19.25 

58.31 

-3.71 

-122.66 

H4S 

7.22 

2.89 

7.52 

70.28 

64.34 

57.23 

S4C 


-0.10 

-0.33 

-6.53 

0.43 


S4S 


-0.04 

7.32 

1.31 

0.33 


H8C 


-0.31 

-2.58 

-4.98 

-7.73 


H8S 


0.09 


-2.23 

-2.84 


S8C 


-0.04 



0.44 


S8S 


0.02 


-0.46 

-0.25 


Impedance forces and moments [C]{X/J 

H4C 



1.68 

5.10 

-0.12 


H4S 


'W& 

0.20 

-3.27 

-6.48 


S4C 

0.02 

-0.03 

-0.18 

-0.45 

0.38 


S4S 

-0.09 

0.07 


1.51 

-0.86 

1 

H8C 

0.50 

2.70 


—4.71 

8.49 

■ 

H8S 

0.05 

-0.12 


-2.48 

-5.74 

-17.28 

S8C 

-0.11 

0.08 


1.75 

-0.77 


S8S 

-0.03 

-0.04 

-1.04 


0.77 

| -0.29 


63 



























Actuator 5 length, in. Actuator 3 length, in. Actuator 1 length, in. 



Time, s Time, s 

Figure 6.1. Actuator lengths for 0.01 in. 4P sine actuation of inflexible ARES-IL 
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Actuator 6 displacement from nominal, in. Actuator 4 displacement from nominal, in. Actuator 2 displacement from nominal, in. 









Hub cosine 
Hub sine 

Swashplate cosine 
Swashplate sine 
Total mag 



.2 

Advance ratio, p 
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.4 





.2 .3 

Advance ratio, p 


.4 




Figure 6.3. Fixed-system loads for uncoupled BRR model. 


:hing moment force, lb 


Advance ratio 


100 r 


50 r 



4P sine pitching moment force, lb 4P sine rolling moment 


balance 



with hub interface joint loads. 



7. Conclusions 

A method for coupling an aeroelastically scaled model rotor system to the ARES-II (second generation ver- 
sion of the Aeroelastic Rotor Experimental System) rotor research testbed has been developed. The method 
requires the development of individual analytical models of the ARES-II and the rotor system to be coupled. For 
this study the Basic Research Rotor (BRR), a generic research rotor, is used. The development of individual mod- 
els of the testbed and the rotor system is used to an advantage by selecting analyses that are best suited to model 
each system. For the ARES-II model, the Dynamic Analysis and Design System (DADS) multibody dynamics 
analysis is used. For the rotor system, the second generation version of the Comprehensive Analytical Model of 
Rotorcraft Aerodynamics and Dynamics (CAMRAD-II) is used. Other analyses are suitable for use in place of 
DADS and CAMRAD-II. 

Impedance matching is used to couple the BRR to the ARES-II. Impedance matching requires that the struc- 
tures to be coupled be represented with a linear model. This representation is of no consequence for the ARES-II, 
a linear system. However, since rotor systems are inherently nonlinear systems, care must be exercised to ensure 
the proper application of the impedance-matching approach. For the BRR, fixed-hub forward flight trim condi- 
tions for 1 g steady flight are used with a correction applied to account for the loads due to rotor hub motion. 

The process used to perform coupled rotor-testbed analysis is as follows: 

1. Develop a mobility matrix [A] for the ARES-II by using the DADS model. This matrix is generated by applying 
4P and 8P dynamic loads at the interface joints and characterizing the 4P and 8P motion at the interface joints. 

2. Develop a motion matrix [B] for the ARES-II by using the DADS model. This matrix is generated by prescribing 
motions of the Stewart platform and characterizing the motion at the interface joints. 

3. Develop the fixed-hub load conditions (no motion) at the interface joints for the BRR by using the CAMRAD-II 
model. Then the impedance matrix [C] is generated by prescribing hub motions and characterizing the loads at the 
interface joints. 

4. The resulting mobility, motion, and impedance matrices are coupled through simple matrix equations. 

This study has shown that the impedance-matching approach is a viable method for the analysis of model rotor 
systems coupled to the ARES-II aeroelastic rotor research testbed. The impedance-matching method has been 
used to an advantage by generating independent models of the ARES-II and the BRR system. Thus, to couple 
alternate rotor systems to the ARES-II requires only the generation of the impedance model for the new rotor sys- 
tem. Additionally, the experimental verification of the ARES-II model will be eased because it may be exercised 
in a stand-alone mode for comparison with shake and actuation tests. This study has extended the impedance- 
matching techniques evident in the literature by including the effects of actuated hub motion and by using two 
interface joints between the rotor system and the testbed. 

Based on the experience developed in rotor-body coupling techniques during this study and the examination of 
the results, the following conclusions were found: 

1. Although impedance matching is a linear theory approach, the effect of rotor system nonlinearity must be consid- 
ered throughout the analysis. The rotor impedance matrices for this study have been developed by using both 
cosine and sine harmonic hub motions with checks to ensure that the rotor system response remained in the linear 
range. 

2. The effect of both hub and control system (swashplate) loads is necessary to properly model the coupled rotor- 
testbed system. Control system loads have been shown to be of particular importance for the calculation of normal 
force and pitching and rolling moments. Neglecting the swashplate load path has been shown to result in fixed- 
system load errors of up to 30 percent. 

3. Significant ARES-II body motion coupling and phase lag are evident in the actuated motion results for both the 
stand-alone ARES-II and for the coupled rotor-testbed system. Stewart platform control vectors required to obtain 
pure hub motion uncontaminated by body coupling or phase lag may be developed for the stand-alone ARES-II 
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with simple matrix manipulations. However, no closed-form method of calculating the required control vector is 
possible for the coupled rotor- testbed system. Therefore, the analytical results indicate that a closed-loop control 
system is necessary to generate desired hub motions. 

4. Based on the analytical results, the ARES-II rotating balance is expected to be a more reliable loads measurement 
device than the fixed-system balance. The results suggest that the rotating balance is essentially free of loads con- 
tamination due to system dynamic response. The fixed- system balance results indicate significant errors for all 
load cases examined. 


NASA Langley Research Center 
Hampton, VA 23681-2199 
December 16, 1997 
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Appendix A 


Rotating to Nonrotating Coordinate Transformations 

In the study of helicopter rotor system dynamics, the transformation of displacements and loads 
from the rotating frame to the nonrotating frame is of prime importance. This appendix describes the 
mathematical formulation of such transformations. Special attention is given to four-bladed articulated 
rotor systems to emphasize the equations necessary for analysis of the Basic Research Rotor (BRR). 

Figure A1 illustrates a nonrotating (or fixed-system) coordinate frame represented by the X,T,Z axes 
and a rotating-system coordinate frame represented by the x,y,z axes. Although the nonrotating coordi- 
nate frame is often also called the fixed system, it should not be confused with the ground-fixed inertial 
system. The nonrotating coordinate frame is a body-fixed system placed within the helicopter or testbed 
fuselage. The orientation in use throughout this work is +X-axis aft, +T-axis right side, and -i-Z-axis up 
along the rotor shaft axis with origin placed at the center of rotation. The rotating coordinate frame 
(x,y,z) is a hub-fixed coordinate system that revolves about the Z-axis with constant angular velocity Q 
whose origin is placed at a constant radius e from the Z-axis. For the purposes of this development, the 
origin of the rotating frame is constrained to lie in the X-Y plane at all times and the distance e is defined 
as the rotor flap-lag hinge offset. For the BRR, the flap and lag hinges are coincident; this allows for a 
single value of e . For rotor hubs with noncoincident flap and lag hinges, two values of e are necessary in 
the transformation formulation. The x- and y-axes are constrained to the X-Y plane and oriented such 
that the +x-axis lies toward the trailing edge and the y-axis lies along the radius of the reference blade. 
The +z-axis is oriented up and is parallel to the Z-axis. The azimuthal position of the rotating frame rel- 
ative to the nonrotating frame is given by \| /. The azimuth is indexed such that \\f = 0° when the refer- 
ence blade is over the X-axis. This position corresponds to the one when the +y-axis is aligned with the 
+X-axis. The azimuth angle increases with counterclockwise rotation of the rotating frame. 

By inspection of figure Al, the displacement, force, and moment relationships between the refer- 
ence blade root and the nonrotating frame are evident and are 
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(A3) 


As a simple example, one may define the vertical shear force in the nonrotating frame (F z ) due to verti- 
cal blade root shear (F z ) by using equation (A2). The result is 


F 7 = F 
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z 



Figure Al. Rotating and nonrotating frame coordinate systems. 


However, note that this result only includes the effect of a single blade. For multibladed rotor systems, deter- 
mining the forces and moments in the nonrotating frame by including the effect of all blades is necessary. There- 
fore, the equation for the total nonrotating vertical shear is 





(A4) 


where b is the blade index and N is the number of blades. Assume now a Fourier series to describe the periodic 
force for each blade: 



X (\. cos n Vb + F z ns sin "V*) 
n= 0 


(A5) 


where the index n is harmonic number, the subscripts c and s indicate the cosine and sine components of the force, 
and b is the azimuthal position of the bth blade as described by 


Vb = V + -ff (A6) 

where y is the azimuthal position of the reference blade. Substituting equation (A6) into equation (A5) and the 
result into (A4) yields the equation for the total nonrotating frame vertical shear as 

N 

F Z m = X I ( F z nc cos n Vb + F z ns sin »V*) ( A7 ) 

b = 1 n=0 

Further development, however, will show that equation (A7) may be further reduced. To aid in this development 
several trigonometric identities are used as follows: 

sin(ot + (3) = sin a cos p + cos a sin p (A8) 
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cos(a - P) = cos a cos p + sin a sin p 
cos(a + P) = cos a cos P - sin a sin P 
sin(a - P) = sin a cos P - cos a sin P 


(A9) 
(A 10) 
(All) 


Reduction of equation (A7) is achieved by developing equations for one harmonic at a time, beginning with the 
0th harmonic (mean) and working through to the 5th harmonic. Because the BRR is a four-bladed rotor system, 
N = 4 is assumed throughout the remainder of the development. 


For n = 0, N = 4: 


g„ =0 = I^ 0c cos (°-^)^ Zoc Sin 

b— 1 


(0 • v»>] 


- x KJ 


b= 1 
= 4F z 

z 0c 


(A 12) 


Equation (A 12) shows the expected result: the nonrotating frame vertical shear is four times the vertical blade root 
shear. 


For n = 1 : 


'x.L, = 


b= 1 
4 

• x 

b= 1 


F 7 cos 

Z lc 


( bn' 

„ . f bn\ 

l v + T, 

) + ^„ sm l ¥ + Tj 


(A13) 


Equations (A10) and (A1 1) allow for reduction of the cosine and sine terms depending upon blade number. 
For blade 1 : 


cos 


¥ + = -sin ¥ 

^ = cos \\f 


For blade 2: 


sin \|/ + 


COS(\|/ + 7t) = -COS\|t 

sin(y + tt) = -siny 


For blade 3: 


cos 



= sin y 


sin 



-cos y 
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For blade 4: 

cos(\|/ + 2ji) = cosy 
sin(\|/ + 2;t) = sin \j/ 

Substituting these expressions into equation (A 13) yields 

F 7 I =F (-sin V)/ - cos \|/ + sin \|/ + cos \|/) 

tot ln= 1 Z]c 

+ F (cos \|/- sin \|/- cos \|f + siny) 

Z \s 

= 0 

Each harmonic may be evaluated similarly. 

For n = 2: 



= F (-cos 2\| f + cos 2 \j / - cos 2 w + cos 2 \i/) 

z 2c 


+ F (-sin 2vi/ + sin 2\i/ - sin 2 w + sin 2 \if) 

z 2s 

= 0 

For n = 3: 



= F (sin 3 w - cos 3\|/ - sin 3\i/ + cos 3xi/) 

z 3c 


+ F (-cos 3x | f - sin 3 \\f + cos 3x| / + sin 3 \\f) 

z 3 s 

= o 

For n = 4: 



= F ( cos 4\|/ + cos 4\|f + cos 4\(/ + cos 4\y) 

z 4c 


+ F ( sin 4x|/ + sin 4x | f + sin 4xj/ + sin 4xj/) 

Mi 


Therefore, 


= 41 F cos 4\)/ + F sin 4\|/ 



(A14) 


(A15) 
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For n = 5: 



= F (-sin 5\|/ - cos 5\)/ + sin 5\|/ + cos 5\\i) 

z 5 


+ F (cos 5\j/ - sin 5\|/ - cos 5 \\f + sin 5\|/) 

z 5s 

= 0 


Continued analysis would confirm that rotating frame, vertical blade root shear is converted to vertical nonrotating 
frame shear only if the harmonic number is an integer multiple of the number of rotor system blades. For a four- 
bladed rotor system, this may be represented mathematically by 



4 F 



( n = 0,4, 8, ...) (A 16) 


An analysis similar to the preceding development of the vertical shear loads may be extended to the transfer of 
chordwise blade root shears to yawing moments in the nonrotating frame; that is. 


M 7 = -4 eF Y 

^tot x nc 

nc 

M 7 = -AeF r 


(11 = 0,4,8,...) (A17) 


However, the manner in which axial (fore-and-aft) and side loads are generated in the nonrotating frame is not as 
straightforward. Therefore, the axial load equations are developed to ensure the clarity of the proper procedures. 

Equation (A2) shows that the axial nonrotating frame load F x is generated by a combination of chordwise 
blade root shears F x , and blade root axial loads F y . The total axial load may be represented by 

N 

F Xujt = 2 Sin Vb + F y COS V*) ( A18 > 

b= 1 

where 

oo 

F x = X l F x nc C0S n Vb + F x ns sin "V*) 
n = 0 

oo 

F y= 1 { F y nc cos nx *b + F y ns sin "V*) 

«= 0 


75 



and was defined in equation (A6). Expanding equation (A18) yields 


N 


F x xol -XX 


b — 1 n = 0 


(f^ cos + F x sin n\|/ fc ^sin \f b 


+ [ F y nc cos n\\i h + F y ^ sin m|/ fc jcosy fc 


(A 19) 


Again assuming a four-bladed rotor system and using the trigonometric identities in equations (A8) through (All) 
give the following results: 

For n = 0, N = 4: 


For n = 1 : 


tot ln=0 


- 1 F 


yoc 


b= 1 

= AF y, 

y^c 


(A20) 


N 


L.i - X 


b— 1 


F „ cos 

x \c 


bn} _ . ( bn\ . f bn 

v + Tj +F ',. s,n l v+ Tj S1 "( ¥ + T 


F y„ cos 


bn\ „ . ( bn'] ( bn\ 

¥+ tJ +f *, s " , ( ¥+ tJ co u ¥+ tJ 


= F (-sin \\f cos \\f - sin \|/ cos \|/ - sin \)/ cos \|f + sin \j/ cos \j /) 


l lc 


2 2 2 2 
+ F (cos \|f + sin \|/ + cos \|f + sin \|f) 




2 2 2 2 
+ F (sin \(/ + cos \\f + sin Y|/ + cos \(f) 

y\ c 

+ F (-sin iff cos iff + sin \|/ cos \|/ - sin \|/ cos \|f + sin \\f cos \f/) 

y\s 

= 2F X + IF 

x u y\ s 


(A21) 
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For n- 2: 


tot l «=2 


N 

- 1 

b= 1 


cos 2 ( v + f) + sin 2 ( v + y)l sin ( v + t ) 


cos :(v + f )+ F yu sin 2( ¥ + f )lcos (y + y 


= F (-cos 2xi/ cos xj/ - cos 2\|/ sin \)/ + cos 2\|/ cos \j/ + cos 2xj / sin X|/) 

*2c 

+ F (-sin 2x|/ cos V|/ - sin 2x|/ sin \|/ + sin 2 X|/ cos \|/ + sin 2\\f sinx|/) 

x 2s 

+ F (cos 2\i/ sin \j/ - cos 2x| / cos Xj/ - cos 2x|/ sin X|/ + cos 2xj/ cos X|/) 

Vie 

+ F (-sin 2 x(/ sin X|/ - sin 2x|/ cos X|/ - sin 2x| / sin X|/ + sin 2x}/ cos xj/) 
= 0 


For /i = 3: 



= F (sin 3xi/ cos X|/ + cos 3x|/ sin X| / + sin 3x|/ cos X|/ + cos 3x|/ sin X|/) 

x 3c 

+ F (-cos 3xi/ cos X)/ + sin 3x|/ sin X| / - cos 3x|/ cos X|/ + sin 3 \\i sin X|/) 

x 3s 

+ F (-sin 3x|/ sin xy + cos 3x|/ cos X}/ - sin 3x|/ sin X|/ + cos 3x|/ cos X|/) 

>3c 

+ F (cos 3x|/ sin xj/ + sin 3x(/ cosx)/ + cos 3x|/ sin X(/ + sin3xj/ cos X|/) 

^3j 

= 2 F sin 4xi/ - 2 F cos 4xi/ + 2F„ cos 4xi/ + 2 F„ sin 4x|/ 

x 3c *3* >3 c >3j 

= 2 (F r +F )sin 4x|/ + 2(-F +F )cos4i|/ 

*3c a 3j - y 3c 


(A22) 
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For n~ 4: 


F *c\ n =4 = 2 {K, C0S ( V + f) + F * 4 , Si " + y)] ( V + f) 

/?=1 ^ 

+ K, cos (v + y)+ f, 4 , *ta (v + f)] cos (y + f )} 

= F (cos 4\(/ cos \|/ - cos 4y sin y - cos 4\| / cos y + cos 4\j / sin \|/) 

*4c 

+ F^ ( sin 4\| / cos \j f - sin 4y sin \| / - sin 4\|/ cos \)/ + sin 4\|/ sin y ) 

'Mj 

+ F (-cos 4\|f sin \(/ - cos 4 \)/ cos y + cos 4\|/ sin \|/ + cos 4\(/ cos \j /) 

■Me 

+ F (-sin 4\|/ sin y + sin4y cosy + sin4\)/ siny + sin 4 V|/ cosy) 

= 0 

For n = 5: 


tot ln=5 


N 

- 1 

fr=l 


F^ cos 


bn} r . f bn\ . ( bn\ 

v+ tJ + f -< s , s,n l y tJ s1 H v+ tJ 


bn) „ . ( ( bn\ 

,| ' + yJ + ^ 5 , sln l v + Tj ““(v + yj 


F y sc cos 


= F (-sin 5\| / cos \| t + cos 5 y sin y - sin 5 y cos y + cos 5\j / sin y) 

x 5c 

+ F (cos 5y cos \(/ + sin 5 y sin \|/ + cos 5y cos y + sin 5 y sin y) 
+ F (sin 5y sin y + cos 5y cos y + sin 5y sin y + cos 5y cos y) 

Me 

+ F (-cos 5y sin y + sin 5y cos y - cos 5y siny + sin 5y cos y) 

Ms 

= 2 F r sin(-4\|/) + 2 F r cos 4w + 2 F„ cos 4\i/ + 2 F sin 4\i t 

x 5c x 5 s ysc T yss 

= 2 (-F + F \in 4 \\r + 2(V + F v Icos 4\|/ 

V x 5c y^sj v 5 j Vsc) 


(A23) 


Combining equations (A22) and (A23) yields the following cosine and sine components of the 4-per-rev axial load 
in the nonrotating frame: 




2(F X +F 

x 3c yis 



(A24) 


Finally, the following equations (eqs. (A25) through (A36)) provide a set of general results for each component of 
force and moment in the nonrotating frame due to blade root shear in the rotating frame. Note that these equations 
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are valid only for four-bladed rotor systems and that the subscript denoting total has been dropped from the nonro- 
tating frame force and moment terms. 


= 21 -F 


+ F V +F„ 

(n-1 )s -'(n-1 )c 
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(n = 0, 4, 8,...) 
(n — 0, 4, 8,...) 
(n = 0, 4, 8,...) 

(n = 0, 4,8,...) 

(n = 0, 4, 8, ...) 
(n = 0, 4, 8,...) 

(n = 0, 4, 8,...) 
(n = 0, 4, 8, ...) 
(n = 0, 4, 8, ...) 

(n = 0, 4, 8,...) 

(n = 0, 4, 8,...) 
(» = 0, 4, 8,...) 


(A25) 

(A26) 

(A27) 

(A28) 

(A29) 

(A30) 

(A31) 

(A32) 

(A33) 

(A34) 

(A35) 

(A36) 
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Appendix B 


ARES-II Model Data 

This appendix provides component data necessary for modeling the ARES-II with the Dynamic Analysis and 
Design System (DADS) software. Although all data are not specifically provided, an experienced DADS user can 
generate a working model similar to that developed for this study. This appendix is organized into tables of rigid- 
body, joint, elastic element, and miscellaneous element data. All locations are referenced to an inertial coordinate 
system oriented in the manner used throughout the report: +X, aft; +Y, right; +Z, up. The bottom center of the stand 
baseplate is placed at the inertial coordinate system origin. 

Bl. Rigid Bodies 

Table B 1 lists the rigid-body components necessary to model the ARES-II. Included in the table are the body 
center of gravity with respect to the inertial coordinate system, the mass of the body, and the inertia of the body 
with respect to the center of gravity. 


Table B 1 . ARES-II Model Rigid-Body Components 


Body description 

Mass, 

Center of Gravity, in. 

Mass moment of inertia about body c.g., lb-in-s 2 

lb-s 2 /in. 

X 

Y 

z 


*yy 

l zz 

*xy 

I* 

bz 

Stand base 


0 

0 

16.304 

704.22 

706.39 

265.79 

0 

0 

0 

Stand yoke 


-0.510 

0 

68.730 

3.3872 

0.4592 

3.1791 

0 

0 

0 

Longeron 

0.1571 

-4.529 

0 

71.691 

6.4772 

5.7451 

11.5442 

0 

-1.2215 

0 

Fixed-system balance 
bottom 

0.0871 

1.400 

0 

70.796 

0.6633 

0.6979 

1.1559 

0 

0 

0 

Fixed-system balance 
top 

0.1350 

1.400 

0 


1.3665 

1.3670 

2.1458 

0 

0 

0 

Main shaft housing 

0.1266 

1.010 

0 



1.0799 

0.4970 

0 

0 

0 

Main shaft and drive 
sheave 

0.0574 

1.400 

0 


1.4559 

1.4559 

0.5157 

0 

0 

0 

Intermediate shaft 
housing 

0.0179 

-6.574 

0 


0.0957 

0.1138 

0.1037 

0 

-0.0264 

0 

Intermediate shaft and 
drive sheaves 

0.0279 

-5.125 

0 

80.884 


0.3238 

0.1767 

0 

0 

0 

Drive motor 

0.1703 

-13.653 

0 

77.539 

1.8022 

1.8022 

0.8264 

0 

-0.0166 

0 

Drive motor sheave 

0.0018 

-13.795 

0 

83.337 

0.0008 

0.0008 

0.0010 

0 

0 

0 

Rotating balance 
bottom 

0.0020 

1.400 

0 

91.304 


0.0016 

0.0028 

! 0 

0 

0 

Rotating balance top 

0.0028 

1.400 

0 

92.577 


0.0019 

0.0034 

0 

0 

0 

Mast 

0.0064 

1.400 

0 

94.615 

0.0164 

0.0164 

0.0038 

0 

0 

0 


B2. Joints 

Five joint types are used to model the ARES-II: translational, revolute, spherical, universal, and bracket. The 
translational joint eliminates all degrees of freedom between two bodies except one linear dof. The revolute joint 
eliminates all dofs between two bodies except one angular dof. The spherical joint eliminates linear dofs but 
allows all angular dofs, and the universal joint eliminates all linear and one angular dof. The bracket joint elimi- 
nates all degrees of freedom between two bodies effectively making them one larger rigid body. Table B2 lists all 
joints, their location with respect to the inertial coordinate system, and the bodies that the joints connect. 
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Table B2. ARES-II Model Joints 


Joint type 

Joint location, in. 

Body 1 

Body 2 

X 

Y 

z 

Revolute 

-0.600 

0 

74.040 

Stand yoke 

Longeron 

Universal 

-7.675 

-6.580 

74.340 

Longeron 

Actuator 1 body 

Universal 

-7.675 

6.580 

74.340 

Longeron 

Actuator 2 body 

Universal 

2.402 

-7.100 

74.345 

Longeron 

Actuator 3 body 

Universal 

2.402 

7.100 

74.345 

Longeron 

Actuator 4 body 

Universal 

8.513 

-0.625 

75.790 

Longeron 

Actuator 5 body 

Universal 

8.513 

0.625 

75.790 

Longeron 

Actuator 6 body 

Spherical 

-5.806 

0 

70.353 

Actuator 1 Rod 

Fixed-system balance bottom 

Spherical 

-5.421 

0 

70.560 

Actuator 2 Rod 

Fixed-system balance bottom 

Spherical 

8.748 

-5.724 

69.961 

Actuator 3 Rod 

Fixed-system balance bottom 

Spherical 

8.748 

5.724 

69.961 

Actuator 4 Rod 

Fixed-system balance bottom 

Spherical 

8.513 

-5.450 

69.728 

Actuator 5 Rod 

Fixed-system balance bottom 

Spherical 

8.513 

5.450 

69.728 

Actuator 6 Rod 

Fixed-system balance bottom 

Bracket 

1.400 

0 

75.590 

Fixed-system balance top 

Main shaft housing 

Bracket 

-2.725 

0 

80.140 

Main shaft housing 

Intermediate shaft housing 

Revolute 

1.400 

0 

77.710 

Main shaft housing 

Main shaft and drive sheave 

Bracket 

1.400 

0 

90.915 

Main shaft and drive sheave 

Rotating balance bottom 

Revolute 

-5.125 

0 

79.941 

Intermediate shaft housing 

Intermediate shaft and drive sheaves 

Bracket 

-9.125 

0 

78.710 

Intermediate shaft housing 

Drive motor 

Revolute 

-13.795 

0 

82.471 

Drive motor 

Drive motor sheave 

Bracket 

1.400 

0 

92.915 

Rotating balance top 

Mast 


B3. Elastic Elements 

Three types of elastic elements are used to model the ARES-II: beam, bushing, and spring-damper. The bush- 
ing is a six-degree-of-freedom elastic element with damping. The nature of the beam and spring-damper are 
self-explanatory. Table B3 provides a list of all elastic elements, their connection locations in inertial coordinates, 
and the bodies that they connect. 

Table B3. ARES-II Model Elastic Elements 


Elastic 
element type 

Connection location for body 1 , 
in. 

Connection location for body 2, 
in. 

Body 1 

Body 2 

X 

Y 

Z 

X 

Y 

Z 



Spring damper 1 

-17.854 

-17.854 

0 

-17.854 

-17.854 

0 

Pad 1 

Stand base 

Spring damper 1 

-17.854 

-17.854 

0 

-17.854 

17.854 

0 

Pad 2 

Stand base 

Spring damper 1 

17.854 

17.854 

0 

17.854 

17.854 

0 

Pad 3 

Stand base 

Spring damper 1 

17.854 

-17.854 

0 

17.854 

-17.854 

0 

Pad 4 

Stand base 

Beam 2 

0 

0 

30.0 

0 

0 

66.840 

Stand base 

Stand yoke 

Bushing 3 

1.400 

0 

72.590 

1.400 

0 

72.590 

Fixed-system balance 
bottom 

Fixed-system balance 
top 

Bushing 4 

1.400 

0 

91.915 

1.400 

0 

91.915 

Rotating balance 
bottom 

Rotating balance top 


Sprint-damper properties: K z = 6000.0 lb/in.; C z = 20.0 lb- s/in. 

2 Beam properties: EI^ = 1.163 x 10 9 lb-in 2 ; EI^, = 1.163 x 10 9 lb-in 2 ; GJ = 8.95 x 10 8 lb-in 2 
3 Bushing properties: 

K x = 862,000.0 lb/in.; K y = 512,800.0 lb/in.; K z = 93,400.0 lb/in. 

K Qjc = 18.75 x 10 6 in-lb/rad; K^= 17.44 x 10 6 in-lb/rad; tf 0 , = 26.30 x 10 6 in-ib/rad 
C x - 8.6 lb-s/in.; C v = 5.1 lb-s/in.; C z = 0.9 lb-s/in. 

Cq x = 1875.0 in-lb- s/rad; C 0y = 1744.00 in-lb-s/rad; C 0 ^ 2630.0 in-lb-s/rad 
4 Bushing properties: 

K x = 709,000.0 lb/in.; K y = 709,000.0 lb/in.; K z = 2.78 x 10 6 Ib/in. 

Kq x =2.16 x 10 6 in-lb/rad; K« = 2.76 x 10 6 in-lb/rad; K e , = 2.82 x 10 6 in-lb/rad 
C x - 7. 1 lb-s/in.; C y = 7. 1 lb-s/in.; C z = 27.8 lb-s/in. 

C Qx = 276.0 in-lb-s/rad; Cq = 276.0 in-lb-s/rad; C Qr = 282.0 in-lb-s/rad 
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B.4. Miscellaneous Elements 


The miscellaneous elements consist of constraint and point-of-interest elements. The constraint elements are 
used to control the length of the ARES-II actuators, hold the longeron in a fixed position relative to the stand, con- 
strain the stand degrees of freedom that are not restrained by the mounting pads (i.e., lateral displacement and 
yaw), and develop the relationships for the drive train speed reduction ratio. The point-of-interest elements allow 
for the output of displacement information at a point on a rigid body. Table B4 provides the location of the point- 
of-interest elements used for the hub and swashplate interface joints. 


Table B4. ARES-II Model Point-of-interest Elements 


Body 

Location, in. 

Purpose 

X 

Y 

z 

Main-shaft housing 
Mast 

1.400 

1.400 

0 

0 

87.908 

96.128 

Swashplate interface joint 
Hub interface joint 
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Appendix C 

Basic Research Rotor Model Data 

The tables in this appendix provide the input necessary to develop the rotorcraft shell and core models of the 
Basic Research Rotor. Note that all CAMRAD-II coordinate frames are not necessarily the same as those used in 
this report. Table Cl provides the basic input data to the rotorcraft shell. Table C2 provides run-time modifications 
to the shell model and the core model inputs. 

During the research program it was noted that results obtained when modeling the R-12 test medium were 
inconsistent with expected trends. This inconsistency was eliminated if the model and environment parameters 
were scaled for air. The most obvious differences were for the high-speed flight regime (p > 0.30); however, some 
differences were noted at the lower speeds. The reason for the inconsistency is unknown as CAMRAD-II executes 
internally in a nondimensional format. An errant dynamic viscosity term is suspected of affecting the high-speed 
aerodynamic environment when executing models in the R-12 medium. 

To obtain the best possible results the scaled model for R-12 was converted to execute as a scaled model in air. 
Therefore, the tables reflect the inputs necessary to execute an air-scaled model of the BRR. The geometric scaling 
remained the same. However, the mass and stiffness scaling has been modified according to the relations 

m air = 0.396m R| 2 

EI,„ - 1.973EI R l2 

The mass scaling applies to mass moment of inertia scaling and the stiffness scaling applies to the torsional 
stiffness scaling as well. Because the geometric scaling is the same, displacement results are not scaled. However, 
forces and moments must be scaled according to the relations 

F r . 12 = 0.5054 

Mr- 12 = 0.5054M air 

Table Cl. BRR CAMRAD-II Shell Inputs 


! Basic Research Rotor base shell input deck 
1 

&NLDEF class= ' CASE ’ &END 

&NLVAL TITLE=' Basic Research Rotor’, 

OPUNIT= 1 , TMTASK= 1 , 

OPDENS=l , ALTMSL=0 . 

&END 

&NLDEF class- ' TRIM 1 &END 

&NLVAL VTIPIN=3 , MTIP=0.65, OPTRIM=0, MPSI=24, MPSIAV=4, 

DOFA= 6 * 0 , D0FD= 8*0, 

ScEND 

&NLDEF class= ’ TRIM ROTOR’, name=’ROTOR 1’ &END 

&NLVAL DOFG = 0 , DOFS=0, DOFM=7 *1 , 33 *2 , DOFB = 8*l , 4*0 , OPMODE = l, 

&END 

&NLDEF class= ' FLUTTER ' &END 

&NLVAL OPFLUT=0 , OPMEAN=l, OPBLD=0, OPSTAB^O, OPAERO=0, 

DOFA=6 * 0 , DOFD=8 * 0 , TASK=1,3*0, OPEIGN=l, MEIGN=2 , 

&END 

&NLDEF class= ' FLUTTER ROTOR ' , name =’ ROTOR 1' &END 
&NLVAL OPAERO^O, OPMODE=l, DOFM=7 * 1 , 33 *2 , OPTRAN=l, 

&END 

&NLDEF class= ' AIFRFRAME ’ , type= ’ STRUCTURE 1 &END 
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ScNLVAL TITLE= 1 ARES-2 Testbed - Simple Model ' , 

CONFIG=0, OPFREE= 0 , OPTRAN=0, OPAERO=0, 

WEIGHT=1 . 0 , MASSR=0.29984, IXXR=0., IYYR=0., IZZR=2 . 193752 , 

IXYR=0 . , IXZR=0 . , IYZR=0 . 

FSCG =-0.0500, BLCG =0., WLCG =6.17, 

FSRTR= 0.1167, BLRTR=0 . , WLRTR=8.01, 

ASHAFT=0 . , HSP=0 . 685 , 

ScEND 

&NLDEF class= ' AIRFRAME * , type= ' CONTROL ’ ScEND 
ScNLVAL K0 = 1 . 0 , KC = 1.0, KS=1.0, 

Sc END 


ScNLDEF class= ' ROTOR 1 , type= * STRUCTURE * , name= * ROTOR 1 1 ScEND 
ScNLVAL TITLE=' Basic Research Blades V.8 on Articulated Hub', 
RADIUS=4 . 685, NBLADE=4 , ROTATE =1 , SIGMA=0 . 0821 , 
CONFIG=l , OPAERO=l , GIMBAL=0 , CONTRL=2 , 

HINGE=3 , EFLAP=0 . 0534 , ELAG=0.0534, DLAG=15.96, 
PITCH=1 , LOCKP=l , EPITCH=0 . 08 , 

OPTPP=2 , ETPP=0 . 0534 , 

LOCKPL=0 , EPH=0 . 08894 , 

XSP=0. 03113, YSP=0 .01939, ZSP=-0 . 14621 , 

XPH=0. 02490, YPH=0. 04002, ZPH= 0., 

TWIN= 1 , TWISTL=0 . , 


&END 

ScNLDEF 


KNODE=4, RNODE=0 .2, 0.4, 0.6, 0.8, OPWING=l, OPPOS=l, NRPOS=0, 
GDAMPU= 0 . 01 , GDAMPV=0 . 01 , GDAMPW=0.01, GDAMPT=0.01 
KNODE=5 , RNODE=0 . 1221 , 0.2, 0.4, 0.6, 0.8, DRELST=0.07, 
NPROP=l 0 , 

RPROP= 0 . , 0.05339, 0.05341, 0.08, 0.1220, 0.1222, 

0.16, 0.18, 0.20, 1.0, 


KP=5*0. 02043, 2*0.01350, 0.00451, 2*0.00669, 

EIFLAP = 2 * 9 . 86 5E5 , 3*189244.0, 2*13179.0, 309.2, 

EILAG= 2*11 . 83 8E5 , 3*189244.0, 2*24051.0, 15145.0, 
GJ= 2*9. 8 65E5 , 

EA= 10*414330.0, 

MASS= 2*0.1188, 

ITHETA=2 * 0 . 003168, 


3*151394.0, 2*11178.0, 429.7, 


3*0.05576, 2*0.02895, 

3*0.000170, 2*0 . 000075 , 


0.00653, 

0.000012 


IPOLAR=2*0 . 003168, 3*0.000170, 2*0.000075, 0.000012 


2*197.3, 

2*5919.0, 

2*268.3, 


2*0.00634, 

, 2*0.000053, 
, 2*0.000053 


class= 'ROTOR* , type =’ AERODYNAMICS ’ , name =' ROTOR 1 • ScEND 


ScNLVAL OPSCEN=l , CHORD=0 . 3 02 1 , NPANEL = 19, REDGE= 0 .19, 0.28, 0.35, 

0.42,0.48,0.54,0.59,0.64,0.69,0.73,0.77,0.81,0.84,0.87, 

0.90,0.92,0.94,0.96,0.98,1.0 


ScEND 

ScNLDEF class= ' ROTOR ’ , type =• INFLOW ' , name =' ROTOR 1* ScEND 
ScNLVAL OPSCEN=l ScEND 


ScNLDEF class= 1 TABLES 1 ScEND 
ScNLVAL ScEND 


ScNLDEF action='end of shell 
ScNLDEF action='end of core' 


ScEND 

ScEND 
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Table C2. BRR CAMRAD-I1 Run-Time Inputs and Core Changes 

1 

I BRR model w/ hub motion 

t 

setenv BLADEAIRFOILl lr0012.tab 
setenv SHELLINPUT thesisl.dat 
camrad > thesisl.out << ’ endof input 1 
ScNLJOB NCASES = 1 ScEND 


ScNLDEF class= ' CASE 1 ScEND 
ScNLVAL TMTASK = 1, 

0PDENS=1 , ALTMSL=0 . 

ScEND 

ScNLDEF class= ' TRIM ' ScEND 
ScNLVALRELAXC = 0 - 1 , WINDING, WVEL=0.05, 

OPTRIM=l , MTRIM=4 , 

MNAME ( 1 ) = ' CL/S 1 , VNAME ( 1 ) = ' COLL ' , CLTRIM= 0.07065, 

MNAME { 2 ) = ' CX / S ' , VNAME < 2 ) = ’ PITCH 1 , CXTRIM--0 . 00 013 , 
MNAME ( 3 ) = ’ BETAS ’ , VNAME ( 3 ) = 1 LATCYC ’ , BSTRIM= 0.0, 

MNAME ( 4 ) = ' BETAC 1 , VNAME ( 4 ) = 1 LNGCYC ' , BCTRIM= 0.0, 
COLL-11 . 14 , PITCH=0 . 15 , LATCYC=-3 . 37 , LNGCYC=0.89, 
ITERU=1, ITERP=1 , NWPRNT=1 , TOLERT=0.5, LEVEL =3 , 

&END 

ScNLDEF class= 1 TRIM ROTOR’, name = ' ROTOR 1' ScEND 
ScNLVAL MHSEN=3 , MCSEN=3, MPSEN=3, 

GDAMPM=4 0*0 . 01 

ScEND 

ScNLDEF class= ' FLUTTER ' ScEND 
ScNLVAL OPEIGN= 0 , MPSIAV=1 ScEND 

ScNLDEF class= ' FLUTTER ROTOR', name = ’ ROTOR 1' ScEND 
ScNLVAL ScEND 

ScNLDEF class= ’AIRFRAME ' , type = ’ STRUCTURE 1 ScEND 
ScNLVAL ScEND 

ScNLDEF class= ' ROTOR ' , type = * STRUCTURE ’ , name = ’ ROTOR 1* ScEND 
ScNLVAL OPAERO=l, OPPOS = 4 , NRPOS=2 , 

ScEND 

ScNLDEF class= ' ROTOR’ , type = ’ AERODYNAMIC ' , name= ' ROTOR 1’ ScEND 
ScNLVAL ScEND 

ScNLDEF class= ' ROTOR ' , type=’WAKE', name =’ ROTOR 1' ScEND 
ScNLVAL OPFW=l, OPNW=l, OPLL=l, OPAX=0 , 

KNW=4 , KRU = 6 , KFW=96, KDW=96, 

OPRWG= 1 , FK2TWG=0 . 9 , TWIST-0., 

OPCORE=0 , CORE ( 1 ) = 0 . 2 , 

OPFWG=2 , MFWG- 6 , COREWG=0.2, MBWG=4 , DQWG=2 * 0 . 00 03 , 
ITERWG=2 , RLXWG=0 . 5 , RFWG=0 . 04 , 0 . 5 , RNWG=0 . 04 , 0 . 12 5 , 
RUDWG=0 . 5 , OPDWG=l , MPSIWG=24, OPDISP=3,l 

ScEND 


ScNLDEF action= ' end of shell' ScEND 


Begin CORE input 


! Set in proper hub mass and inertia properties 
! Add Hub Displacement and Hub Rotation sensors and output 
! Add Non-rotating force and moment sensors and outputs which include 
! hub mass 

ScNLDEF class= ’ COMPONENT 1 , type = ' RIGID BODY', name= ' ROTOR 1 HUB' ScEND 
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fitNLVAL MASS=0 . 663 , IXX=0.0377, IYY=0.0377, IZZ=0.0647, 

I X Y = 0 . , IXZ = 0. , IYZ = 0 . 

NSEN=8, KSEN= 8 , SENV ( 5 ) =5 , 6 , 7 , 8 , 

SNAME ( 5 ) = ' HUB DISPLACEMENT ’ , SLABL ( 5 ) = • HUB DISP ' , 

SNAME ( 6 ) = ' HUB ROTATION' , SLABL ( 6 ) = ' HUB ROT' , 

SNAME (7) = 'NONROTATING HUB FORCE (w/ HUB)' , 

SNAME ( 8 )=' NONROTATING HUB MOMENT (w/ HUB)’, 

KINDQ(5)=4*1, IDENTQ(5) =4*1, QUANT {5)=1,7,41,42, 

KINDR(5) =2*6, FRAMER ( 5 ) =2 * ' ' , SCALE { 5 ) =4 * 1 . 0 

St END 

fitNLDEF class= 'OUTPUT' , type= ' ' , name =' HUB DISPLACEMENT' &END 
ScNLVAL KINDY=4 , YNAMEC= ' ROTOR 1 HUB', YNAMEV= ' HUB DISPLACEMENT', 

TMPART= ’ROTOR 1 HUB LOAD SENSOR', 

TNPART= ’ROTOR 1 HUB LOAD SENSOR', 

FLPART= ' ROTOR 1 HUB ’ , 

NTPRNT = 1, NHPRNT=1 , MTIME=24, MH ARM = 10, OPHARM=l 

&END 

i 

&NLDEF class= ' OUTPUT ’ , type= ' ' , name= ' HUB ROTATION’ &END 
fitNLVAL KINDY=4, YNAMEC= ' ROTOR 1 HUB', YNAMEV= ’ HUB ROTATION', 

TMPART= ’ROTOR 1 HUB LOAD SENSOR' 

TNPART= ’ROTOR 1 HUB LOAD SENSOR’, 

FLPART= ’ ROTOR 1 HUB ' , 

NTPRNT=1, NHPRNT=1 , MTIME=24, MHARM= 10 , OPHARM=l 

Sc END 

i 

&NLDEF class= ’OUTPUT’ , type= ' ' , name =' NONROTATING HUB FORCE’ &END 
fitNLVAL KINDY=4, YNAMEC= ’ ROTOR 1 HUB', YNAMEV= ' NONROTATING HUB FORCE (w / HUB)', 
TMPART= 'ROTOR 1 HUB LOAD SENSOR', 

TNPART=’ ROTOR 1 HUB LOAD SENSOR', 

FLPART= ’ ROTOR 1 HUB ’ , 

NTPRNT=1 , NHPRNT=1 , MTIME=24, MHARM=10 , OPHARM=l 

fit END 

ficNLDEF class =' OUTPUT' , type= ' ' , name =' NONROTATING HUB MOMENT' ficEND 
fitNLVAL KINDY=4, YNAMEC= 1 ROTOR 1 HUB', YNAMEV= ' NONROTATING HUB MOMENT (w/ HUB)', 
TMPART= ' ROTOR 1 HUB LOAD SENSOR' , 

TNPART= ' ROTOR 1 HUB LOAD SENSOR ’ , 

FLPART= ' ROTOR 1 HUB ’ , 

NTPRNT=1 , NHPRNT=1 , MTIME=24, MHARM=10, OPHARM=l 

St END 

fitNLDEF class= ’ OUTPUT ' , type= ' * , name =' NONROTATING SWASHPLATE LOAD’ &END 
fitNLVAL KINDY=4 , YNAMEC= ' PYLON ' , YNAMEV= ' NONROTATING SWASHPLATE LOAD', 

TMPART= ' ROTOR 1 HUB LOAD SENSOR', 

TNPART= ' ROTOR 1 HUB LOAD SENSOR ' , 

FLPART= 1 ROTOR 1 HUB ’ , 

NTPRNT=1, NHPRNT=1 , MTIME=24, MHARM=10 , OPHARM=l 

fit END 


! Delete unused components 

fitNLDEF class =’ COMPONENT ' , type= ’ RIGID BODY', name =' AIRFRAME ROTOR 1', 
action= ’ delete ' 

fit END 


! Build PYLON ROTOR 1 component from AIRFRAME ROTOR 1 input 

&.NLDEF class= 'COMPONENT' , type= ' RIGID BODY', name =' PYLON ROTOR 1* &END 

fitNLVAL LABEL = ' PY RTR 1', DEBUG=0 , 

TMPART= ’AIRFRAME 1 , 
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TNPART= * ROTORCRAFT 1 , 

FLPART= * AIRFRAME 1 , 

FRAME= ’ AIRFRAME ’ , OPMTRX=0, 

KINDFC = 0 , KINDRM=2 , 

RNAME= ' PYLON ROTOR 1 RIGID BODY’, RLABL= ' P R 1 RB ' , 
JNTDOF = 0 , 

NLOC=l, CEREP=0, ZCJ=3*0., 

NJOINT=0, 

NCNXN=1, CCREP=0 , ZCJ=3*0., CNXJNT=0, CNXLOC=l, 
NSDINT=2 , SDICNX=1, 1, 

SDNAME ( 1 ) = ’ PYLON/ PYLON ROTOR 1 ' # SDLABL { 1 ) = ’ P /PR 1 1 , 
SDNAME (2) = ’ PYLON ROTOR 1/HUB 1', SDLABL { 2 )= 1 PRl /H 1’, 
MASS=0 . 

Sc END 


! Modify existing rigid body components 

&NLDEF class =' COMPONENT ’ , type= ' RIGID BODY’, name = ’ ROTOR 1 HUB’ ScEND 
&NLVAL SDNAME ( 1 ) = ' PYLON ROTOR 1 /HUB 1 1 , SDLABL ( 1 ) = ' P /H 1 ' ScEND 

i 

&NLDEF c lass =' COMPONENT ’ , type=’RIGID BODY’, name='ROTOR 1 SWASHPLATE ' 
&NLVAL SDNAME ( 1) =’ PYLON/ SWASHPLATE 1*, SDLABL = 1 P /SP 1’ &END 


! Add PYLON component 

&NLDEF c lass = ’COMPONENT ' , type= * RIGID BODY ’ , name^’ PYLON' ScEND 
&NLVAL LABEL = ' PYLON ’ , DEBUG =0 , 

TMPART= ’ AIRFRAME ’ , 

TNPART= ' ROTORCRAFT ' , 

FLPART= ’ AIRFRAME ’ , 

FRAME= 1 AIRFRAME ’ f OPMTRX=0 , 

KINDFC= 0 , KINDRM = 2 , 

RNAME= ’ PYLON RIGID BODY ' , RLABL= ' PYLON R’ , 

JNTDOF=6, 

NLOC=3 , CEREP=3*0, ZEB= 0., 0., 0., 

0., 0., 1.0067, 

0., 0., 0.3217, 

NJOINT=l, CJREP=0, ZJE=3*0 . , JNTLOC=l, KINDJ-3, KINDJL=3, KINDJA: 
AXISL = 1 ,2,3, GAINL=3 *1.0 , TYPEL=3*1, JNTEL=1,2,3, 

AXISA^l ,2,3, GAINA=3*0. 0174533, TYPEA-3 *1 , JNTEA=4,5,6, 

NSPRNG=6 , KTYPE = 6 *1 , CTYPE=6*1, 

KLIN=3 *1 . E2 0 , 3 *1 . E8 , KEQUIV=3 *1 . E2 0 , 3 *1 . E8 , 

CLIN=6 *0 . , CEQUIV=6 *0 . , 

OFFSET = 6* 0 . , BIAS = 6 * 0 . , OPTENS = 0, 

NCNXN=3 , CCREP=3 *0 , ZCJ=9*0., CNXJNT=1 , 2 *0 , CNXLOC=l , 2 , 3 , 

NSDINT=3 , SDICNX=1 ,2,3, TORQ=3*0, 

SDNAME (1)=’ AIRFRAME /PYLON' , SDLABL ( 1 )=' AF/P 1 , 

SDNAME ( 2 ) = ’ PYLON/ PYLON ROTOR 1 ' , SDLABL ( 2 ) = ’ P /PR 1 ’ , 

SDNAME (3) =’ PYLON/ SWASHPLATE 1’, SDLABL (3 )=’P /SP 1', 

NSEN=2, OPDEFN=2 *1 , 

SNAME ( 1 ) = ' HUB MOTION SENSOR’, SLABL ( 1 ) = 1 HUB SEN', 

SNAME ( 2 ) = ’ NONROTATING SWASHPLATE LOAD ' , SLABL ( 2 ) = ' NR SP LD ' , 


LENS =6 , 6, 

KSEN=12 , SENV=6*1, 6*2, SENE=1 , 2 , 3 , 4 , 5 , 6 , 1 , 2 , 3 , 4 , 5 , 6 , 
KINDQ=6*2, 6*1, IDENTQ=1, 2 ,3, 4, 5, 6,6*3, 

QUANT=6 *21 ,3*41,3*42, AXIS=1 , 2 , 3 , 1 , 2 , 3 , 1 , 2 , 3 , 1 , 2 , 3 , 
SCALE=3 *1.0, 3*0. 0174533, 6*1.0, OPSCL=12*0, 


SELABL ( 1 ) = ’ HUB LX ’ , 
SELABL ( 2 ) = ’ HUB LY ' , 
SELABL ( 3) = ’HUB LZ ’ , 
SELABL ( 4 ) — ' HUB AX 1 , 
SENAME ( 5 ) = ’ HUB MOTION AY’, SELABL ( 5 ) = ’ HUB AY', 


SENAME ( 
SENAME ( 
SENAME ( 


1 ) 

2 ) 

3) 


SENAME ( 4 ) = 


• HUB MOTION LX 1 
' HUB MOTION LY ' . 
’HUB MOTION LZ' 
’HUB MOTION AX' 


&END 
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SENAME ( 6) =' HUB MOTION AZ ' , SELABL ( 6 ) = ’ HUB AZ ’ , 

SENAME ( 7)=’ NR SP LOAD LX' , SELABL ( 7)=’ NR SP LX ' , 

SENAME ( 8)=' NR SP LOAD LY 1 , SELABL ( 8)=’ NR SP LY ' , 

SENAME ( 9)=' NR SP LOAD LZ 1 , SELABL ( 9 ) = ’ NR SP LZ 1 , 

SENAME ( 10) = ’NR SP LOAD AX' , SELABL ( 1 0 ) = ' NR SP AX', 

SENAME ( 11 )= 'NR SP LOAD AY 1 , SELABL { 11 ) = 1 NR SP AY 1 , 

SENAME { 12 )= 'NR SP LOAD AZ ’ , SELABL { 12 )=’ NR SP AZ ' , 

NCON=l , CNAME= 1 HUB MOTION INPUT’, CLABL= ' HUB MOTN' , LENC=6, 

OPACT = 6 *2 , CONVKC=6*l, CONEKC=l, 2 , 3, 4 , 5 , 6, 

MASS=0 . 

& END 

StNLDEF class= ' RESPONSE 1 , type= ' RIGID ' # name =' DOF PYLON RIGID' StEND 

StNLVAL KINDR=1 , 

RNAMEC= ' PYLON’ , RNAMEV= ’ PYLON RIGID BODY * , 

OPFILT=l, KINDRM=2 , REPLIN=1, RE PANG =2 , 

WTCLIN= . 05 , WTCANG=.02, WTPLIN=.05, WTPANG=.02, 

KINDN=3 , KINDDN=1, CNREP=4 , 

CNOM (1,1)= -1. , 0. # 0. # 

CNOM (1,2)= 0. f 1 . , 0 . , 

CNOM (1,3)= 0., 0., -1., 

ZNOM=-0 . 1667 , 0., -0.8333, 

QREF = 3*0, QINIT=3 *0 , QOPSYM=3*0, QTMRED=3*1, QTNRED=3*1, QFLRED=3*1 , 
PREF = 3*0, PINIT=3*0 , POPSYM=3*0, PTMRED=3*1, PTNRED=3*1, PFLRED=3*1, 

StEND 


! Modify AIRFRAME normal modes component 

StNLDEF c 1 as S= ' COMPONENT ' , type= ' LINEAR NORMAL MODES', name =' AIRFRAME ’ StEND 
StNLVAL ZEB (1,1) =-0 .1667, 0., -0.8333, 

NSDINT=l f SDICNX=1, 

SDNAME (1) = 'AIRFRAME /PYLON' , SDLABL ( 1 ) = ' AF/P 1 

& END 


! Delete old structural dynamic interfaces 

StNLDEF class= ' INTERFACE' , type='SD', name = 1 AIRFRAME /AIRFRAME ROTOR 1', 
action^ 'delete ' 

& END 

StNLDEF class= ' INTERFACE ' , type= ' SD ’ , name= ' AIRFRAME ROTOR 1 /HUB 1 ' , 
action^ ' delete ' 

St END 

StNLDEF class= ' INTERFACE’ , type= ' SD* , name = ' AIRFRAME/ SWASH PLATE 1 ’ , 
action= ' delete ’ 

St END 


I Build new structural dynamic interfaces 

StNLDEF class= * INTERFACE ' , type='SD', name =' AIRFRAME /PYLON ' StEND 
StNLVAL LABEL= ’ AF/P ' , DEBUG=0 , 

TMPART= 'AIRFRAME ’ , 

TNPART= ' ROTORCRAFT ’ , 

FLPART= ' AIRFRAME ' , 

AN AMEC = ' AIRFRAME ’ , ANAMEV= ' AIRFRAME/ PYLON ' , 

BNAMEC= ' PYLON ' , BNAME V = ' AIRFRAME /PYLON ' , 

KINDSD=1 , AXESC=0 , OPELIM=l, NXELIM=0 

StEND 

StNLDEF class= ' INTERFACE’ , type='SD' f name =' PYLON/ PYLON ROTOR 1* StEND 
StNLVAL LABEL= ' P /PR 1’, DEBUG=0 , 

TMPART= ’ AIRFRAME ’ , 

TNPART= ’ ROTORCRAFT ’ , 

FLPART- ’ AIRFRAME ’ , 

ANAMEC= ’ PYLON' , ANAMEV= ' PYLON/ PYLON ROTOR 1', 
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BNAMEC= ' PYLON ROTOR 1', BNAMEV= ' PYLON/ PYLON ROTOR 1', 
KINDSD=1, AXESC=0 , OPELIM=l, NXELIM-0 

& END 


ScNLDEF class=' INTERFACE' , type='SD', name - 1 PYLON ROTOR 1/HUB 1’ St. END 
ScNLVAL LABEL=’PR 1/H 1', DEBUG=0 , 

TMPART= ' ROTOR 1 HUB ' , 

TNPART= ' ROTORCRAFT ’ , 

FLPART= ’ ROTOR 1 HUB 1 , 

ANAMEC= ' PYLON ROTOR 1', ANAMEV= ' PYLON ROTOR 1/HUB 1', 

BNAMEC= ‘ROTOR 1 HUB’ , BNAMEV= ’ PYLON ROTOR 1/HUB l‘ f 
KINDSD=1, AXESC=0 , OPELIM=l, NXELIM-0 

& END 


ScNLDEF class= ' INTERFACE ' , type='SD', name =' PYLON/ SWASH PLATE 1‘ ScEND 
ScNLVAL LABEL- * P /SP 1 ' , DEBUG=0, 

TM P ART = 'ROTOR 1 HUB 1 , 

TNPART= ’ ROTORCRAFT ' , 

FLPART= * ROTOR 1 HUB 1 , 

ANAMEC= ' PYLON ' , ANAMEV= ’ PYLON /SWASH PLATE 1*, 

BNAMEC= 'ROTOR 1 SWASH PLATE ’ , BNAMEV= ' PYLON/ SWASH PLATE 1’, 
KINDSD=1, AXESC=0 , OPELIM=l, NXELIM=0 

ScEND 


! Delete old RESPONSES 

ScNLDEF class= ‘RESPONSE’ , type = ’ RIGID ' , name = ' DOF AIRFRAME ROTOR 1 RIGID BODY', 


ScEND 

action^ ' delete ’ 



ScNLDEF 

ScEND 

class= ' RESPONSE ‘ , 
action^ ' delete 1 

type= 

•VARIABLE’ 

ScNLDEF 

ScEND 

class= 'RESPONSE' , 
action^ ' delete 1 

type = 

•VARIABLE’ 

&NLDEF 

ScEND 

class= ’RESPONSE’ , 
action^ ' delete 1 

type= 

•VARIABLE' 


name= ’ INT AIRFRAME /SWASH PLATE 1 ' , 
name= 1 INT AIRFRAME /AIRFRAME ROTOR 1 ‘ , 
name= 1 INT AIRFRAME ROTOR 1/HUB l' # 


! Create new RESPONSES 

ScNLDEF class= ' RESPONSE 1 , type = ' RIGID ' , name =' DOF PYLON ROTOR 1 RIGID BODY' 
ScEND 

ScNLVAL DEBUG=0 , KINDR=1, 

RNAMEC= ' PYLON ROTOR 1 ’ , 

RNAMEV= 1 PYLON ROTOR 1 RIGID BODY', 

OPFILT=l, KINDRM=2 , REPLIN=1, RE PANG =2 , 

WTCLIN= . 05 , WTCANG=.02, WTPLIN=.05, WTPANG=.02, 


KINDN=3 , KINDDN= 

: 1 , 

CNREP=4 

CNOM (1,1)= -1. , 

0 

. , 0 . , 

CNOM (1,2)= 0 . , 

1 

. , 0 . , 

CNOM (1,3)= 0 . , 

0 

- 1 ., 

ZNOM=-0 .1667, 0. 


-1.84, 


QREF = 3 *0 , QINIT = 3*0, QOPSYM=3*0, QTMRED=3*1, QTNRED=3 *1 , QFLRED=3 * 1 , 
PREF = 3 *0 , PINIT=3*0 , POPSYM=3*0, PTMRED=3*1, PTNRED=3*1, PFLRED=3*1, 

ScEND 

i 

ScNLDEF class= ’ RESPONSE ' , type= 1 VARIABLE ’ , name = ' INT PYLON/ SWASHPLATE 1' ScEND 
ScNLVAL DEBUG = 0 / KINDR=2, 

RNAMEC= ’ ' , 

RNAMEV= ’ PYLON/ SWASHPLATE 1 ’ , 

OPFILT-O, MVAR=6 , KINDN=3, 
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& END 
1 

ScNLDEF 

ScNLVAL 


ScEND 

i 

&NLDEF 

ScNLVAL 


ScEND 


! Hub 
ScNLDEF 
ScNLVAL 


&END 

i 

ScNLDEF 

ScNLVAL 


ScEND 

i 

ScNLDEF 

ScNLVAL 


ScEND 

i 

ScNLDEF 

ScNLVAL 


ScEND 


XN0M=6 * 0 . , WTCONV =3*1. 3, 3*0. 6, WTPERT=3 *1 . 3 , 3 *0 . 6 , WTMODE=6*0. # 
XREF=6*0 . , XINIT=6*0 . , TMRED=6 *1 , TNRED=6 * 1 , FLRED=6*1, OPSYM=6*0 


class= 'RESPONSE ' , type= 'VARIABLE 1 , name= ' INT PYLON/PYLON ROTOR 1' ScEND 
DEBUG=0 , KINDR=2, 

RNAMEC= * 1 , 

RNAMEV= * PYLON/ PYLON ROTOR l 1 , 

OPFILT=0, MVAR=6 , KINDN=3, 

XNOM= 6*0., WTCONV =3*1. 3,3*0. 6, WTPERT=3 *1 . 3 , 3 *0 . 6 , WTMODE=6*0., 

XREF = 6*0., XINIT=6*0., TMRED=6 * 1 , TNRED=6*1 , FLRED=6*1, OPSYM=;6*0 


class= ' RESPONSE ' , type = 'VARIABLE ’ , name = • INT PYLON ROTOR 1/HUB 1' ScEND 
DEBUG =0 , KINDR=2, 

RNAMEC= * ' , 

RNAMEV= * PYLON ROTOR 1/HUB 1', 

OPFILT=0, MVAR=6 , KINDN=3 , 

XNOM=6 * 0 . , WTCONV= 3 * 1 . 3 , 3 * 0 . 6 , WTPERT=3 *1 . 3 , 3 *0 . 6 , WTMODE=6*0. , 
XREF=6*0., XINIT=6*0., TMRED=6*1, TNRED=6*1, FLRED=6*1, OPSYM=6*0 


motion input /output and connections 

class=* OUTPUT’ , type= * * , name=’HUB MOTION SENSOR* ScEND 
K INDY =4 , YNAMEC= * PYLON 1 , YNAMEV= * HUB MOTION SENSOR ' , 
TMPART= 1 ROTOR 1 HUB LOAD SENSOR * , 

TNPART= 'NO SOLUTION' , 

FLPART- 1 NO SOLUTION ’ , 

NTPRNT=1, NHPRNT=1 , MTIME=24, MHARM=10 , OPHARM=l 


class= 1 RESPONSE * , type =* VARIABLE ' , name =* OUT HUB MOTION SENSOR' ScEND 
KINDR=4 , RNAMEV= ' HUB MOTION SENSOR', 

MVAR-6 , KINDN=3 , XNOM=6*0. , WTMODE=6*0., XREF^6*0., XINIT=6*0., 
OPSYM=6 *0 , WTCONV=3 * 0 . 05,3*1.0, WTPERT=3 *0 . 05 , 3 *1 . 0 , 

TMRED=6 * 1 , TNRED= 6*1, FLRED=6 *1 


class= ' INPUT’ , type= * 1 , name=*HUB MOTION LX HHC INPUT* ScEND 
LABEL= ' LX HHC*, NCNXN=1, 

FNAMEC= ’ HUB MOTION OSC AMP LX 1 , 

FNAMEV = ’ HUB MOTION LX HHC INPUT ' , 

LENGTH=17 , 


UENAME= ' MEAN ’ , 

*1C* , 

'IS* , 

* 2C ' , 

’ 2S* , 

' 3C * , 

*3S* , 

' 4C ' , 

* 4S 


' 5C ' , 

' 5S ' , 

* 6C* , 

’ 6S* , 

1 7C ' , 

' 7S ' , 

*8C' , 

' 8S 

UELABL= ’MEAN* , 

*1C* , 

* IS’ , 

' 2C ’ , 

' 2S ' , 

' 3C * , 

' 3S ' , 

* 4C ' , 

* 4S 


* 5C * , 

* 5S * , 

*6C* , 

' 6S* , 

* 7C ' , 

1 7S ’ , 

*8C* , 

’ 8S 


class= ' RESPONSE ’ , type='VAR', name= ' INP HUB MOTION LX HHC INPUT' ScEND 
KINDR=3, RNAMEV = 1 HUB MOTION LX HHC INPUT’, 

MVAR=17 , KINDN=3 , XNOM=17*0., WTMODE=17 * 0 . , XREF=17*0., 

Axial hub motions are set in here — MEAN, 1PC, IPS, 2 PC, 2 PS, 

XINIT=0 . , 6*0., 0., 0., 6*0., 0., 0., 

OPSYM=17* 0 , WTCONV=17 *0 . 05 , WTPERT = 17 *0 . 05 , 

TMRED=17*1 , TNRED=17*1 # FLRED=17*1 
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ScNLDEF 

ScNLVAL 


ScEND 

i 

ScNLDEF 

ScNLVAL 


ScEND 

i 

ScNLDEF 

ScNLVAL 


&. END 

I 

ScNLDEF 

ScNLVAL 


& END 

ScNLDEF 

ScNLVAL 


& END 

i 

ScNLDEF 

ScNLVAL 


Sc END 

ScNLDEF 

ScNLVAL 


c las s= ' COMPONENT 1 , type =’ FOURIER ' , name = ' HUB MOTION OSC AMP LX' ScEND 
LAB EL = ' HUB MOTN ' , MHARM=8 , NBHARM=0, NEHARM=8 , TYPEA=2 , 

PERIOD- ' ROTOR 1', AZREF=0 . , GAIN=1.0 # OPDRVF=0, OPDRVX=0, 

F0NAME = ' HUB MOTION LX HHC INPUT • , F0LABL= * LX HHC', 

XNAME = ' HUB MOTION OSC AMP LX' , XLABL= ' LX HHC' 


class= ‘ INTERFACE ' , type='IO', name='HUB MOTION OSC AMP LX' ScEND 
LAB EL = ' AMP LX 1 , 

TMPART= ’ROTOR 1 CONTROL’, TMPASS=1, 

TNPART= ’ NO SOLUTION 1 , 

FLPART= ' NO SOLUTION 1 , 

XNAMEC = ’ HUB MOTION OSC AMP LX ' , 

XNAMEV= ' HUB MOTION OSC AMP LX ' , 

NCNXN = 1 , 

FNAMEC= ' HUB MOTION * , FNAMEV = * HUB MOTION OSC AMP LX' 


class= * RESPONSE 1 , type= 1 VARIABLE ' , name = ' INT HUB MOTION OSC AMP LX* ScEND 
KINDR=2 , RNAMEV = * HUB MOTION OSC AMP LX*, 

MVAR=1 , KINDN=3 , XNOM=0., WTMODE=0 . , XREF=0., XINIT=0., 

OPSYM=0 f WTCONV=0 . 05 , WTPERT=0.05, 

TMRED=1, TNRED=1 , FLRED=1 


class=* INPUT* , type= * ' , name='HUB MOTION LY HHC INPUT’ ScEND 
LABEL = ’ LY HHC ’ , NCNXN=1 , 

FNAMEC= * HUB MOTION OSC AMP LY*, 

FNAMEV = ’ HUB MOTION LY HHC INPUT’ 

LENGTH =17 , 


UENAME= ’MEAN ' , 

'1C , 

’IS’ , 

* 2C 1 , 

' 2S 1 , 

' 3C* , 

' 3S ' , 

* 4C ' , 

' 4S 


’ 5C* , 

’ 5S ’ , 

' 6C' , 

' 6S’ , 

1 7C' , 

' 7S ' , 

' 8C' , 

' 8S 

UELABL= 1 MEAN ’ , 

' 1C , 

' IS ' , 

' 2C ' , 

’ 2S ' , 

’ 3C ’ , 

' 3S ' , 

' 4C ’ , 

' 4S 


’ 5C' , 

* 5S ’ , 

' 6C’ , 

* 6S ' , 

' 7C ' , 

* 7S ’ , 

' 8C’ , 

’ 8S 


class= ' RESPONSE ’ , type= ’VAR * , name= * INP HUB MOTION LY HHC INPUT' ScEND 
KINDR=3 , RNAMEV='HUB MOTION LY HHC INPUT', 

MVAR=17 , KINDN=3 , XNOM=17*0., WTMODE=17 *0 . , XREF=17*0., 

Side hub motions are set in here -- MEAN, 1PC, IPS, 2 PC, 

XINIT=0 . , 6*0., 0., 0., 6*0., 0., 0., 

OPSYM=17 *0 , WTCONV=17*0 .05, WTPERT = 17*0 .05, 

TMRED=17 *1 , TNRED=17*1, FLRED=17*1 


class= 'COMPONENT' , type= * FOURIER ’ , name='HUB MOTION OSC AMP LY ' ScEND 
LABEL= 'HUB MOTN* , MHARM= 8 , NBHARM= 0 , NEHARM= 8 , TYPEA=2 , 

PERI OD= ’ ROTOR 1*, AZREF=0 . , GAIN=1.0, OPDRVF=0, OPDRVX=0, 

F0NAME= ' HUB MOTION LY HHC INPUT', F0LABL= ' LY HHC', 

XNAME = * HUB MOTION OSC AMP LY ' , XLABL= ' LY HHC' 


class= 1 INTERFACE 1 , type=’IO', name='HUB MOTION OSC AMP LY ’ ScEND 
LABEL = ' AMP LY ’ , 

TMPART= ' ROTOR 1 CONTROL’, TMPASS=1, 

TNPART= ' NO SOLUTION ' , 

FLPART= ’ NO SOLUTION', 

XNAMEC = ’ HUB MOTION OSC AMP LY 1 , 


2 PS , 
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ScEND 

j 

ScNLDEF 

ScNLVAL 


ScEND 

j 

ScNLDEF 

ScNLVAL 


& END 

ScNLDEF 

ScNLVAL 


Sc END 

I 

ScNLDEF 

ScNLVAL 


Sc END 

i 

ScNLDEF 

ScNLVAL 


Sc END 

i 

ScNLDEF 

ScNLVAL 


ScEND 


XNAMEV= ' HUB MOTION OSC AMP LY 1 , 

NCNXN=1 , 

FNAMEC= ' HUB MOTION' , FNAMEV= ' HUB MOTION OSC AMP LY ' 


class= ’RESPONSE’ , type = ' VARIABLE ’ , name= ' INT HUB MOTION OSC AMP LY 1 ScEND 
KINDR=2, RNAMEV= ’ HUB MOTION OSC AMP LY ’ , 

MVAR= 1 , KINDN=3 , XNOM=0., WTMODE=0., XREF=0., XINIT=0. f 
OPSYM=0 , WTCONV=0.05, WTPERT=0.05, 

TMRED=1 , TNRED=1 / FLRED=1 


class= ' INPUT' # type= ,, # name =' HUB MOTION LZ HHC INPUT' ScEND 
LABELS LZ HHC, NCNXN=1, 

FNAMEC= ' HUB MOTION OSC AMP LZ ' , 

FNAMEV= ' HUB MOTION LZ HHC INPUT ’ , 

LENGTH=17, 


UENAME= ' MEAN ' , 

•1C’ , 

'IS' , 

’ 2C ' , 

' 2S ' , 

' 3C ' , 

' 3S ' , 

' 4C’ , 

' 4S 


' 5C’ , 

' 5S' , 

’ 6C' , 

' 6S ' , 

' 7C ' , 

' 7S ' , 

' 8C’ , 

1 8S 

UELABL= ' MEAN ’ , 

’ 1C’ , 

'IS' , 

' 2C ' , 

* 2S ' , 

1 3C' , 

' 3S ’ , 

’ 4C' , 

' 4S 


' 5C ' , 

' 5S ' , 

1 6C' , 

’ 6S’ , 

' 7C ’ , 

’ 7S ' , 

’ 8C’ , 

1 8S 


class= 1 RESPONSE ’ , type=’VAR', name = ' INP HUB MOTION LZ HHC INPUT’ ScEND 
KINDR=3 , RNAMEV= ’ HUB MOTION LZ HHC INPUT’, 

MVAR= 17 , KINDN=3 , XNOM=17*0., WTMODE =17*0. , XREF=17*0., 


Normal hub motions are set in here -- MEAN, 1PC, IPS, 2 PC, 

XINIT=0., 6*0., 0., 0., 6*0., 0., 0., 

OPSYM-17 *0 , WTCONV =17*0.05, WTPERT = 17*0 . 05 , 

TMRED=17*1 , TNRED=17 * 1 , FLRED=17*1 


class= ' COMPONENT ' , type= ’ FOURIER ' , name=’HUB MOTION OSC AMP LZ ’ ScEND 
LABEL = ' HUB MOTN ' , MHARM=8 , NBHARM= 0 , NEHARM-8 , TYPEA=2 , 

PERIOD^ 'ROTOR 1’, AZREF=0 . , GAIN=1 . 0 , OPDRVF=0, OPDRVX=0, 

F0NAME= ' HUB MOTION LZ HHC INPUT’, F0LABL= ' LZ HHC 1 
XNAME = ’ HUB MOTION OSC AMP LZ ’ , XLABL= ’ LZ HHC’ 


class= ' INTERFACE’ , type= ’ IO 1 , name='HUB MOTION OSC AMP LZ ' ScEND 
LABEL= ' AMP LZ ' , 

TMPART= ' ROTOR 1 CONTROL', TMPASS=1, 

TNPART= 'NO SOLUTION’ , 

FLPART= ' NO SOLUTION’, 

XNAMEC= ' HUB MOTION OSC AMP LZ ’ , 

XNAMEV = ' HUB MOTION OSC AMP LZ ' , 

NCNXN=1, 

FNAMEC= ' HUB MOTION', FNAMEV= ’ HUB MOTION OSC AMP LZ ' 


class= ’RESPONSE’ , type= ’ VARIABLE ' , name =’ INT HUB MOTION OSC AMP LZ ' ScEND 
KINDR=2 , RNAME V = ' HUB MOTION OSC AMP LZ ' , 

MVAR=1, KINDN=3, XNOM=0 . , WTMODE=0., XREF=0 . , XINIT=0., 

OPSYM=0 , WTCONV=0 .05, WTPERT=0.05, 

TMRED=1 , TNRED=1 , FLRED=1 


2 PS, 
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ScNLDEF 

ScNLVAL 


ScEND 

i 

ScNLDEF 

ScNLVAL 


& END 

ScNLDEF 

ScNLVAL 

ScEND 

I 

ScNLDEF 

ScNLVAL 


ScEND 

t 

ScNLDEF 

ScNLVAL 

ScEND 

ScNLDEF 

ScNLVAL 


ScEND 


class= ' INPUT ' , type= * ’ , name=’HUB MOTION AX HHC INPUT' ScEND 
LABEL = ' AX HHC ' , NCNXN=1, 

FNAMEC = ' HUB MOTION OSC AMP AX', 

FNAMEV = ' HUB MOTION AX HHC INPUT ' , 

LENGTH=17 / 


UENAME= ' MEAN ’ , 

' 1C’ , 

' IS 1 , 

'2C', 

’ 2S 1 , 

*3C' , 

' 3S ’ , 

' 4C ' , 

’ 4S ' 


' 5C ’ , 

’ 5S ' , 

’ 6C ' , 

' 6S ' , 

' 7C ' , 

' 7S ’ , 

' 8C ' , 

' 8S ' 

UELABL- ’ MEAN ’ , 

'1C’ , 

'IS' , 

' 2C 1 , 

' 2S 1 , 

' 3C ' , 

' 3S ' , 

' 4C' , 

' 4S ' 


' 5C * , 

’ 5S ’ , 

’ 6C' , 

' 6S ' , 

' 7C ’ , 

' 7S ' , 

’ 8C' , 

' 8S ' 


class= ' RESPONSE ’ , type='VAR', name= ' INP HUB MOTION AX HHC INPUT' ScEND 
KINDR-3 , RNAMEV = ' HUB MOTION AX HHC INPUT’, 

MVAR=17, KINDN-3 , XNOM=17*0., WTMODE^17 *0 . , XREF=17*0., 

Roll hub motions are set in here — MEAN, 1PC, IPS, 2 PC, 
Must input some Side motion if want pure roll at hub 

XINIT=0., 6*0., 0., 0., 6*0., 0., 0., 

OPSYM=17 * 0 , WTCONV =17*0.05, WTPERT = 17*0 . 05, 

TMRED=17*1, TNRED=17 *1 , FLRED=17*1 


class= ' COMPONENT ' , type= ’ FOURIER ' , name='HUB MOTION OSC AMP AX’ ScEND 
LAB EL = ’ HUB MOTN 1 , MHARM-8, NBHARM= 0 , NEHARM= 8 , TYPEA=2 , 

PERIOD- ' ROTOR 1', AZREF=0 . , GAIN=1.0, OPDRVF=0, OPDRVX=0, 

F0NAME= ' HUB MOTION AX HHC INPUT’, F0LABL= ’ AX HHC’, 

XNAME= * HUB MOTION OSC AMP AX’ , XLABL= 1 AX HHC’ 


class- 1 INTERFACE ' , type-' IO', name='HUB MOTION OSC AMP AX' ScEND 
LABEL= 1 AMP AX ’ , 

TMPART= 1 ROTOR 1 CONTROL', TMPASS-1, 

TNPART= ' NO SOLUTION ’ , 

FLPART= 'NO SOLUTION' , 

XNAMEC = 1 HUB MOTION OSC AMP AX ' , 

XNAMEV- 1 HUB MOTION OSC AMP AX ' , 

NCNXN=1, 

FNAMEC = ' HUB MOTION’, FNAMEV =’ HUB MOTION OSC AMP AX' 


class= ' RESPONSE ’ , type= ' VARIABLE ' , name= ' INT HUB MOTION OSC AMP AX’ ScEND 
KINDR=2 , RNAMEV = ' HUB MOTION OSC AMP AX', 

MVAR= 1 , KINDN=3 , XNOM=0 . , WTMODE=0., XREF=0., XINIT = 0 . , 

OPSYM=0 , WTCONV=0 . 05 , WTPERT=0.05, 

TMRED=1, TNRED=1 , FLRED=1 


class= ’ INPUT ' , type= ’ ' , name=’HUB MOTION AY HHC INPUT’ ScEND 
LABEL = ' AY HHC ' , NCNXN=1 , 

FNAMEC = ’ HUB MOTION OSC AMP AY', 

FNAMEV = 1 HUB MOTION AY HHC INPUT’, 

LENGTH=17 , 


UENAME= 

' MEAN ’ , 

’1C , 

'IS' , 

' 2C , 

' 2S ' , 

’ 3C ' , 

' 3S ' , 

’4C , 

' 4S ' 



' 5C' , 

' 5S ' , 

•6C , 

' 6S ’ , 

’ 7C ' , 

' 7S ' , 

’80 , 

' 8S 1 

UELABL= 

' MEAN ’ , 

’ 1C’ , 

' IS’ , 

' 2C , 

' 2S ' , 

' 3C ' , 

’3S' , 

'40, 

' 4S ’ 



' 5C' , 

' 5S’ , 

’ 6C ' , 

' 6S ' , 

•7C' , 

’ 7S’ , 

' 8C' , 

’ 8S ' 


2PS, 
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&NLDEF 

ScNLVAL 


&END 

i 

ScNLDEF 

ScNLVAL 


ScEND 

i 

ScNLDEF 

ScNLVAL 


ScEND 

[ 

ScNLDEF 

ScNLVAL 


ScEND 

i 

ScNLDEF 

ScNLVAL 


&END 

j 

ScNLDEF 

ScNLVAL 


ScEND 


class=' RESPONSE ' , type= ' VAR 1 , name= ' INP HUB MOTION AY HHC INPUT’ ScEND 
KINDR=3 , RNAMEV= ' HUB MOTION AY HHC INPUT', 

MVAR=17 , KINDN=3, XNOM=17*0., WTMODE=17 * 0 . , XREF=17*0., 

Pitch hub motions are set in here -- MEAN, 1PC, IPS, 2 PC 
Must input some Axial motion if want pure roll at hub 

XINIT=0 . , 6*0., 0., 0., 6*0., 0., 0., 

OPSYM=17 * 0 , WTCONV =17*0 .05, WTPERT = 17*0 . 05 , 

TMRED=17 *1 , TNRED=17 *1 , FLRED=17*1 


class =’ COMPONENT' , type= ' FOURIER ' , name=’HUB MOTION OSC AMP AY’ ScEND 
LABEL= ' HUB MOTN ' , MHARM=8 , NBHARM=0 , NEHARM=8 , TYPEA=2 , 

PERIOD= 'ROTOR 1', AZREF=0 . , GAIN=1.0, OPDRVF^O, OPDRVX=0, 

F0NAME= ’ HUB MOTION AY HHC INPUT ' , F0LABL=’AY HHC', 

XNAME= ' HUB MOTION OSC AMP AY’ , XLABL= 1 AY HHC' 


class= ’ INTERFACE’ , type=’IO', name =’ HUB MOTION OSC AMP AY’ ScEND 
LABEL= ' AMP AY ' , 

TMPART= ’ROTOR 1 CONTROL', TMPASS=1, 

TNPART= ' NO SOLUTION’, 

FLPART= ' NO SOLUTION ’ , 

XNAMEC= ' HUB MOTION OSC AMP AY', 

XNAMEV= ’ HUB MOTION OSC AMP AY ' , 

NCNXN = 1 , 

FNAMEC= ' HUB MOTION', FNAMEV= ' HUB MOTION OSC AMP AY' 


class= ' RESPONSE ’ , type=‘ VARIABLE' , name = ’ INT HUB MOTION OSC AMP AY' ScEND 
KINDR=2 , RNAMEV='HUB MOTION OSC AMP AY', 

MVAR=1, KINDN=3 , XNOM=0., WTMODE=0 . , XREF=0., XINIT=0., 

OPSYM-0, WTCONV=0 .05, WTPERT=0.05, 

TMRED=1, TNRED-1, FLRED=1 


class= 1 INPUT ' , type= ' 1 , name='HUB MOTION AZ HHC INPUT 1 ScEND 
LABEL= 1 AZ HHC ’ , NCNXN=1 , 

FNAMEC- ' HUB MOTION OSC AMP AZ 1 , 

FNAMEV= ' HUB MOTION AZ HHC INPUT ’ , 

LENGTH=17 , 


UENAME= 'MEAN' , 

’1C , 

' IS’ , 

' 2C , 

1 2S ' , 

’ 3C’ , 

’ 3S ' , 

' 4C' , 

' 4S 


' 5C , 

' 5S ’ , 

' 6C , 

' 6S ' , 

'7C , 

' 7S 1 , 

' 8C , 

1 8S 

UELABL= ' MEAN ’ , 

’1C , 

' IS’ , 

' 2C , 

’ 2S ' , 

' 3C , 

’ 3S ' , 

' 4C 1 , 

' 4S 


’ 5C' , 

' 5S ' , 

’ 6C , 

’ 6S ' , 

’ 7C’ , 

’ 7S ' , 

’ 8C' , 

1 8S 


class= 'RESPONSE' , type= ' VAR ' , name= ' INP HUB MOTION AZ HHC INPUT' ScEND 
KINDR=3 , RNAME V = * HUB MOTION AZ HHC INPUT 1 , 

MVAR=17 , KINDN=3 , XNOM=17*0., WTMODE=17 *0 . , XREF=17*0., 

Yaw hub motions are set in here -- MEAN, 1PC, IPS, 2 PC, 

XINIT=0 . , 6*0., 0., 0., 6*0., 0., 0., 

OPSYM=17 * 0 , WTCONV=17 *0 . 05 , WTPERT=17 *0 . 05 , 

TMRED=17*1, TNRED=17*1, FLRED=17*1 


, 2 PS , 


2 PS , . 
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ScNLDEF 

ScNLVAL 


ScEND 

i 

&NLDEF 

ScNLVAL 


ScEND 

i 

ScNLDEF 

ScNLVAL 


& END 

j 

5cNLDEF 

5cNLVAL 


5c END 
! 

ScNLDEF 

ScNLVAL 


5c END 

i 

ScNLDEF 

ScNLVAL 


5c END 

t 

ScNLDEF 

ScNLVAL 


c lass = 'COMPONENT* , type= ' FOURIER ' , name='HUB MOTION OSC AMP AZ 1 ScEND 
LABEL= ' HUB MOTN 1 , MHARM=8 , NBHARM=0, NEHARM= 8 , TYPEA-2 , 

PERIOD= ' ROTOR 1*, AZREF=0 . , GAIN=1.0, OPDRVF=0, OPDRVX^O, 

F ONAME = ’ HUB MOTION AZ HHC INPUT ' , F0LABL= * AZ HHC ' , 

XNAME= ' HUB MOTION OSC AMP AZ 1 , XLABL= 1 AZ HHC' 


class= ' INTERFACE' , type='IO’, name= 1 HUB MOTION OSC AMP AZ ' ScEND 
LABEL= ’ AMP AZ ' , 

TMPART= 'ROTOR 1 CONTROL', TMPASS=1, 

TNPART='NO SOLUTION', 

FLPART= ' NO SOLUTION', 

XNAMEC= ' HUB MOTION OSC AMP AZ ' , 

XNAMEV= ' HUB MOTION OSC AMP AZ ' , 

NCNXN = 1 , 

FNAMEC= ' HUB MOTION', FNAMEV = ' HUB MOTION OSC AMP AZ ’ 


class= ' RESPONSE 1 , type= 'VARIABLE ' , name= ' I NT HUB MOTION OSC AMP AZ ' ScEND 
KINDR=2 , RNAMEV= ' HUB MOTION OSC AMP AZ ' , 

MVAR=1 , KINDN=3 , XNOM=0 . , WTMODE=0 . , XREF=0., XINIT=0., 

OPSYM= 0 , WTCONV=0 . 05 , WTPERT=0.05, 

TMRED-1 , TNRED=1 , FLRED=1 


c 1 ass =' COMPONENT' , type='DE', name =' HUB MOTION' ScEND 
LABEL= ' HUB MOTN ' , OPMTRX = 1 , FORMEQ=0 , INPTEQ=0, INPTDS=1, 
XNAME= 1 HUB MOTION ’ , XLABL= ’ HUB MOTN 1 , LENX=6 , 

NCON=6 , 


CNAME ( 1 ) = 

’HUB 

MOTION 

OSC 

AMP 

LX' , 

CLABL(l) = 

’LX' , 

LENC(l) =1 

CNAME ( 2 ) = 

'HUB 

MOTION 

OSC 

AMP 

LY' , 

CLABL (2 ) = 

’ LY’ , 

LENC ( 2 ) = 1 

CNAME ( 3 ) = 

'HUB 

MOTION 

OSC 

AMP 

LZ' , 

CLABL ( 3 } = 

’LZ' , 

LENC (3 ) =1 

CNAME ( 4 ) = 

' HUB 

MOTION 

OSC 

AMP 

AX' , 

CLABL { 4 ) = 

’ AX’ , 

LENC (4 ) =1 

CNAME ( 5 ) = 

'HUB 

MOTION 

OSC 

AMP 

AY' , 

CLABL ( 5 ) = 

' AY’ , 

LENC (5) =1 

CNAME ( 6 ) = 

' HUB 

MOTION 

OSC 

AMP 

AZ’ , 

CLABL { 6 ) - 

' AZ ' , 

LENC (6) =1 


OPCON^l , LENV=6 , CONV=l , 2 , 3 , 4 , 5 , 6 , CONE=6*l 


class= ' INTERFACE ' , type= ' IO ' , name='HUB MOTION INPUT' ScEND 
LABEL= 'HUB MOTN' , 

TMPART- 1 ROTOR 1 CONTROL’, TMPASS = 2 , 

TNPART- ' NO SOLUTION ' , 

FLPART= ' NO SOLUTION ' , 

XNAMEC= ' HUB MOTION ' , XNAMEV= 1 HUB MOTION ' , 

NCNXN =1 , 

FNAMEC= ' PYLON ' , FNAMEV = ' HUB MOTION INPUT ' 


class= ' RESPONSE ' , type= ' VARIABLE * , name= ’ HUB MOTION INPUT 1 ScEND 
KINDR=2 , RNAMEV= ' HUB MOTION INPUT', 

MVAR=6 , KINDN=3, XNOM=0., WTMODE=0., XREF=0 . , XINIT=0., 

OPSYM=6 *0 , WTCONV=3 *0.05,3*0.02, WTPERT = 3*0 . 05 , , 3*0 . 02 , 

TMRED=6*1 , TNRED=6*1 , FLRED=6*1 


class = ' OUTPUT ’ , type= ' ' , name=’OUT HUB MOTION INPUT' ScEND 
KINDY=4 , YNAMEC= ' HUB MOTION’ , YNAMEV= ’ HUB MOTION ’ , 

TMPART= 'ROTOR 1 HUB LOAD SENSOR', 

TNPART= 'NO SOLUTION’ , 

FLPART= ' NO SOLUTION', 
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NTPRNT=1 , NHPRNT=1, MTIME=24, MHARM=10 , 

&END 


&NLDEF action='end of core' &END 
' endof input ' 


OPHARM-1 
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